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VOLUME  III 

PRINCIPLES  OF  DYNAMIC  ANALYSIS 


INTRODUCTION 


3-1  Purpose 

The  purpose  of  this  six  volume  manual  Is  to  present  methods  of  design  for 
protective  construction  used  In  facilities  for  development,  testing,  produc¬ 
tion,  maintenance,  modification.  Inspection,  disposal  and  storage  of  explosive 
materials. 

3-2  Objectives 

The  primary  objectives  are  to  establish  design  procedures  and  construction 
techniques  whereby  propagation  of  explosion  (from  one  building  or  part  of  a 
building  to  another)  or  mass  detonation  can  be  prevented  and  protection  for 
personnel  and  valuable  equipment  will  be  provided. 

The  secondary  objectives  are: 

(1)  Establish  the  blast  load  parameters  required  for  design  of 
protective  structures; 

(2)  Provide  methods  for  calculating  the  dynamic  response  of  struc¬ 
tural  elements  Including  reinforced  concrete,  structural  steel, 
etc.; 

(3)  Establish  construction  details  and  procedures  necessary  to 
afford  the  required  strength  to  resist  the  applied  blast  loads; 

(4)  Establish  guide  lines  for  siting  explosive  facilities  to  obtain 
maximum  cost  effectiveness  in  both  the  planning  and  structural 
arrangements;  providing  closures,  and  preventing  damage  to 
interior  portions  of  structures  due  to  structual  motion,  shock, 
and  fragment  perforation. 


3-3  Background 

For  the  first  60  years  of  the  20th  century  criteria  and  methods  based  upon  the 
results  of  catastrophic  events  have  been  used  for  the  design  of  explosive 
facilities.  The  criteria  and  methods  did  not  include  a  detailed  or  reliable 
quantitative  basis  for  assessing  the  degree  of  protection  afforded  by  the 
protective  facility.  In  the  late  1960 ' s  quantitative  procedures  were  set 
forth  in  the  first  edition  of  the  present  manual,  "Structures  to  Resist  the 
Effects  of  Accidental  Explosions."  This  manual  was  based  on  extensive 
research  and  development  programs  which  permitted  a  more  reliable  approach  to 
design  requirements.  Since  the  original  publication  of  this  manual,  more 
extensive  publication,  more  extensive  testing  and  development  programs  have 
taken  place.  This  additional  research  was  directed  primarily  towards 
materials  other  than  reinforced  concrete  which  was  the  principal  construction 
material  referenced  in  the  initial  version  of  the  manual. 
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Modern  methods  for  the  manufacture  and  storage  of  explosive  materials,  which 
include  many  exotic  chemicals,  fuels,  propellants,  etc.,  required  less  space 
for  a  given  quantity  of  explosive  material  than  was  previously  needed.  Such 
concentrations  of  explosives  increase  the  possibility  of  the  propagation  of 
accidental  explosions  (one  accidental  explosion  causing  the  detonation  of 
other  explosive  materials).  It  is  evident  that  a  requirement  for  more 
accurate  design  techniques  has  become  essential.  This  manual  describes 
rational  design  methods  to  provide  the  required  stuctural  protection. 

These  design  methods  account  for  the  close-in  effects  of  a  detonation 
including  associated  high  pressures  and  nonuniformity  of  the  blast  loading  on 
protective  structures  or  barriers  as  well  as  intermediate  and  far-range 
effects  which  are  encountered  in  the  design  of  stuctures  which  are  positioned 
away  from  the  explosion.  The  dynamic  response  of  structures,  constructed  of 
various  materials,  or  combination  of  materials,  can  be  calculated,  and  details 
have  been  developed  to  provide  the  properties  necessary  to  supply  the  required 
strength  and  ductility  specified  by  the  design.  Development  of  these 
procedures  has  been  directed  primarily  towards  analyses  of  protective 
structures  subjected  to  the  effects  of  high  explosive  detonation.  However, 
this  approach  is  general  and  is  applicable  to  the  design  of  other  explosive 
environments  as  well  as  other  explosive  materials  as  numerated  above. 

The  design  techniques  set  forth  in  this  manual  are  based  upon  the  results  of 
numerous  full-  and  small-scale  structural  response  and  explosive  effects  tests 
of  various  materials  conducted  in  conjunction  with  the  development  of  this 
manual  and/or  related  projects. 

3-4  Scope  of  Manual 

This  manual  Is  limited  only  by  variety  and  range  of  the  assumed  design 
situation.  An  effort  has  been  made  to  cover  the  more  probable  situations. 
However,  sufficient  general  information  on  protective  design  techniques  has 
been  included  in  order  that  application  of  the  basic  theory  can  be  made  to 
situations  other  than  those  which  were  fully  considered. 

This  manual  is  generally  applicable  to  the  design  of  protective  structures 
subjected  to  the  effects  associated  with  high  explosive  detonations.  For 
these  design  situations,  this  manual  will  generally  apply  for  explosive 
quantities  less  than  25,000  pounds  for  close-in  effects.  However,  this  manual 
is  also  applicable  to  other  situations  such  as  far  or  intermediate  range 
effects.  For  these  latter  cases  the  design  procedures  as  presented  are 
applicable  for  explosive  quantities  up  to  500,000  pounds. 

Because  the  tests  conducted  so  far  in  connection  with  this  manual  have  been 
directed  primarily  towards  the  response  of  structural  steel  and  reinforced 
concrete  elements  to  blast  overpressures,  this  manual  concentrates  on  design 
procedures  and  techniques  for  these  materials.  However,  this  does  not  imply 
that  concrete  and  steel  are  the  only  useful  materials  for  protecitve 
construction.  Tests  to  establish  the  response  of  wood,  brick  blocks, 
plastics,  etc.  as  well  as  the  blast  attenuating  and  mass  effects  of  soil  are 
contemplated.  The  results  of  these  test  may  require,  at  a  later  date, the 
supplementation  of  these  design  methods  for  these  and  other  materials. 
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Other  manuals  are  available  which  enable  one  to  design  protective  structures 
against  the  effects  of  high  explosive  or  nuclear  detonations.  The  procedures 
In  these  manuals  will  quite  often  complement  this  manual  and  should  be 
consulted  for  specific  applications. 

Computer  programs,  which  are  consistent  with  the  procedures  and  techniques 
contained  in  the  manual,  have  been  approved  by  the  appropriate  representative 
of  the  U.S.  Army,  the  U.S.  Navy,  the  U.S.  Air  Force  and  the  Department  of 
Defense  Explosive  Safety  board  (DDESB) .  These  programs  are  available  through 
the  following  repositories: 

1.  Department  of  the  Army 

Commander  and  Director 
U.S.  Army  Engineer 
Waterways  Experiment  Station 
Post  Office  Box  631 
Vicksburg,  Mississippi  39180 

Attn:  WESKA 

2.  Department  of  the  Navy 

Officer-In-Charge 
Civil  Engineering  Laboratory 
Naval  Battalion  Construction  Center 
Port  Hueneme,  California  93043 

Attn:  Code  L51 

. ;  3.  Department  of  the  Air  Force 

Aeorspace  Structures 
Information  and  Analysis  Center 
Wright  Paterson  Air  Force  Base 
Ohio  45433 

Attn:  AFFDL/FBR 


Limited  number  of  copies  of  the  above  program  are  available  for  each 
repository  upon  request.  The  individual  programs  are  identical  at  each 
repository.  If  any  modifications  and/or  additions  to  these  programs  are 
required,  they  will  be  submitted  by  the  organization  for  review  by  DDESB  and 
the  above  services.  Upon  concurrence  of  the  revisions,  the  necessary  changes 
will  be  made  and  notification  of  these  changes  will  be  made  by  the  individual 
repositories. 
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3-5  Foraat  of  Manual 

This  manual  Is  subdivided  Into  six  specific  volumes  dealing  with  various 
aspects  of  design.  The  titles  of  these  volumes  are  as  follows: 


Volume 

I 

-  Introduction 

Volume 

II 

-  Blast,  Fragment  and  Shock  Loads 

Volume 

III 

-  Principles  of  Dynamic  Analysis 

Volume 

IV 

-  Reinforced  Concrete  Design 

Volume 

V 

-  Structural  Steel  Design 

Volume 

VI 

-  Special  Considerations  in  Explosive  Facility  Design 

Appendix  A  pertinent  to  a  particular  volume  and  containing  illustrative 
examples  of  explosive  effects  and  structural  response  problems  appears  at  the 
end  of  each  volume. 

Commonly  accepted  symbols  have  been  used  as  much  as  possible.  However, 
protective  design  involves  many  different  scientific  and  engineering  fields, 
and,  therefore,  no  attempt  has  been  made  to  standardize  completely  all  the 
symbols  used.  Each  symbol  has  been  defined  where  it  is  first  introduced,  and 
a  list  of  the  symbols,  with  their  definitions  and  units,  is  contained  in 
Appendix  B  of  each  volume. 


VOLUME  CONTENTS 


3-6  General 

This  Volume  contains  the  procedures  for  analyzing  structural  elements  subject 
to  blast  overpressures.  These  procedures  are  contained  in  the  next  eleven 
sections;  Section  3-7  deals  with  a  simplified  discussion  of  the  basic 
principles  of  dynamics  as  well  as  the  procedures  for  calculating  the  various 
components  used  to  perform  the  dynamic  analyses.  Presented  in  Sections  3-8 
through  3-15  are  resistance-deflection  functions  for  various  elements 
including  both  one-  and  two-way  panels  as  well  as  beam  elements.  These 
functions  include  the  elastic,  el  asto-plastic,  and  plastic  ranges  of 
response.  In  addition,  a  discussion  of  dynamic  equivalent  systems  is 
presented  in  Sections  3-16  and  3-17.  These  include  single-  and  multi -degree- 
of-freedom  systems.  Presented  in  this  Section  also  are  methods  for 
calculating  load  and  mass  factors  required  to  perform  the  dynamic  analyses. 

Sections  3-18  through  3-?0  include  both  a  step-by-step  numerical  integration 
of  an  element's  motion  under  dynamic  loads  utilizing  the  Acceleration-Impulse- 
Extrapolation  Method  or  the  Average  Acceleration  Method  and  design  charts  for 
idealized  loads.  Presented  also  in  these  Sections  are  methods  for  analyzing 
elements  subjected  to  impulse  type  loadings;  that  is,  loadings  whose  durations 
are  short  in  comparison  to  the  time  to  reach  maximum  response  of  the  elements. 
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BASIC  PRINCIPLES 


3-7  General 

The  principles  used  In  the  analysis  of  structures  under  static  load  will  be 
reviewed  briefly,  since  the  same  principles  are  used  In  the  analysis  and 
design  of  structures  subjected  to  dynamic  loads.  Two  different  methods  are 
used  either  separately  or  concurrently  In  static  analysis:  one  is  based  on 
the  principle  of  equilibrium,  and  the  other  on  work  done  and  internal  energy 
stored. 

Under  the  application  of  external  loads,  a  given  structure  is  deformed  and 
internal  forces  developed  in  its  members.  In  order  to  satisfy  static 

equilibrium,  the  vecto'  sum  of  all  the  external  and  internal  forces  acting  on 

any  free  body  portion  of  the  stucture  must  be  equal  to  zero.  For  the 

equilibrium  of  the  structure  as  a  whole,  the  vector  sum  of  the  external  forces 
and  the  reactions  of  the  foundation  must  be  equal  to  zero. 

The  method  based  on  work  done  and  energy  considerations  is  sometimes  used  when 
It  is  necessary  to  determine  the  deformation  of  a  structure.  In  this  method, 
use  is  made  of  the  fact  that  the  deformation  of  the  structure  causes  the  point 
of  application  of  the  external  load  to  be  displaced.  The  force  then  does  work 
on  the  structure.  Meanwhile,  because  of  the  structural  deformations, 
potential  energy  is  stored  in  the  structure  in  the  form  of  strain  energy.  By 

the  principle  of  energy  conservation,  the  work  done  by  the  external  force  and 

the  energy  stored  in  the  members  must  be  equal.  In  static  analysis, 

simplified  methods  such  as  the  method  of  virtual  work  and  the  method  of  the 

unit  load  are  derived  from  the  general  principle  of  energy  conservation. 

In  the  analysis  of  statically  indeterminate  structures,  in  addition  to 
satisfying  the  equations  of  equilibrium,  it  is  necessary  to  include  a 
calculation  of  the  deformation  of  the  strucuture  in  order  to  arrive  at  a 

complete  solution  of  the  Internal  forces  in  the  structure.  The  methods  based 
on  energy  considerations  such  as  the  method  of  least  work  and  the  method  based 
on  Castigliano's  theorems  are  generally  used. 

For  the  analysis  of  structures  under  dynamic  loading,  the  same  two  methods  are 
basically  used;  but  the  load  changes  rapidly  with  time  and  the  acceleration 
velocity  and,  hence,  the  inertia  force  and  kinetic  energy  are  of  magnitudes 
requiring  consideration.  Thus,  in  addition  to  the  internal  and  external 
forces,  the  equation  of  equilibrium  includes  the  inertia  force  and  the 
equation  of  dynamic  equilibrium  takes  the  form  of  Newton's  equation  of  motion: 

F  -  R  =  Ma  3-1 

where  F  =  total  external  force  as  a  function  of  time 
R  =  total  internal  force  as  a  function  of  time 
M  =  total  mass 

a  =  acceleration  of  the  mass 

As  for  the  principle  of  conservation  of  energy,  the  work  done  must  be  equal  to 
the  sum  of  the  kinetic  energy  and  the  strain  energy: 
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WD  =  KE  +  SE 


3-2 


where  WD  =  work  done 

KE  *  kinetic  energy 
SE  =  strain  energy 

and  the  strain  energy  Includes  both  reversible  elastic  strain  energy  and  the 
Irreversible  plastic  strain  energy.  Thus,  the  difference  between  structures 
under  static  and  dynamic  loads  Is  the  presence  of  inertial  force  (Ma)  in  the 
equation  of  dynamic  equilibrium,  and  of  kinetic  energy  In  the  equation  of 
energy  conservation.  Both  terms  are  related  to  the  mass  of  the  structure; 
hence,  the  mass  of  the  structure  becomes  an  important  consideration  in  dynamic 
analysis. 

In  the  dynamic  analysis  of  structures,  both  the  energy  balance  equation  and 
the  force  balance  equation  are  applied  with  explicit  description  of  the 
external  forcing  function  F,  and  the  internal  resisting  forcing  function  R. 
The  difference  between  these  forcing  functions  is  the  inertia  force  as 
described  above.  The  following  Is  a  discussion  of  the  details  of  how  these 
forces  are  utilized  In  the  design  of  structures  which  respond  in  the  ductile 
mode. 

In  the  design  of  a  structure  to  resist  the  blast  from  an  HE  explosion,  the 
total  external  force  acting  on  the  structure  can  be  obtained  by  the  principles 
discussed  Volume  II.  The  design  method  also  consists  of  the  determination  of 
the  total  internal  force,  l.e.  the  resistance  of  the  structure  required  to 
limit  calculated  deflections  of  the  individual  members  and  the  structure  as  a 
whole  under  the  external  force  (blast  loading),  to  within  prescribed  maximum 
values.  The  determination  of  the  resistance  of  the  individual  members  of  the 
structure  is  presented  in  Sections  3-8  through  3-17.  Subsequent  sections  of 
this  manual  present  the  principles  and  methods  of  dynamic  analysis  and 
equations,  charts,  and  procedures  for  design. 


RESISTANCE  -  DEFLECTION  FUNCTIONS 

3-8  Introduction 

Under  the  action  of  external  loads,  a  structural  element  is  deformed  and 
internal  forces  set  up.  The  sum  of  these  internal  forces  tending  to  restore 
the  element  to  its  unloaded  static  position  is  defined  as  the  resistance.  The 
resistance  of  a  structural  element  is  a  reactive  force  associated  with  the 
deflection  of  the  element  produced  by  the  applied  load.  It  is  convenient  to 
consider  the  resistance  as  an  equivalent  load  in  the  same  manner  as  the 
applied  load,  but  opposite  in  direction.  The  variation  of  the  resistance  vs. 
displacement  is  expressed  by  a  resistance-deflection  function  and  may  be 
represented  graphically.  An  idealized  resistance-deflection  function  for  an 
element  spanning  in  two  directions  and  covering  in  the  complete  flexural  range 
to  incipient  failure  is  shown  in  figure  3-1. 

As  load  is  applied  to  a  structural  element,  the  element  deflects  and,  at  any 
instant,  exerts  a  resistance  to  further  deformation,  which  is  a  function  of 
its  units  stiffness  K,  until  the  ultimate  unit  resistance  ru  (total  resistance 
is  ry A  where  A  is  the  element  area)  of  the  element  is  reached  at  deflection 
Xp.  The  intial  portion  of  the  resistance-deflection  diagram  is  composed  of 


ELASTO- PLASTIC  RANGE 


the  elastic  and  elasto-plastlc  ranges,  each  with  its  corresponding  stiffness, 
the  transition  from  one  range  to  another  occurring  as  plastic  hinges  are 
formed  at  points  of  maximum  stress  (yield  lines).  The  number  of  elasto- 
plastic  ranges  required  before  the  ultimate  resistance  of  a  particular  element 
is  reached  depends  upon  the  type  and  number  of  supports,  and  the  placement  of 
reinforcing  steel  (in  the  case  of  reinforced  concrete  elements).  For  example 
a  beam  with  simple  supports  subjected  to  uniformly  distributed  loads  needs 
only  one  plastic  hinge  to  develop  the  ultimate  resistance  (or  full  plastic 
strength)  of  the  element;  whereas  for  the  same  beam  fixed  at  both  ends,  more 
than  one  plastic  hinge  is  required. 

In  subsequent  paragraphs,  various  procedures,  equations  and  illustrations  are 
presented  to  enable  the  designer  to  determine  the  resistances  of  both  one-  and 
two-way  elements.  The  piocedures  outlined  apply  mainly  to  reinforced  concrete 
elements  and  so  do  the  equations  appearing  in  the  text,  unless  the  equations 
are  given  as  part  of  an  illustrative  example.  However,  they  can  also  be  used 
for  structural  steel  elements  as  well  as  other  structural  elements  such  as 
aluminum,  plastics,  etc.  Equations  have  been  derived  for  specific  cases  most 
often  encountered  in  practice.  These  are  applicable  for  structural  steel  and 
reinforced  concrete  elements  of  uniform  thickness  in  both  the  horizontal  and 
vertical  directions.  Before  the  equations  and  figures  can  be  used  for 
reinforced  concrete  element,  however,  the  reinforcing  steel  across  any  yield 
line  must  have  a  uniform  distribution  in  both  the  vertical  and  horizontal 
directions;  however,  the  reinforcement  across  the  positive  yield  lines  can  be 
different  from  that  across  the  negative  yield  lines  and  the  reinforcing 
pattern  in  the  vertical  direction  different  from  that  in  the  horizontal 
direction. 

Regardless  of  whether  it  is  reinforced  concrete  or  structural  steel  element, 
any  opening  in  the  element  must  be  compact  in  shape  and  small  in  area, 
compared  to  the  total  area  of  the  element. 

3-9  Ultimate  Resistance 

3-9.1  Gene-pal 

The  ultfmate  resistance  of  an  element  depends  upon: 

(1)  The  distribution  of  the  applied  loads. 

(2)  The  geometry  of  the  element  (length  and  width). 

(3)  The  number  and  type  of  supports. 

(4)  The  distribution  of  the  moment  capacity  or  reinforcement  in  the 
case  of  reinforced  concrete  elements. 

The  distribution  of  the  loads  depends  upon  the  design  range  of  the  element; 
i.e.,  high,  intermediate  or  low  pressure.  For  intermediate  and  low  pressure 
ranges,  it  can  be  assumed  that  the  pressure  is  uniform  across  the  surface  of 
the  element  although  it  varies  with  time.  At  high  pressure  ranges,  however, 
the  blast  loads  are  variable  across  the  surface  of  the  element.  However,  for 
structural  steel  elements  and  concrete  elements  utilizing  laced  rei nforcement, 
or  for  concrete  elements  with  standard  shear  reinforcement  which  sustain 
relatively  small  deflections,  a  good  estimate  of  the  resulting  deflections  can 
be  made  using  the  resistance  functions  conforming  to  those  of  uniformly  loaded 
elements. 


The  other  factors  that  affect  the  ultimate  resistance  of  an  element  are 
predetermined  by  the  requirements  of  the  protective  structure  (where  the 
element  is  used)  and  the  magnitude  of  the  blast  output. 

3-9.2  One-Way  Elements 

The  ultimate  resistance  of  a  one-way  reinforced  concrete  element  with  an 
elastic  distribution  of  Its  reinforcing  steel  is  based  on  the  moment  capacity 
at  first  yield  since  all  critical  sections  yield  simultaneously.  For  one-way 
reinforced  concrete  elements  (such  as  beams  or  slabs)  with  non-elastic 
distribution  of  reinforcing  steel  and  for  structural  steel  elements,  the 
ultimate  resistance  is  a  function  of  the  moment  capacity  at  the  first  yield 
plus  the  added  moment  capacity  due  to  subsequent  yielding  at  other  critical 
sections. 

Values  of  the  ultimate  resistance  for  one-way  elements  are  shown  in  table  3-1 
where  the  following  symbols  are  used: 

Mn  =  ultimate  negative  unit  moment  capacity  at  the  support. 

M_  -  ultimate  positive  unit  moment  capacity  at  midspan. 

L  =  length 

ru  =  ultimate  unit  resistance 
Ru  =  total  ultimate  resistance 

Table  3-1  applies  to  both  beams  and  slabs.  However,  special  attention  must  be 
paid  to  the  units  used  for  the  respective  element.  The  moment  capacity  of  a 
slab  is  expressed  for  a  unit  strip  of  the  slab  (inch-pounds  per  inch)  whereas 
the  total  moment  capacity  (inch-pounds)  is  considered  for  a  beam. 
Consequently,  the  resistance  of  a  slab  is  expressed  in  load  per  unit  area 
(psi)  where  the  resistance  of  a  beam  is  expressed  in  load  per  length  along  the 
beafn  (pounds  per  inch). 

3_9.3  Two-Way  Elements 

The  amount  of  data  available  on  the  limit  analysis  of  retangular  steel  plates 
is  very  limited.  However,  an  elementary  approach  imagines  a  mechanism  formed 
of  straight  yield  lines,  as  is  customary  in  reinforced  concrete.  This 
approach  for  reinforced  concrete  elements  will  be  considered  appropriate  for 
structural  steel  elements. 

In  the  design  of  two-way  reinforced  concrete  elements,  it  is  not  necessary  to 
define  accurately  the  stress  distribution  during  the  initial  and  intermediate 
stages  of  loading  since  the  ultimate  load  capacity  can  be  readily  determined 
by  the  use  of  yield  line  procedures.  The  yield  line  method  assumes  that  after 
inital  cracking  of  the  concrete  at  points  of. maximum  moment,  yielding  spreads 
until  the  full  moment  capacity  is  developed  along  the  length  of  the  cracks  on 
which  failure  will  take  place.  Several  illustrative  examples  of  the 
simplified  yield  or  crack  lines  for  two-way  elements  are  illustrated  in  figure 
3-2. 

In  using  the  yield  line  solution,  the  inital  step  is  to  assume  a  yield  line 
pattern  (as  shown  in  figure  3-2)  applying  the  following  rules: 

(1)  To  act  as  plastic  hinges  of  a  collapse  mechanism  made  up  of 


plane  segments,  yield  lines  must  be  straight  lines  forming  axes 
of  rotation  for  the  movements  of  the  segments. 


(2)  The  supports  of  the  slabs  will  act  as  axes  of  rotation.  A  yield 
line  may  form  along  a  fixed  support  and  an  axis  of  rotation  will 
pass  over  a  column. 

(3)  For  compatibility  of  deformations,  a  yield  line  must  pass 
through  the  intersection  of  the  axes  of  rotation  of  the  adjacent 
slab  segments. 


Tests  indicate  that  the  actual  location  and  extent  of  these  lines  on 
reinforced  concrete  elements  differ  only  slightly  at  failure  from  the 
theoretical  ones.  Use  of  the  idealized  yield  lines  results  in  little  error  in 
the  determination  of  the  ultimate  resistance  and  the  error  is  on  the  side  of 
safety. 

The  corner  sections  of  two-way  elements  are  stiff  in  comparison  to  the 
remainder  of  the  member;  therefore,  straining  of  the  reinforcement  which  is 
associated  with  the  reduced  rotations  at  these  sections  will  be  less.  To 
account  for  the  corner  effects,  the  design  of  any  one  particular  section  of  a 
two-way  element  should  consider  a  variation  of  the  moment  capacity  along  the 
yield  lines  rather  than  a  uniform  distribution. 

This  variation  is  approximated  by  taking  the  full  moment  capacity  along  the 
yield  lines,  except  in  the  corners  where  two-thirds  of  the  moment  capacity 
over  the  lengths  described  in  figure  3-3  are  used.  The  variation  applies  to 
both  the  negative  moments  along  the  supports  and  the  positive  moments  at  the 
interior. 

The  ultimate  unit  resistance  can  be  determined  from  the  yield  line  pattern 
using  either  the  principle  of  virtual  work  or  the  equations  of  equilibrium. 
Each  approach  has  its  advantages;  in  general,  the  virtual  work  method  is 
easier  in  principle  but  difficult  to  manipulate  algebraically  since  it 
involves  differentiating  a  usually  complex  mathematical  expression  for  a 
minimum  value  of  resistance.  The  equilibrium  method,  which  is  used  in  this 
manual,  also  has  its  disadvantages.  Since  equilibrium  requires  that  the  shear 
forces  acting  on  each  side  of  a  yield  line  have  to  be  equal  and  opposite, 
correction  forces  (also  known  as  nodal  forces)  have  to  be  introduced  around 
openings  in  two-way  members  and  at  free  edges,  and  these  correction  forces  may 
not  be  available  from  simple  analysis.  However,  in  three  of  the  six  cases 
shown  in  figure  3-2,  (cases  c,  e,  f),  nodal  forces  exist;  but  their  effects 
are  negligible. 

In  order  to  calculate  the  ultimate  unit  resistance  r  of  a  two-way  element, 
the  equation  of  equilibrium  of  each  sector  formed  ny  the  yield  lines  is 
expressed  in  terms  of  the  moments  produced  by  the  internal  and  external 
forces.  The  sum  of  the  resisting  moments  acting  along  the  yield  lines  (both 
positive  and  negative)  of  each  sector  is  equated  to  the  moment  produced  by  the 
applied  load  about  the  axis  of  rotation  (supoort  of  the  sector),  assuming  that 
the  shear  forces  are  zero  alonq  the  positive  yield  lines. 

2Mn  +  2Mp  =  Rc  =  ruAc  3-3 
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Figure  3-3  Determination  of  ultimate  unit  resistance 
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where 


Mn  =  sum  of  the  ultimate  unit  resisting  moments 

acting  along  the  support  (negative  yield  lines) 

M_  =  sum  of  the  ultimate  unit  resisting  moments 
p  acting  along  the  interior  failure  lines  (positive 
yield  lines) 

R  =  total  ultimate  resistance  of  the  sector 
c  *  distance  from  the  centroid  of  the  load  to 
the  line  of  rotation  of  the  sector 
r.  =  ultimate  unit  resistance  of  the  sector 
A  =  area  of  the  sector 

Once  the  equations  of  equilibrium  are  known  for  all  sectors,  the  ultimate 
resistance  is  obtained  either  by  solving  the  equations  simultaneously  or  by  a 
trial  and  error  procedure  noting  that  the  unit  resistance  of  all  sectors  must 
be  equal. 

To  illustrate  the  above  procedure  (equation  3-3),  consider  the  two-way 
concrete  element  shown  in  figure  3-3  which  is  fixed  on  three  edges  and  free  on 
the  fourth,  and  where  the  nomenclature  is  as  follows: 

L  =  length  of  element 

H  =  height  of  element 

x  =  yield  line  location  in  horizontal  direction 

y  =  yield  line  location  in  vertical  direction 

MyN  =  ultimate  unit  negative  moment  capacity  in 
the  vertical  direction 

My p  =  ultimate  unit  positive  moment  capacity 
in  the  vertical  direction 
M^j  =  ultimate  unit  negative  moment  capacity 
in  the  horizontal  direction 
MHP  =  ultimate  unit  positive  moment  capacity 
in  the  horizontal  direction 

The  nomenclature  as  stated  in  the  paragraph  above  is  strictly  applicable  to 
two-way  elements  which  are  used  as  walls.  However,  when  roof  slabs  or  other 
horizontal  elements  are  under  consideration,  the  preceding  nomenclature  will 
also  be  applicable  if  the  element  is  treated  as  being  rotated  into  a  vertical 
position. 

The  first  step  in  the  solution  is  to  assume  the  location  of  the  yield  lines  as 
defined  by  the  coordinates  x  and  y.  It  should  be  noted  that  in  some  cases, 
because  of  geometry,  the  value  of  x  and  y  will  be  known  and  therefore  need  not 
be  evaluated.  In  this  example,  the  negative  reinforcement  in  the  horizontal 
direction  at  opposite  supports  is  assumed  to  be  equal;  therefore,  the  vertical 
yield  line  is  located  at  the  center  of  the  span  and  the  value  of  x  is 
numerically  equal  to  L/2  (a,  fig.  3-3).  However,  in  other  cases,  neither  the 
location  of  x  nor  y  will  be  known,  and  the  solution  will  require  the 
determination  of  both  coordinates. 

Once  the  yield  lines  have  been  assumed,  the  distribution  of  the  resisting 
moments  along  the  yield  lines  is  determined.  In  the  case  at  hand,  the  reduced 
moments,  as  a  result  of  the  increased  stiffness  at  the  corners,  act  over 
lengths  equal  to  x/2  and  y/2  in  the  horizontal  and  vertical  directions, 
respectively  (a,  fig.  3-3).  The  equations  of  equilibrium  are  then  written  for 


each  sector  with  the  use  of  the  free  body  diagrams  (b,  fig.  3-3).  For  the 
triangular  sector  I: 


Mvn  =  (2/3)MVN(L/4  +  L/4)  +  Myjgd/2) 

=  (5/6)MvNL 

3-4 

Myp  =  ( 2/3 )Myp (L/4  +  L/4)  +  Mvp(L/2) 

=  (5/6)MypL 

3-5 

Cj  ■  y/3 

3-6 

Rj  =  (  Myj^  +  Myp)/Cj 
=  [5L(Mvn  +  Myp)]/2y 

3-7 

Aj  *  Ly/2 

3-8 

ru( Sector  I)  =  Rj/Aj 

=  [5(Mvn  +  Mvp)]/y2 

3-9 

For  the  trapezoidal  sector  II,  a  similar  procedure  gives 

Mhn  =  ( 2/3 )MHN (y/2 )  +  Mhn(H  -  y/2) 

=  Mhn(H  -  y/6) 

3-10 

mHP  =  (2/3)MHP(y/2)  +  MHP(H  -  y/2) 

=  MHp(H  -  y/6) 

3-11 

Cn  =  (1/3) (L/2)[2(H-y)  +  H]/(H  +  H  -  y) 

- 

=  [L(3H  -  2y)]/6(2H  -  y) 

3-12 

RII  =  {  MHN  +  MHP)/CII 

=  [(6H  -  y)(2H  -  y)(MHN  +  MHp)]/L(3H  -  2y) 

3-13 

An  =  ( 1/2 )  ( L/2 )  ( H  +  H  -  y) 

=  [L(2H  -  y) ]/4 

3-14 

ru(Sector  II)  =  Rji/Ajj 

=  C4(Mhn  +  MHp)(6H  -y)]/L2(3H  -  2y) 

3-15 

Equations  3-9  and  3-15  are  the  equations  of  equilibrium  for  the  triangular  (I) 
and  the  trapezoidal  (II)  sectors,  respectively.  As  mentioned  previously, 
these  equations  can  be  solved  simultaneously  or  by  a  trial  and  error 
procedure.  In  the  latter  method,  values  of  y  are  substituted  into  both 
equations  until  ru  (sector  I)  is  equal  to  ru  (sector  II). 
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If  a  numerical  solution  based  on  the  above  procedure  (equation  3-3)  yields 
negative  values  for  either  x,  y  or  r,,  then  the  assumed  yield  line  location  is 
wrong.  In  this  example,  the  only  otner  possible  yield  line  pattern  (x^L/2) 
would  be  as  shown  In  figure  3-2c. 

The  solution  of  equation  3-3  Is  universally  applicable  for  any  two-way 
element.  If  the  negative  reinforcement  In  the  horizontal  direction  had  been 
unequal  at  the  opposing  supports,  the  value  of  x  =  L/2  would  have  changed,  and 
all  three  sectors  wold  have  had  to  be  considered  to  determine  x,  y  and  hence, 


Simultaneous  solution  of  equations  3-9  and  3-15  reveals  that  the  locations  of 
the  yield  lines  are  a  function  of  the  ratio  of  the  spans  L/H  and  the  ratio  of 
the  sum  of  the  unit  vertical  to  horizontal  moment  capacities  as  follows: 

ru( Sector  I)  =  ru( Sector  II)  3-16 

5(Mvn  +  MypJ/y2  =  C4(Mhn  +  MHp)(6H  -  -  2y) 


l2(mvn  +  mvp)/h2(mhn  +  mhp)  = 

C 4y 2 ( 6  -  y/H)]/[5H2(3  -  2y/H)3 

( L/H )  [  ( Myjyj  +  Myp)/(Mjjm  +  M^p)]^2  = 

(y/H)[(4(6  -  y/H)/5(3  -  2y/H)]1/2 


3- 1 7a 


3-17b 


3-17c 


Equation  3-17c,  which  relates  the  location  of  the  yield  lines  to  the  moment 
capacity  of  the  element.  Is  used  to  plot  figure  3-6.  Knowing  the  location  of 
the'  yield  lines,  the  resistance  of  the  two-way  element  can  be  obtained  from 
either  equation  3-9  or  3-15  which  are  also  presented  in  table  3-2. 

Using  the  procedure  outlined  above,  the  values  of  the  ultimate  unit 
resistances  for  several  two-way  elements  with  various  support  conditions  are 
given  in  tables  3-2  and  3-3,  the  nomenclature  confirming  to  that  previously 
listed.  Table  3-2  covers  the  special  cases  where  opposite  supports  provide 
the  same  degree  of  restraint  thus  resulting  in  symmetrical  yield  line 
patterns.  Table  3-3  deals  with  the  general  cases  when  the  yield  line  patterns 
are  not  symmetrical  (that  is,  when  opposite  supports  provide  different 
restraints).  Yield  line  location  ratios  x/L  and  y/H  for  the  same  elements  are 
depicted  in  figures  3-4  through  3-20. 

Figures  3-4  and  3-5  show  the  location  of  the  yield  lines  for  two-way  elements 
with  two  adjacent  edqes  supported  and  the  other  two  free.  In  each  of  these 
figures  eight  curves  are  shown  which  represent  different  ratios  of  the 
positive  to  the  negative  moment  capacities  in  both  the  vertical  and  horizontal 
directions.  Figures  3-6  through  3-16  illustrate  the  yield  line  location  for 
two-way  elements  with  three  edges  supported  and  one  edge  free.  Figures  3-6 
and  3-11  covers  the  case  when  the  yield  line  pattern  is  symmetrical  (opposite 
supports  provide  the  same  degree  of  restraint).  Figures  3-17  through  3-20 
show  the  yield  line  location  for  two-way  elements  with  four  sides  supported. 


.  %  .  ’  «  ’  ,  ^  , 


.  v  *  «  '  •  * . 
.  .  .  ■  -  '  > 


.NV-'.V.' 


.'•V  ' 


™  : 


of  yield  lines  for  two-way  element  with  two  adjacent 
pported  and  two  edges  free  (values  of  x) 


n  of  yield  lines  for  two-way  element  with  two  adjacent 
upported  and  two  edges  free  (values  of  y) 


>  > 

2  2 


ztm 

m 

him 

'  "  7  f'  -  M 

■Ipr  miiieml 
■aum 
/j 

19 

■  ■  ■  WMMMaman 

w far  MX  **i 

'jf/ffA'saVi 

fi 

msm 

,  wB 

wm ■ 


■■■w^n  HiHiiiimminmimnnmirai 
lilMlSHusMiiiiiuBSafli'ia  aiERliiisgamai^mi 


cai  y  i  el 


unsymmetrical  yield  lines  for  two-way  element 
dges  supported  and  one  edge  free  (X2/Xj=0.5) 


n  of  unsymmetrical  yield  lines  for  two-way  element 
ree  edges  supported  and  one  edge  free  (X2/Xi=0.75) 


el ement 


cal 


iiiiip2iiii=i=m=3na 


=====“:  ; 


,!|frr 


z—MM:;:; 


IMpjBj 


wn§ 


:=Hr//; 


7/i 


nr«f 


Yj/MS&iiM 
W/SMiMMm 
W/f/iMifMm 

ji/xSss^^a 


mf/iSr 


»  ~ 


'/J P-^rjS 

HOOP 


is-s-si 

~~  — >  — 

/  *  MV  »*WMB 

'tftariavM  _ 

fiwrju=2\ 


aBigisagaBj 

— llllBMBnMiMiln 


laiiiHiaiifiiiij 


N 

Q.  Z 
>  > 
2  2 


>  + 

2  z 
z 

>  2  i 


nniHminwffTgiii«j‘  ■■i  iwii»»u)»n» 

wihRU ^  ....  .  ...  ..  ... .  ... 

B8HII  nBRBBB&sssassBa 

jcp  |Sp;Eli  Si  SS  SfiJS 


Location  of  unsymmetrical  yield  lines  for  two-way  element  with 
three  edges  supported  and  one  edge  free  (X2/X1  =  1.75) 
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Location  of  symmetrical  yield  lines  for  two-way  element  with  four  edges 


nsymmetrical  yield  lines  for  two-way  element  with 
pported  (values  of  X]) 
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Figure  3-19  Location  of  unsymmetrical  yield  lines  for  two-way  elemen 
four  edges  supported  (values  of  Yj) 
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Figure  3-20  Location  of  unsymmetri cal  yield  lines  for  two-way  element  with  four  edges 
supported  (values  of  x/l  and  •VH) 


Figure  3-17  covers  the  special  case  when  opposite  supports  provide  the  same 
degree  of  restraint  thus  resulting  In  a  symmetrical  yield  pattern.  An  example 
Illustrating  the  use  of  some  of  these  figures  Is  provided  in  Appendix  A. 

3-9.4  Openings  in  Two-Way  .Elements 

The  use  of  openings  in  two-way  elements,  whether  for  access  as  a  door  opening 
or  for  visual  communication  as  in  the  case  of  observation  ports,  is 
permissible  with  certain  reservations.  It  is  difficult  to  state  exact  rules 
concerning  openings,  but  their  effect  on  the  design  is  generally  a  function  of 
location,  size  and  shape. 

Small  compact  openings  with  approximate  areas  of  less  than  5  percent  of  the 
panel  area  and  located  away  from  regions  of  high  stress  can  usually  be  ignored 
In  the  design.  However,  as  in  the  case  of  conventional  design,  reinforcement 
at  least  equal  to  the  amount  interrupted  should  be  placed  adjacent  to  the 
opening.  For  example,  in  figure  3-21,  the  openings  r.ho*  n  (a)  and  (b)  can 
be  disregarded.  If  the  opening  in  (b)  were  made  more  rxt  rqular  as  in  (c), 
then  the  design  must  be  modified  to  account  for  the  cha<  t  ii  the  y;eld  lines 
and,  hence,  the  change  in  the  resistance.  This  char;;'  in  resistance  is  a 
function  of  both  the  shape  and  the  location  of  the  opening. 

Door  openings  invariably  require  special  analysis  because  of  t!  .•  r  size.  As 
depicted  in  (d),  (e)  and  (f),  figure  3-21,  the  presence  of  door  openings 
causes  gross  relocations  of  the  yield  lines  which  generally  propagate  from  the 

corners  of  the  openings.  Since  the  door  also  sustains  the  blast  loading, 

concentrated  line  loads  are  present  around  the  periphery  of  such  openings. 

These  concentrated  loads  must  be  included  in  the  analysis  since  they  change 
the  resistance.  As  previously  outlined  for  solid  elements  and  for  this  case 
also,  the  yield  line  locations  are  assumed  and  each  sector  is  individually 
analyzed.  The  presence  of  line  loads  modifies  equation  3-3  to 

+  ZMp  =  ruAc  +  Mc  3-18 

where  Mc  is  the  moment  of  the  concentrated  loads  about  the  line  of  rotation  of 
the  sector  being  considered.  Solution  of  elements  with  openinns  is  most 
easily  accomplished  through  a  trial  and  error  procedure  by  setting  up  the 

simultaneous  equations  for  each  sector  and  assuming  various  values  of  x  and  y 
until  the  several  values  of  ru  agree  to  within  a  few  percent. 

3-10  Post-Ultimate  Resistance 

In  general,  the  two-way  elements  described  in  this  manual  exhibit  a  post- 
ultimate  resistance  after  initial  failure  occurs  as  indicated  in  figure  3-1. 
Prior  to  this  partial  failure,  the  element  is  spanning  in  two  directions  with 
a  resistance  equal  to  the  ultimate  resistance  ru.  At  a  particular  deflection, 
denoted  as  X^,  failure  occurs  along  one  side  or  two  opposite  sides,  and  the 
element  then  spans  in  one  direction  with  the  reduced  post-ultimate  unit 
resistance  rUp  until  complete  failure  occurs  at  deflection  Xu.  One-way 
elements  do  not  exhibit  this  behavior. 

The  location  of  the  yield  lines  determines  the  presence  or  absence  of  this 
range.  If  the  yield  lines  emanating  from  the  corners  of  the  elements  bisect 
the  90-deqree  corner  angle,  then  all  supports  fail  simultaneously  and  there  is 


no  post-ultimate  range.  As  previously  shown,  the  location  of  the  yield  lines 
for  a  particular  element  is  a  function  of  L/H  and  the  ratio  of  the  unit 
vertical  to  horizontal  moment  capacities.  Post-ultimate  resistances  for  two- 
way  elements  are  shown  in  table  3-4. 

3-11  Partial  Failure  and  Ultimate  Deflection 

Partial  failure  deflection  Xj,  for  two-way  elements  and  ultimate  deflections 
Xu  for  both  one-way  and  two-way  elements  are  a  function  of  the  angle  of 
rotation  of  the  element  at  its  supports  and  the  geometry  of  the  sectors  formed 
by  the  position  yield  lines. 

Once  the  ultimate  resistance  r,.  is  reached  (full  moment  capacity  developed 
along  the  yield  lines),  the  structural  element  becomes  a  mechanism  which 
rotates  with  no  further  increase  in  either  the  moment  or  curvature  between  the 
hinges.  For  one-way  elements,  the  rotation  continues  and  the  deflection 
increases  until  either  the  maximum  deflection  Xm  is  reached  or  failure  occurs 
at  0  max.  The  equations  for  the  maximum  deflection  Xm  in  the  range  0  <  X_  < 
X  for  several  one-way  elements  as  a  function  of  the  rotation  angle  0  arid  tne 
ultimate  deflection  Xu  are  given  in  table  3-5,  when  the  values  for  Xu  are 
based  on  the  development  of  a  maximum  support  rotation,  0  max.,  prior  to 
failure. 

Actually,  tie  maximum  support  rotation  will  vary  with  the  material  type  and 
geometry  of  the  element.  The  criteria  for  partial  and  incipient  failure  for 
concrete  and  structural  steel  elements  can  be  found  in  Volumes  IV  and  V 
respectively. 

For  two-way  elements,  the  rotations  of  all  the  sectors  must  be  considered  in 
order  to  define  the  deflections  of  partial  and  incipient  failure.  Prior  to 
partial  failure  (0  <  Xj ) ,  the  maximum  deflection  is  a  function  of  the 

larger  angle  of  roTation  formed  along  either  the  vertical  or  horizontal 
supports.  At  deflection  Xj,  this  larger  angle  equals  0  max,  and  failure 
occurs  along  this  support.  Beyond  this  point,  the  element  spans  in  one 
direction  until  the  angle  of  rotation  at  the  adjacent  supports  (in  the 
direction  opposite  to  that  at  which  failure  has  already  occured)  reaches  Q  at 
which  time  total  collapse  occurs  (Xm  =  Xu). 

To  illustrate  the  above,  consider  a  two-way  element  (fig.  3-22)  which  is  fully 
restrained  on  four  edges  and  whose  positive  yield  lines  are  defined  by  H/2  <_x 
<  L/2  and  y  =  H/2.  Denoting  0u  as  the  angle  of  rotation  in  the  horizontal 
direction  (along  vertical  supports)  and  0y  as  the  vertical  angle  of  rotation 
(along  horizontal  supports),  the  maximum  deflection  Xm  at  the  center  of  the 
element  prior  to  reaching  the  deflection  is 

Xm  =  (H  tan  0v)/2  3-19 

and  at  the  partial  failure  deflection  X^  where  0y  max 

Xj  =  (H  tan  0max)/2  3-20 
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Table  3-4  Post-Ultimate  Unit  Resistances  for  Two-Way  Elements 


y  >  L/2  8  (  Mhn  *  Mhp)  /  H 
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Table  3-5  General  and  Ultimate  Deflections  for  One-Way  Elements 


Edg<  Conditions  and  Loading  Diagrams  0eJfS7Xm  D.flecTior'x, 


Ultimate 


-  tan  9  max 


j  tan  0 


(L/2)tan0  (L/2)tan  0  max 


Y  tan  0 


(L/2)  tan  0  (L/2)tan  0  max 


L/3  i  L/3 


-  tan  0  max 


(L/2)  tan  0 

(L/2)tan  0  max 

L  tan  0 

L  tan  0  max 

L  tan  9 

L  tan  0  max 

(L/3)  tan  0 

(L/3)  tan  0  max 

ELEVATION 


y  =  H/z 


—  PARTIAL  FAILURE 

—  INCIPIENT  FAILURE 


Oeflection  of  two-way  element 


Referring  to  figure  3-22,  the  deflected  shape  at  deflection  is  indicated  by 

the  solid  line  and  has  value  3  which  is  defined  as 

3  3  tan  “^ { X^/x )  =  tan  _1{H  tan  9max/2x)  3-21 

As  the  element  continues  to  deflect  the  angle  of  rotation  increases,  its 

magnitude  becoming  equal  to 

9h  =  A  +  3  3-22 

where  \  is  the  angular  rotation  in  excess  of  3  •  For  a  two-way  element  which 
undergoes  partial  failure  but  does  reach  incipient  failure,  the  maximum 
deflection  in  the  range  X^  ±  Xm  _<  Xu  becomes 

Xm  =  x  tan  9h  +  [ ( L/2)  -  x)]  tan  (9H  -  3)  3-23 

When  9h  equal  9  max,  the  ultimate  deflection  Xu  at  incipient  failure  is 

Xu  =  x  tan  9  max  +  [(L/2  -  x) ]  tan  [9max  -  tan_1(H  tan  9  max/2x)] 


Equations  3-19  through  3-24  are  specifically  for  two-way  elements  described  in 
figure  3-2  and  will  vary  for  other  two-way  elements  with  different  material 
properties  and  geometry. 

The  maximum  deflection  X_  for  several  two-way  elements  in  the  ranges  0  <  Xm 


■given  in  table  3- 
(X  )  deflections, 
indicated. 


urn  deflection  Xm  for  several  two-way  elements  in  the  ranges  0  <_  Xm 
i  X,  <  X_  <  Xu  as  a  function  of  the  rotation  angles  9H  and  9V  are 
table  3-o  aTong  with  the  values  of  partial  failure  (X^)  and  ultimate 
ections.  The  support  which  fails  at  partial  deflection  X1  is  also 


3-12  Elasto-Plastic  Resistance 

As  stated  in  Section  3-8,  the  initial  portion  of  the  resistance  function  (fig. 
3-1)  generally  is  composed  of  an  elastic  and  one  or  more  elasto-plastic 
ranges.  The  elastic  unit  resistance  rg  is  defined  as  the  resistance  at  which 
first  yield  occurs;  similarly,  the  elasto-plastic  unit  resistance  rg  is  the 
resistance  at  which  second  yields  subsequently  occur.  Where  all  hinges  form 
in  a  member  at  one  time,  re  will  be  equal  to  the  ultimate  unit  resistance  ru; 
where  two  or  more  hinges  are  formed  at  separate  times,  the  maximum  value  of 
rg  will  be  equal  to  ru  (depending  UDon  the  hinges  formed,  one  or  more  values 
or  rep  may  exist). 

These  resistances  for  one-way  elements  are  listed  in  table  3-7.  In  those 
cases  where  the  elasto-plastic  resistance  is  eaual  to  the  ultimate  resistance, 
the  value  can  be  determined  from  table  3-1. 


The  determination  of  the  elasto-plastic  resistances  of  two-way  elements  is 
more  complicated  than  that  for  one-way  elements,  since  the  resistance  varies 
with  the  span  ratio  and,  in  the  case  of  reinforced  concrete  elements,  with  the 
placement  of  the  reinforcement.  Data  for  calculating  the  resistances  of  two- 
way  elements  during  the  el asto-ol  astic  ranges  (graphically  summarized  in 
Figure  3-23)  are  presented  in  fiaures  3-24  through  3-38. 


Table  3-7  Elastic  and  Elasto-Plastic  Unit  Resistances  for  One-Way  Elements 


Edge  Conditions  and  Loading  Diagrams 


Elastic 

Resistance,^ 


Elasto  -  Plastic 
Resistance,  r>p 


FIG.  3-36  FIG.  3-37  FIG.  3-38 

LEGEND:  EDGE  CONDITIONS 
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Moment  and  deflection  coefficients  for  uniformly- loaded 
two-way  element  with  two  adjacent  edges  fixed  and  two 
edges  free 


Figure  3-24 
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Figure  3-25  Moment  and  deflection  coefficients  for  uniformly-loaded, 
two-way  element  with  one  edge  fixed,  an  adjacent  edge 
simply-supported  and  two  edges  free 
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Figure  3-26  Moment  and  deflection  coefficients  for  uniformly-loaded, 
two-way  element  with  two  adjacent  edges  simply-supported 
and  two  edges  free 
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Moment  and  deflection  coefficients  for  uniformly-loaded, 
two-way  element  with  three  edges  fixed  and  one  edge  free 
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Moment  and  deflection  coefficients  for  uniformly-loaded 
two-way  element  with  two  opposite  edges  fixed,  one  edge 
simply-supported  and  one  edge  free 


Figure  3-28 
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Moment  and  deflection  coefficients  for  uniformly-loaded, 
two-way  element  with  two  opposite  edges  simply-supported 
one  edge  fixed,  and  one  edge  free 
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Figure  3-30 


Moment  and  deflection  coefficients  for  uni formly-loaded , 
two-way  element  with  three  edges  simply-supported  and  one 
edge  free 
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Figure  3-31 


Moment  and  deflection  coefficients  for  uniformly-loaded 
two-way  element  with  two  adjacent  edges  fixed, 
one  edge  simply-supported,  and  one  edge  free 
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Moment  and  deflection  coefficients  for  uniformly-loaded, 
two-way  element  with  two  adjacent  edges  simply-supported 
one  edge  fixed,  and  one  edge  free 
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Moment  and  deflection  coefficients  for  uniformly-loaded 
two-way  element  with  all  edges  fixed 
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Moment  and  deflection  coefficients  for  uniformly-loaded 
two-way  element  with  two  opposite  edges  fixed  and  two 
edges  simply-supported 
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Figure  3-35  Moment  and  deflection  coefficients  for  uniformly-loaded, 

two-way  element  with  three  edges  fixed  and  one  edge  simply 
supported 


1.0 

0.7 

0.9 


(XI 

0.2 

0.1 

.07 

.09 

.03 

.02 

.01 

.007 

.009 

.003 

.002 


Figure  3-36  Moment  and  deflection  coefficients  for  uniformly-loaded,  two-way 
element  with  all  edges  simply-supported 
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Figure  3-37 


Moment  and  deflection  coefficients  for  uniformly-loaded ,  two-way 
element  with  two  adjacent  edges  fixed,  and  two  edges  simply- 
supported 
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Figure  3-38 


Moment  and  deflection  coefficients  for  um formly-loaded,  two-way 
element  with  three  edges  simply-supported  and  one  edge  fixed 


Figures  3-24  through  3-26  are  for  a  two-way  element  supported  on  two  adjacent 
sides  and  free  at  the  others.  Figures  3-27  through  3-32  are  for  a  two-way 
element  supported  on  three  sides  and  free  on  the  fourth,  while  figures  3-33 
through  3-38  are  for  elements  fixed  on  four  sides.  The  resistances  in  each 
range  can  readily  be  determined  using  the  coefficients  3,  the  subscript 
referring  to  the  points  listed  In  the  accompanying  illustration. 


For  example.  In  figure  3-27,  for  an  element  fixed  on  three  sides  and  free  on 
the  fourth.  If  the  ultimate  unit  resisting  moments  NL  are  known  for  points  1, 
2  and  3,  a  resistance  r  for  each  point  can  be  calculated  from 


r  = 


3-25 


where  the  values  of  3  are  found  in  figure  3-27.  The  smallest  value  of 
resistance  r,  ( say  at  point  2)  corresponds  to  the  first  yield  and  is  equal  to 
rA.  Next,  the  moments  at  the  remaining  two  points  are  computed  for  this  value 


>?' 


and  the  differences  between  these  and  the  ultimate  values  are 


determined.  These  differences  represent  the  remaining  moment  capacities 


available  for  additional  load.  At  rot  the  elements  supports  become  free. 


fixed,  and  simply-supported  on  opposite  sides.  Using  the  moment  differences 
and  entering  figure  3-26  two  values  of  the  change  in  resistance  can  be 
calculated  as  above,  the  smaller  being  Ar  and  therefore: 


ep 


=  re  +  A  r 


3-26 


In  similar  fashion,  the  resistance  at  the  end  of  each  ranqe  can  be  determined 
until  the  ultimate  unit  resistance  ru  is  reached. 

3-13  Elasto-Plastlc  Stiffnesses  and  Deflections 


The- slopes  of  the  elastic  and  elasto-plastic  ranges  of  the  resistance  function 
are  defined  by  the  stiffness  K  of  the  element: 


K  *  r/X 


3-27 


where  r  is  the  unit  resistance  and  X  is  the  deflection  corresponding  to  the 
value  of  r.  The  elastic  range  stiffness  is  denoted  as  Ke,  the  elasto-plastic 
range  as  Kep,  while  in  the  plastic  range  the  stiffness  is  zero. 


Typical  resistance-deflection  functions  used  for  design  are  shown  in  Figure 
3-39.  One-  and  two-step  systems  are  generally  used  for  one-way  elements  while 
two-  and  three-step  systems  are  used  for  two-way  elements.  Two  way  elements 
fixed  on  all  four  sides  will  exhibit  a  four  step  system.  As  can  be  seen  from 
the  figure,  the  elastic  range  stiffness 


K, 


re/Xe 


3-28 


the  elasto-plastic  stiffness  for  a  two-step  system 


'ep 


(ru  -  re)/(Xp  -  Xe) 
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and  for  a  three-step  system,  the  elasto-plastic  stiffnesses 


K 


ep 


(r 


ep 


re)/(Xep  -  xe) 
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RESISTANCE  —  RESISTANCE 


(a)  ONE  STEP  ELASTO-PLASTIC  SYSTEM 


TWO  STEP  ELASTO-PLASTIC  SYSTEM 


(c)  THREE  STEP  ELASTO-PLASTIC  SYSTEM 


Figure  3-39 


Resistance-deflection  functions  for  limited  deflections 


3-31 


K  ep  =  <ru  ‘  rep^xp  ”  xep^ 

The  elastic  and  elasto-plastlc  stiffnesses  of  one-way  elements  are  given  in 
table  3-8  as  a  function  of  the  modulus  of  elasticity  E,  moment  of  inertia  I, 
and  span  length.  Knowing  the  resistances  and  stiffnesses,  the  corresponding 
elastic  and  elasto-plastlc  deflections  can  be  computed  from  the  above 
equations. 

The  determination  of  the  elasto-olastlc  stiffnesses  and  deflections  of  two-way 
elements  Is  more  complicated  than  for  one-way  elements  since  another 
variable,  namely,  the  aspect  ratio  L/H,  must  be  considered.  For  two-way 
elements,  the  deflections  at  the  end  of  each  range  of  behavior  is  obtained 
from  the  y  coefficients  presented  in  figures  3-24  through  3-38.  The 
deflection  for  each  range  of  behavior  is  obtained  from 

XD  =  yrH4  3-32 

where  D,  the  flexural  rigidity  of  the  element  is  defined  as 

D  =  EI/ll-v2)  3-33 

E  Is  the  modulus  of  elasticity,  I  is  the  moment  of  inertia,  and  v  is  Poisson's 
ratio.  It  must  be  realized  that  except  for  the  elastic  range,  the  values  of  X 
(the  displacement)  and  r  In  Equation  3-32  represent  change  in  deflection  and 
resistance  from  one  range  of  behavior  to  another.  Therefore,  for  two-way 
members  the  change  in  deflection  and  resistance  (as  previously  explained)  is 
obtained  from  figures  3-24  through  3-38  and  the  stiffnesses  are  computed  from 
equations  3-28  through  3-33. 

3-14  Resistance-Deflection  Functions  for  Design 

3-14.1  General 

The  resistance-deflection  function  used  for  design  depends  upon  the  maximum 
permitted  deflection  according  to  the  design  criteria  of  the  element  being 
considered.  This  maximum  deflection  Xm  can  be  categorized  as  either  limited 
or  large.  In  the  limited  deflection  range,  the  maximum  deflection  of  the 
system  is  limited  to  the  elastic,  elasto-plastic  and  plastic  ranges.  When  the 
maximum  deflection  falls  in  the  large  deflection  range,  the  response  of  the 
system  is  principly  within  the  plastic  range  and  the  elastic  and  elasto- 
plastlc  ranges  need  not  be  considered.  The  error  resulting  from  the  omission 
of  the  elastic  and  elasto-plastic  portions  in  this  analysis  is  negligible. 

The  support  rotation  that  corresponds  to  limited  deflection  varies  for  the 
different  materials  used  in  protection  design.  The  response  criteria  for  each 
material  is  obtained  from  the  volume  that  describes  the  design  procedures  for 
that  material. 

3-14.2  Limited  Deflections 

When  designing  for  limited  deflections,  the  maximum  deflection  of  the 

element  is  kept  within  the  elastic,  elasto-plastic,  and  limited  plastic 
ranges,  and  the  resistance-deflection  function  for  design  takes  the  form  shown 


Table  3 
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In  figure  3-39  a,  b,  and  c  for  a  one-step  system,  a  two-step  system,  and  a 
three-step  system,  respectively.  The  design  charts  presented  In  Section  3- 
19.3  were  established  for  a  one-step  system;  for  two-  and  three-step  systems, 
these  charts  can  be  used  If  the  resistance-deflection  functions  are  replaced 
with  equivalent  elastic  resistance-deflection  functions  defined  by  KE  and  XE 
as  Indicated  by  the  dotted  lines  In  figure  3-39.  The  equivalent  elastic 
stiffness  KE  and  the  eqlvalent  maximum  elastic  deflection  XE  are  calculated 
such  that  the  area  under  the  dotted  curve  Is  equal  to  the  area  under  the  solid 
curve,  thereby  producing  the  same  potential  energy  in  each  system.  The 
equivalent  maximum  elastic  deflection  Xr  for  the  two-step  and  three-step 
systems  shown  Is  expressed  by  equations  3-34  and  3-35,  respectively. 

XE  =  Xe  +  Xp(l  -  re/ru)  3-34 

Xj:  =  xe( rep/ru^  +  xep^  ”  re/ru^  +  “  rep^ru^ 

3-35 

The  equivalent  elastic  stiffness  KE  in  each  case  is  equal  to 

Ke  =  ru/XE  3-36 

One-way  elements  exhibit  one-  and  two-step  resistance  deflection  curves 
depending  on  the  type  of  supports.  Consequently,  the  equivalent  elastic 
stiffness  Kc  is  given  for  one-way  elements  in  table  3-8.  The  equivalent 
elastic  deflection  can  then  be  calculated  from  equation  3-36.  Two-way 
elements  generally  exhibit  two-  and  three-step  resistance-deflection  curves 
which  are  a  function  of  not  only  the  type  of  supports  but  also  of  the  aspect 
ratio  L/H  of  the  element.  The  equivalent  elastic  deflection  XE  of  the  element 
under  consideration  must  be  calculated  from  equations  3-34  and  3-35  for  two- 
and, three-step  systems,  respectively.  The  value  of  KE  for  the  system  can  then 
be  obtained  from  equation  3-36. 

3-14.3  Large  deflections 

When  designing  for  a  large  deflection,  it  can  be  assumed  without  significant 
error  that  the  resistance  rises  instantaneously  from  zero  to  its  ultimate 
value  r,,  at  the  onset  of  the  blast  loading  thus  neglecting  the  elastic  and 
elasto-piastic  ranges.  The  design  resistance  function  for  one-way  elements  is 
approximated  by  a  constant  plastic  range  resistance  as  shown  in  figure  3-40a, 
while  for  two-way  elements,  the  resistance  function  becomes  that  shown  in 
figure  3-40b,  where  the  maximum  deflection  Xm  can  be  either  smaller  or  larger 
than  the  partial  failure  deflection  X^. 

3-15  Support  Shears  or  Reactions 

Support  shears  or  reactions  are  a  function  the  applied  load  and  the  maximum 
resistance  attained  by  an  element,  its  geometry  and  yield  line  location. 
However,  for  short  duration  blast  loads,  the  support  shears  can  be  reasonably 
estimated  by  neglecting  the  applied  load.  Therefore,  the  ultimate  support 
shear  can  be  assumed  to  be  developed  when  the  resistance  reaches  the  ultimate 
value,  ru. 
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3-40  Resistance-deflection  functions  for  large  deflections 
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Equations  for  the  ultimate  support  shears  V$  for  one-way  elements  are  given  in 
table  3-9.  For  those  cases  where  an  element  does  not  reach  its  ultimate 
resistance,  the  support  shears  are  obtained  based  upon  elastic  theory  for  the 
actual  resistance  r  attained  by  the  element. 

For  two-way  elements,  the  ultimate  shears  acting  at  each  section  are  calcu¬ 
lated  by  use  of  the  "yield  line  procedure"  previously  outlined  for  the  deter¬ 
mination  of  the  ultimate  resistance  ry.  The  shear  along  the  support  is 
assumed  to  vary  in  the  same  manner  as  the  moment  varies  (2/3  V  at  the  corners 
and  V  elsewhere)  to  account  for  the  higher  stiffness  of  the  corners  (fig.  3- 
41).  Since  the  shear  is  assumed  to  be  zero  along  the  interior  (usually 
positive)  yield  lines,  the  total  shear  at  any  section  of  a  sector  is  equal  to 
the  resistance  ru  times  the  area  between  the  section  being  considered  and  the 
positive  yield  lines.  Referring  to  fiqure  3-41,  the  support  shear  V$v  for  the 
triangular  sector  I  is 

( 2/3 ) V  sV ( L/4  +  L/4)  +  VsV(L/2)  = 


ru(Ly/2)  3-37 

VSV  «  3^/5  3-38 

and  for  the  trapezoidal  sector  II 

( 2/3 ) V sH (y/2 )  +  V sH ( H  -  y/2) 

ru(L/2)[(H  -  y  +  H)/2]  3-39 

VsH  *  [3ruL(2-y/H)]/2(6  -  y/H)  3-40 


Values  of  the  ultimate  support  shears  VsH  and  VsV  for  several  two-way  elements 
derived  as  above  are  presented  in  table's  3-10  and  3-11.  Table  3-10  gives  the 
ultimate  support  shears  for  the  case  where  opposite  supports  provide  the  same 
degree  of  restraint  resulting  in  symmetrical  yield  line  patterns.  Table  3-11 
is  for  completely  general  situations  where  opposing  supports  provide  different 
degrees  of  restraint  and  result  in  the  formation  of  unsymmetrical  yield 
lines.  .In  these  cases,  the  ultimate  support  shears  are  not  equal  at  opposing 
supports*. 

For  the  situations  where  the  ultimate  resistance  of  an  element  is  not 
attained,  the  support  shears  are  less  than  the  ultimate  value.  The  proportion 
of  the  total  load  along  each  support  as  well  as  the  distribution  of  the  shears 
along  the  support  is  assumed  to  be  the  same  as  previously  cited  for  the 
ultimate  support  shear.  Therefore,  the  support  shears  corresponding  to  the 
resistance  r  of  the  element  is  obtained  from  the  equations  presented  in  tables 
3-10  and  3-11  by  replacing  ru  with  the  actual  resistance  attained  (rg,  rg 
etc. ) . 

DYNAMICALLY  EQUIVALENT  SYSTEMS 

3-16  Introduction 

In  dynamic  analysis,  there  are  only  three  quantities  to  be  considered:  (1) 
the  work  done,  (2)  the  strain  energy,  and  (3)  the  kinetic  energy.  To  evaluate 
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Table  3-9  Support  Shears  for  One-Way  Elements 
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the  work  done,  the  displacement  at  any  point  on  the  structure  under  externally 
distributed  or  concentrated  loads  must  be  known.  The  strain  energy  is  equal 
to  the  summation  of  the  strain  energies  in  all  the  structural  elements  which 
may  be  in  bending,  tension,  compression,  shear  or  torsion.  The  kinetic  energy 
Involves  the  energy  of  translation  and  rotation  of  all  the  masses  of  the 
structure.  The  actual  evaluation  of  these  quantities  for  a  given  structure 
under  dynamic  load  would  be  complicated.  However,  for  practical  problems, 
this  can  be  avoided  by  using  appropriate  assumptions. 

In  most  cases,  a  structure  can  be  replaced  by  an  idealized  (or  dynamically 
equivalent)  system  which  behaves  timewise  in  nearly  the  same  manner  as  the 
actual  structure.  The  distributed  masses  of  the  given  structure  are  lumped 
together  into  a  number  of  concentrated  masses.  The  strain  energy  is  assumed 
to  be  stored  in  several  weightless  springs  which  do  not  have  to  behave 
elastically;  similarly,  the  distributed  load  is  replaced  by  a  number  of 
concentrated  loads  acting  on  the  concentrated  masses.  Therefore,  the 
equivalent  system  consists  merely  of  a  number  of  concentrated  masses  joined 
together  by  weightless  SDrings  and  subjected  to  concentrated  loads  which  vary 
with  time.  This  concentrated  mass-spring-load  system  is  defined  as  an 
equivalent  dynamic  system. 

Although  all  structures  possess  many  degrees  of  freedom,  one  mode  usually 
predominates  in  the  response  to  short  duration  loads;  thus,  for  all  practical 
purposes,  this  one  mode  may  be  considered  to  define  the  behavior  of  the 
structure  and  the  problem  can  be  simplified  by  considering  a  single-degree-of- 
freedom  system  whose  properties  are  those  of  the  fundamental  mode  of  the 
structures.  A  single-degree-of-freedom  system  is  defined  as  one  in  which  only 
one  type  of  motion  is  possible  or,  in  other  words,  only  one  coordinate  is 
required  to  define  its  motion.  Such  a  system  is  shown  in  figure  3-42  which 
consists  of  a  concentrated  mass  on  a  frictionless  surface  attached  to  a 
weightless  spring  and  subjected  to  a  concentrated  load.  A  single  displacement 
variable  x  is  sufficient  to  describe  its  motion.  The  material  presented 
herein  will  be  limited  to  single-degree-of-freedom  systems  only. 

Many  structures,  however,  exist  which  can  not  be  adequately  described  by  the 
first  vibration,  particularly  when  consideration  of  dynamic  strain  is 
necessary.  For  these  situations,  a  two  (or  more)-degree-of-freedom  analysis 
may  be  required.  Procedures  are  presented  in  section  3-19.2  for  performing 
these  analyses. 

Two  fundamental  methods  are  available  for  treatinq  simple  systems  subjected  to 
dynamic  forces.  The  first  of  these  methods  is  concerned  with  solving  the 
differential  equations  of  the  system  by  either  classical,  numerical  or 
graphical  means.  The  second  method  of  analysis  and  chart  solutions  depend  on 
solutions  which  have  been  determined  by  the  use  of  the  first  method  and  are 
approximate  solutions  to  the  problems  in  hand.  There  is  a  third  method  which 
may  also  be  very  useful  in  vibrational  problems;  this  involves  the  energy 
equation. 

In  u'.c  following  paragraphs  the  design  factors  used  to  determine  the 
equivalent  dynamic  system  that  is  used  to  analyze  structural  elements  are 
defined  and  the  methods  used  to  obtain  these  factors  are  discussed. 


3-17  Dynamic  Design  Factors 
3-17.1  Introduction 

The  values  of  the  mass  and  external  force  used  in  the  equation  of  motion 
(equation  3-1)  are  the  actual  values  only  if  all  the  elements  of  the  mass  of 
the  structure  experience  the  same  force  and,  consequently,  move  as  a  unit,  in 
which  case,  the  entire  mass  may  be  assumed  to  be  concentrated  at  its  center  of 
gravity.  In  other  cases,  the  assumption  of  uniform  motion  of  the  entire  mass 
cannot  be  made  without  introducing  serious  error.  This  is  true  for  members 
with  uniformly  distributed  mass  which  bend  or  rotate  under  load.  In  such 
cases,  the  motion  of  the  particles  of  mass  varies  along  the  length  of  the 
member.  Beams,  slabs,  etc.,  with  such  distributed  mass  actually  have  an 
infinite  number  of  degrees  of  freedom;  however,  such  structural  elements  can 
be  represented  by  an  equivalent  single-degree-of-freedom  system. 

In  order  to  define  an  equivalent  one-degree  system,  it  is  necessary  to 
evaluate  the  parameters  of  that  system;  namely,  the  equivalent  mass  M^,  the 
equivalent  spring  constant  l<£  and  the  equivalent  load  F^.  The  equivalent 
system  is  selected  usually  so  that  the  deflection  of  the  concentrated  mass  is 
the  same  as  that  for  a  significant  point  of  the  structure.  The  single-degree- 
of-freedom  approximation  of  the  dynamic  behavior  of  the  structural  element  may 
be  achieved  by  assuming  a  deflected  shape  for  the  element  which  is  usually 
taken  as  the  shape  resulting  from  the  static  application  of  the  dynamic 
loads.  The  assumption  of  a  deflected  shape  establishes  an  equation  relating 
the  relative  deflection  of  all  points  of  the  element. 

The  accuracy  of  the  computed  deformation  of  the  single-degree-of-freedom 
system  is  dependent  on  the  assumed  deflected  shape  and  the  loading  reqion. 
For  example,  in  the  impulsive  loadinq  realm,  if  the  assumed  deformed  shape  for 
a  'simply  supported  beam  is  taken  to  bo  a  parabola,  then  the  exact  solution 
will  be  obtained  for  the  deflection  of  the  beam.  In  the  quasi-static  loading 
realm,  however,  the  static  deformed  shape  gives  the  exact  answer.  For  blast 
resistance  design,  however,  the  differences  in  either  a  static  deformed  shape 
or  a  first  mode  approximation  is  negligible  and  thus  either  way  is 
permissible.  Usually,  though,  the  static  deformed  shape  is  easier  to  use  as 
it  covers  both  symmetric  and  asymmetric  deformation. 

To  permit  rapid  design  of  structural  elements  subjected  to  dynamic  loads,  it 
is  convenient  to  introduce  design  or  transformation  factors.  These  factors 
are  used  to  convert  the  real  system  to  the  equivalent  system.  To  obtain  them, 
it  is  necessary  to  equate  the  relations  expressing  the  kinetic  enerqy,  strain 
and  work  done  on  the  actual  structural  element  deflecting  according  to  the 
assumed  deflected  shape,  with  the  correspondi ng  relations  for  the  equivalent 
mass  and  spring  system. 

3-17.2  Load ,  Mass  and  Resistance  Factors 

3-17.2.1  Load  Factor.  The  load  factor  is  the  desiqn  or  transformation  factor 
by  which  the  total  load  applied  on  the  structural  element  is  multiplied  to 
obtain  the  equivalent  concentrated  load  for  the  equivalent  single-degree-of- 
freedom  system.  If  the  actual  total  load  on  the  structure  is  F  and  the 
equivalent  load  is  F^,  the  load  factor  is  defined  by  the  equation 


3-41 


kl  -  fe/f 

The  load  factor  Is  derived  by  setting  the  external ^  JJJ^e^y  the  actual 
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and  is  called  the  shape  function. 

NOTE.  The  shape  function,  ^x)  Is  different  for ^ L  different*"96 
plastic  range  and  therefore  the  load  factor,  KL,will  be  different. 

For  example;  the  shape  factor  for  a  ^ly  supporteo  beam  with  a  uh.formly 
distributed  load,  In  the  elastic  rangers  Is  defined 
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where  the  shape  factor  Is  defined  as 

i  3-46 

*r  *  2  6r/5max 
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Again  there  Is  a  <L  for  the  elastic  3*1  Seated  ^at  f  the  ^po  In  maximum 

However,  for  a  single  concentrated  load  located  v 
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deflection  (I.e.  the  center  of  a  simply  supported  or  fixed-fixed  beam,  or  the 
end  of  a  cantilever  beam)  <Pr  is  equal  to  one  for  both  the  elastic  and  plastic 
ranges  and  therefore  KL  is  equal  to  one  for  both  ranges. 

Values  for  Kj_  for  one-way  elements  are  presented  in  table  3-12.  Equation  3-41 
through  3-4b  were  used  to  calculate  these  values.  An  example  calculating  KL 
is  shown  in  Appendix  A. 

3-17.2.2  Mass  Factor.  The  mass  factor  is  the  design  or  transformation  factor 
by  which  the  total  distributed  mass  of  an  element  is  multiplied  to  obtain  the 
equivalent  lumped  mass  of  the  equivalent  si ngle-degree-of- freedom  system.  If 
the  total  mass  of  the  actual  element  is  M  and  the  mass  of  the  equivalent 
system  is  ME,  the  mass  factor  KM  is  defined  by  the  equation 

Km  =  Me/M  3-47 

Km  can  be  obtained  by  setting  the  kinetic  energy  of  the  equivalent  system 
equal  to  the  kinetic  energy  of  the  actual  structure  as  determined  from  its 
deflected  shape. 

For  a  structure  with  continuous  mass 

KE  =  l/2ME((J5max)2  =  l/2^m(x)[u>6(x)]2dx  3-48a 

where  u  =  natural  circular  frequency 

m(x)  =  distributed  mass  per  unit  length 

rearranging  terms 

ME  =^on’(x)(p2(x)dx  3-48b 

where  the  shape  function  <p(x)  is  based  on  the  deflected  shape  of  the  element 
due  to  the  applied  loading  and  not  to  the  distributionof  the  mass.  Since  the 
deflected  shape  of  the  element  is  different  for  the  elastic  and  plastic 
ranges,  <p(x),  and  therefore  KM,  will  also  be  different. 

Using  equations  3-47  and  3-48,  values  for  were  calculated  for  one-way 
elements  with  constant  mass.  These  values  are  shown  in  table  3-12. 

For  a  concentrated  mass  system, 

KE  =  l/2ME(u^,ax)2  =  1/22  Mr(w5r)2  3-49a 

r=l 

where 

Mr  =  rth  mass 

5r  =  deflection  of  mass  r 

i  =  number  of  lumped  masses 


rearranging  terms 


M 


E 


I  M  (C2 

i  ir^r 

r=l 


3-48b 
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Table  3-12  Transformation  Factors  for  One-Way  Elements 


Edge  Conditions  and  Loading  Diagrams 


L/2  i  L/2 


Range 

of 

Behavior 

Load 

Factor 

Kl 

Mass 

Factor 

Km 

Load-Mass 

Factor 

Klm 

Elastic 

0.64 

0.50 

0.78 

Plastic 

0.50 

0.33 

0.66 

Elastic 

1.0 

0.49 

0.49 

Plastic 

1.0 

0.33 

0.33 

Elastic 

0.58 

0.45 

0.78 

Elasto- 

Plastic 

0.64 

0.50 

0.78 

Plastic 

0.50 

0.33 

0.66 

Elastic 

1.0 

0.43 

0.43 

Elasto- 

Plastic 

1.0 

0.49 

0.49 

Plastic 

1.0 

0.33 

0.33 

Elastic 

0.53 

0.41 

0.77 

Elasto- 

Plastic 

0.64 

0.50 

0.78 

Plastic 

0.50 

0.33 

0.66 

An  example  of  calculations  for  finding  KM  can  be  found  in  Appendix  A. 

3-17.2.3  Resistance  Factor  .  Resistance  factor  is  the  design  factor  by  which 
the  resistance  of  the  actual  structural  element  must  be  multiplied  to  obtain 
the  resistance  of  the  equivalent  si ngl e-degree-of- freedom  system.  To  obtain 
the  resistance  factor,  it  is  necessary  to  equate  the  strain  energy  of  the 
structural  element,  as  computed  from  the  assumed  deflection  shape,  and  the 
strain  energy  of  the  equivalent  single-degree-cf-freedom  system.  If  the 
computed  total  resistance  of  the  structural  element  is  R  ana  the  equivalent 
total  resistance  of  the  equivalent  system  is  RE,  then  the  resistance  factor  is 
defined  by  the  equation 

Kr  =  Re/R  3-50 

Since  the  resistance  of  an  element  is  the  internal  force  tending  to  restore 
the  element  to  its  unloaded  static  position,  it  can  be  shown  that  the 
resistance  factor  KR  must  always  equal  the  load  factor  K|_. 

3-17.3  Load-Mass  Factor 

The  load-mass  factor  is  a  factor  formed  by  combining  the  two  basic 
transformation  factors,  and  KM.  It  is  merely  the  ratio  of  the  mass  factor 
to  the  load  factor,  and  it  is  convenient,  since  the  equation  of  motion  may  be 
written  in  terms  of  that  factor  alone.  The  equation  of  motion  of  the  actual 
system  is  given  as 

F  -  R  =  Ma 

and  for  the  equivalent  system 
KlF  -  KlR  =  KMMa 
which  can  be  re-written  as 

F  -  R  =  (KM/KL)Ma 

or  F  -  R  =  KLMMa  =  M0a 

where  KLM  =  load-mass  factor 

3  km/kl 

and  Me  =  effective  total 

when  expressed  in  terms  of  the  unit  area  of  the  element,  equation  3-52c  can  be 
written  as 

f  -  r  =  K|_Mma  =  m0a  3-54 

Values  of  load,  mass  and  load-mass  factors  are  presented  in  table  3-12. 
Equation  3-53  was  used  to  calculate  the  load  mass  factor. 

Instead  of  computing  the  several  factors  above,  the  load-mass  factors  in  the 
elastic  and  elasto-plastic  ranaes  can  be  determined  by  relatinq  the  primary 
mode  of  vibration  of  the  member  to  that  of  an  equivalent  si nql e-deqree-of- 
freedom  system.  In  the  elastic  range,  it  can  be  assumed  that  neither  the 


3-51 

3-52a 

3- 52b 
3-52c 

3-53 

mass  of  the  equivalent  system 


moment  nor  the  curvature  changes  between  the  plastic  hinges  under  increasing 
deflection.  This  behavior  results  in  a  linkage  action,  consideration  of  which 
can  be  used  to  evaluate  the  effective  plastic  mass  as  follows: 

In  figure  3-43,  a  portion  of  a  two-way  element  bounded  by  the  support  and  the 
yield  line  is  shown.  The  equation  of  angular  motion  for  this  section  is 

2M  =  Im0  3-55 

where  2m  =  summation  of  moments 

Im  =  mass  moment  of  inertia  about  the  axis  of  rotation 

9=  angular  acceleration 

Substituting  in  equation 

Fc  -  (2mn  +  2Mp)  =  (I^/l^la  3-56a 

where  a  Is  the  acceleration.  Dividing  through  by  c 

F  -  (2MN  +  2Mp)/c  =  ( I(n/cL1 )  a  3-56b 

Since  the  second  term  is  the  resistance  R 

F  -  R  =  (Im/cL^a  =  Mfia  3-56c 

and  the  load-mass  factor  KLM  for  the  sector  shown  is 

KLH  >  Im/cLlM  3-57 

where  c  Is  the  distance  from  the  resultant  applied  load  to  the  axis  of 
rotation,  is  the  total  length  of  the  sector  normal  to  the  axis  of  rotation 

and  M  is  the  total  mass  of  the  sector. 

When  the  element  is  composed  of  several  sectors,  each  sector  must  be 
considered  separately  and  the  contributions  then  summed  to  determine  the  load- 
mass  factor  for  the  entire  element. 

Klm  =  2(1^4  )/M  3-58 

For  elements  of  constant  depth  and,  therefore,  of  constant  unit  mass,  the 
above  equation  can  be  written  in  terms  of  the  area  moment  of  inertia  I  and  the 
areas  A  of  the  individual  sectors. 

Klm  =  SU/cLil/A  3-59 

Table  3-13  gives  the  load-mass  factors  for  the  elastic  and  elasto-plastic 
response  ranges  for  various  two-way  elements.  The  values  given  in  the  table 
apply  for  two-way  elements  of  uniform  thickness  and  supported  as  indicated. 
Also  if  there  are  openings  in  the  element,  they  must  be  compact  in  shape  end 
small  in  area  compared  to  the  total  area  of  the  element. 


Figure  3-44  gives  the  load-mass  factors  for  the  plastic  response  of  various 
two-w^  elements.  The  load-mass  factors  are  a  junction  of  support  condition 
and  yield  line  location.  These  two-way  elements  must  conform  to  the  following 
limitations: 

(1)  Element  must  be  of  constant  thickness. 

(2)  Opposite  support  must  provide  the  same  restraint  so  that  a 
symmetrical  yield  line  pattern  occurs. 

(3)  Any  openings  must  be  compact  In  shape  and  small  in  area  compared 
to  the  total  area  of  the  element. 

3-17.4  Natural  Period  of  Vibration 

To  determine  the  maximum  response  of  a  system  either  by  numerical  methods  or 
design  charts  (both  methods  are  described  in  Section  3-19),  the  effective 
natural  period  of  vibration  Is  required.  This  effective  natural  period  of 
vibration,  when  related  to  the  duration  of  a  blast  loading  of  given  intensity 
and  a  given  structural  resistance,  determines  the  maximum  transient  deflection 
Xm  of  the  structural  element. 

The  effective  natural  period  of  vibration  is 

Tn  =  2n(me/KE)1/2  =  2Tt(KLMm/KE)1/2  3-60 

where  me  =  the  effective  unit  mass 

Ke  =  the  equivalent  unit  stiffness  of  system 

The  values  used  for  the  effective  mass  and  stiffness  for  a  particular  element 
depends  on  the  allowable  maximum  deflection.  When  designing  for  completely 
elastic  behavior,  the  elastic  stiffness  is  used.  In  all  other  cases  the 
equivalent  elasto-pl astic  stiffness,  KE,  is  used.  The  elastic  value  of  the 
effective  mass  is  used  for  the  elastic  range,  while  in  the  el asto-plastic 
range  the  effective  mass  is  the  average  of  the  elastic  and  elasto-plastic 
values.  For  small  plastic  deflections,  the  value  of  the  effective  mass  is 
equal  to  the  average  of  the  plastic  value  and  the  equivalent  elastic  value. 
The  plastic  effective  mass  is  used  for  large  plastic  deflections. 


OYNAMIC  ANALYSIS 

3-18  Introduction 

Structural  elements  must  develop  an  internal  resistance  sufficient  to  maintain 
all  motion  within  the  limits  of  deflection  prescribed  for  the  particular 
design.  The  load  capacity  of  the  member  depends  on  the  peak  strength 
developed  by  the  specific  member  (see  Volumes  IV  and  V)  and  on  the  ability  of 
the  member  to  sustain  its  resistance  for  a  specific  though  relatively  short 
period  of  time.  This  section  describes  the  methods  used  for  the  analysis  of 
structural  elements  subjected  to  dynamic  loads,  and  these  elements  are  divided 
into  three  groups,  namely:  (1)  elements  which  respond  to  the  pressure  only 
(low  pressure  range),  (2)  elements  which  respond  to  pressure-time  relationship 
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( intermedi ate-pressure  design  range)  and  (3)  elements  which  respond  to  the 
impulse  which  corresponds  to  the  high  pressure  design  range. 

The  method  of  analysis  is  different  for  the  groups  of  elements.  Those 
elements  that  respond  to  pressure  only  (this  is  still  a  pressure-time 
relation)  and  those  that  respond  to  pressure-time  relationship  may  be  analyzed 
using  either  numerical  methods  or  design  charts,  while  elements  that  respond 
to  the  impulse  are  analyzed  using  either  design  charts  for  large  deflection  or 
an  Impulse  method  for  designs  with  limited  deflections. 

The  method  of  analysis  employing  numerical  techniques  provides  a  solution  that 
is  obtained  by  a  step-by-step  integration,  and  is  generally  applicable  for  any 
type  of  load  and  resistance  function.  This  method  of  analysis  is  presented  in 
this  section.  Also  presented  is  the  result  of  a  systematic  analysis  of 
single-degree  dynamic  systems  for  several  idealized  loadings  plotted  in  non- 
dimensional  design  charts. 

3-19  Elements  Which  Respond  to  Pressure  Only  and  Pressure-Time  Relationship. 
3-19.1  General 

To  analyze  these  elements,  the  true  magnitudes  for  the  pressure-time 
relationship  of  the  applied  blast  loads  must  be  known.  The  determination  of 
the  dynamic  response  of  these  systems  is  accomplished  using  numerical 
techniques  or  design  charts  which  relate  the  dynamic  properties  of  the  element 
(natural  period  of  vibration,  resistance,  and  deflection)  to  those  of  the 
blast  overpressures. 

3-19.2  Analysis  by  Numerical  Methods 

The  numerical  method  of  analysis  involves  a  step-by-step  integration 
procedure,  starting  at  zero  time,  when  displacement  and  the  velocity  of  the 
system  are  presumably  known.  The  time  scale  is  divided  into  several  discrete 
intervals,  and  one  progresses  by  successively  extrapolating  the  displacement 
from  one  time  station  to  the  next.  As  the  time  interval  in  the  sequence  of 
time  is  reduced,  i . e . ,  as  the  number  of  discrete  intervals  is  increased,  the 
accuracy  of  the  numerical  method  is  improved. 

3-19.2.1  ' Single-Degree-0 f -Freedom  Systems  .  As  stated  previously,  a  single 
degree  of  freedom  system  is  one  in  which  only  one  coordinate  is  essential  to 
define  its  motion.  The  equation  of  motion  for  the  single-degree-of-freedom 
system  shown  in  figure  3-42  can  be  expressed  as: 

F  -  kx  -  cv  =  Ma  3-61 

when  c  =  damping  constant 

There  are  several  methods  that  can  be  used  to  numerically  integrate  equation 
3-61,  but  the  two  procedures  applied  to  this  problem  will  be  the  averaqe 
acceleration,  and  the  acceleration-impulse  extrapol ation  methods. 

3-19.2.1.1  Average  Acceleration  Method  .  In  the  average  acceleration  method, 
the  velocity  and  displacement  at  any  time  t  are  expressed  as 

vt  =  vt-l  +  1/2  (at  +  at-l5  at  3-62 
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3-63 


x  3  +  1/2  (vt  +  vt-1)  At 

Substituting  equation  3-62  Into  equation  3-63, 

x  3  xt-1  +  vt-1  fit  +  1/4  (at  +  a^j)  at2  3-64 

Substituting  equations  3-62  and  3-63  Into  equation  3-61,  we  have 

F  -  k  [xt_i  +  1/2  (vt  +  vt-1)  at]  -  c  Cvt_i  +  1/2  (at  +  at_ ^ )  at]  =  Ma 


Equation  3-65  can  be  simplified  to 


m  +  (l/2)c  at  +  (l/4)k  At* 


[F  -  k  (xt_1  +  vt-1  at  +  l/4at_1  at2) 


-  c  (vt-1  +  1/2  at-1  at)]  3-66 

The  Integrating  procedure  can  be  summarized  as  follows: 

Step  1.  At  t  3  0:  Compute:  a0  using  equation  3-61  and  specified  values  for 
x0(t  =  o),  v0(t  3  o),  and  F(t  3  o).  . 

Step  2.  Increment  time:  t  3  t  +  at 

Step  3.  At  t  3  t  +  at:  Compute  a{t  3  t  +  at)  using  equation  3-66, 

v  (t  =  t  +  at)  using  equation  3-62  and  x  (t  3  t  +  at)  using  equation 
3-63. 


Step  4.  Repeat  steps  2  and  3  until  some  specified  time  or  until  the  maximum 
displacement  is  reached. 


It  should  be  noted  that  if  x  (t  =  0)  and  v  (t  3  0)  are  zero,  then  the  system 
will  not  move  as  long  as  F  (t  =  0)  is  zero;  however,  if  F  (t  =  0)  is  not  zero, 
Equation  3-66  will  give  an  acceleration  and  motion  will  start.  Thus  the  system 
of  equation  is  self-starting  and  requires  no  special  starting  provisions. 

The  procedure  outlined  above  for  an  elastic  single-degree-of-freedom  system  can 
be  easily  extended  to  accomodate  elastic-plastic  behavior.  To  do  so,  equation 
3-61  is  re-written  to  include  the  restoring  force  R  which  replaces  the  spring 
restoring  force  kx,  Equation  3-61  becomes 

F  -  R  -  cv  =  ma  3-67 


Substituting  equations  3-62  and  3-63  into  equation  3-67,  the  following  equation 
for  acceleration  is  obtained: 


m  +  cflt 


[F  -  R  -  c  (vt  i  +  A-at--.1 ] 

2 


Unlike  equation  3-66,  equation  3-68  has  a  non-linear  term  R,  and  it  depends  on 
the  displacement  x  at  time  t.  Therefore,  instead  of  using  the  direct 
Integration  approach  as  in  the  previous  case,  a  predictor-corrector  method  is 
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utilized.  Briefly,  the  displacement  x  at  time  t  Is  estimated  (predicted)  and 
then  corrected  and  covergence  of  this  procedure  can  be  obtained  In  a  single 
Iteration  If  the  value  of  At  Is  small  enough.  The  step-by-step  procedure  is 
outlined  as  follows: 

Step  1.  At  t  *  0:  Compute  at  (t  *  t^)  from  equation  3-67  and  the  Initial 
conditions. 

Step  2.  Increment  time:  t  3  t  +  At 

Step  3.  At  t  3  t  +  At:  Set  at  3  at-1  (t  -  1  3  0,  Initially) 

Step  4.  Compute  and  xt  from  equations  3-62  and  3-63 

Step  5.  Compute  R 

Step  6.  Compute  at  from  equation  3-68 

Step  7.  If  you  are  In  the  predicting  pass,  return  to  Step  4,  if  in  the 

correcting  pass,  set  xt_j  3  xt,  3  vt,  at  3  at.j  and  go  to  Step  2. 

3-19.2.1.2  Acceleration  Impulse  Extrapolation  Method  ,  Suppose  the 

acceleration  of  the  system  Is  defined  by  figure  3-45a.  The  acceleration- 
impulse  extrapolation  method  assumes  that  the  actual  acceleration  curve  can  be 

replaced  by  a  series  of  equally  spaced  impulses  occurring  at  tQ,  t^,  . 

tn,  as  shown  In  figure  3-45b. 

The  magnitude  of  the  acceleration  impulse  at  tn  is  given  by 

I(tn)  3  an(  At)  3-69 

where  t-  t,  -  t  3  to  -  t,  3  .  .  .  3  tn  -  tn  ..  This  is  shown  in  figure  3- 

45b.  Since  an  impulse  is  applied  at  tn,  there  is  a  discontinuity  in  the  value 

of  velocity  at  tn.  In  the  time  Interval  from  t-  to  tn_^,  the  velocity  is 
constant  and  the  displacement  varies  linearly  witn  time.  The  velocity  and 
displacement  thus  obtained  are  shown  In  figures  3-45c  and  3-45d. 

Suppose  tl  and  tt  Indicate  the  time  Immediately  before  and  after  the 
application  of  the  impulse  at  tn,  and  let  vl  and  vn  indicate  respectively  the 
velocity  at  tp  and  t^;  these  two  velocities  are  related  by  the  following 
equation: 


*?!  1  at> 

The  relationship  between  xn_j  and  xn, 

xn  -  V1  =  V  4t> 
xn*l  -  xn  ’ 

Combining  equations  3-70  and  3-71, 
related  by: 

xn+l  =  2xn  -  xn-l  +  an(w>2 


3-70 

and  between  xn  and  xn+1  are  given  by: 

3-71  a 

3-71  b 

the  three  successive  displacements 

3-72 


are 
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TIME  (t) 


(c)  VELOCITY 


(d)  DISPLACEMENT 


Figure  3-45  Acceleration-impulse  extrapolation  method 
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This  Is  the  basic  recurrence  formula  for  the  acceleration  impulse  extrapolation 
method.  Once  the  values  of  x  at  tp_^  and  tn  are  known,  the  value  at  tn+1  can 
be  directly  computed  without  resorting  to  a  trial  and  error  procedure. 

It  Is  necessary,  however,  to  use  a  special  procedure  in  the  first  time  interval 
because,  at  t  =  0,  no  value  of  xn-1  is  available.  Two  different  procedures  may 
be  used.  The  first  assumes  the  acceleration  varies  linearly  up  to  the  first 
time  station,  in  which  case,  the  displacement  of  the  system  at  that  time  can  be 
expressed  as: 

xj  =  (1/6)  (2a0  +  ai)(Z\t)2  3-73 

The  second  procedure  assumes  that  the  acceleration  is  constant  during  the  first 
time  interval  and  is  equal  to  the  initial  value.  This  assumption  results  in 
the  following  expression  for  x^: 

Xj  =  (l/2)a0(At)2  3-74 

It  should  be  noted  that  equation  3-73  has  to  be  solved  by  trial  and  error  since 
aj  depends  on  Xi.  This  trial -and-error  procedure  can  be  avoided  by  taking  the 
value  of  ai  to  be  f(t)/m  instead  of  the  exact  value  of  C f ( t j)  -  R(ti)  ]/m.  In 
figure  3-46a  if  a^  =  0,  aQ  is  equal  to  ag/6  and  in  figure  3-46b  if  a0  and  a^ 
are  approximately  equal  to  each  other,  the  value  a'0/2  may  be  used  for  aQ.  The 
recurrence  formula  is  directly  applicable  for  the  evaluation  of  xp,  X3,  ..., 
x0+l  whenever  the  acceleration  is  a  continuous  function.  When  there  is  a 
discontinuity  in  the  acceleration,  this  equation  can  still  be  used  but  not 
before  some  modifications  are  made  to  reflect  this  discontinuity.  In  figure  3- 
46a  there  is  a  discontinuity  at  t4  which  occurs  at  the  end  of  the  time  interval 
fit.  Under  this  condition,  the  value  of  an  used  in  the  numerical  procedure  is 
the  average  value  at  discontinuity,  namely  1/2  (a^  +  a£)  for  a4.  In  figure  3- 
46b  the  discontinuity  occurs  within  time  interval,  t.  The  correct  value  for 
a3  in  this  procedure  is: 

a3  =  ( 1/2) [a3  +  a$  +  (a$  -  a$)  (At'/flt)2]  3-75 

When  there  are  more  than  two  discontinuities,  this  procedure  becomes  cumbersome 
and  it  is  more  convenient  to  replace  the  given  acceleration  by  a  smooth 
curve.  Equation  3-73  must  be  used  if  there  is  zero  force  (and  hence  zero 
acceleration)  at  zero  time,  for  in  no  other  way  can  x^  be  determined.  If 
acceleration  at  t  =  0  is  not  zero,  equation  3-74  may  be  used  without 
appreciable  error,  provided  the  force  does  not  change  greatly  in  the  first 
interval . 

The  tedious  method  of  numerical  analysis  is  made  somewhat  easier  if  the 
computation  procedure  is  set  up  as  shown  in  table  3-14.  In  the  first  time 
interval,  the  value  used  for  the  acceleration,  aQ,  depends  on  the  forcing 
function.  If  the  value  of  the  forcing  function  is  zero  at  time  zero,  then  aQ 
is  taken  as  one-sixth  of  a^ ,  where  a,  is  the  ratio  of  the  value  of  the  forcinq 
function  at  t,  to  the  mass  of  the  system  [a0  =  f ( t.j_ ) /6m] .  When  the  value  of 
the  forcing  function  is  not  zero  at  time  zero,  and  if  the  values  of  the  forcing 
function  at  times  t=o  and  t=tj  are  approximately  equal,  then  aQ  may  be  taken  as 
(P0-R0)/2.  The  recurrence  formula  (equation  3-72)  is  used  in  the  other 
interval s. 
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Table  3-14  Details  of  computation  by  acceleration  impulse  extrapolation  method. 


Xn+i  =  cXn- 
Atj  =Tn  / 10 
At2  =2At, 


It  may  be  necessary  during  the  computational  procedure,  to  increase  the  value 
of  t  so  as  to  reduce  the  number  of  time  Intervals.  In  figure  3-47,  for 
example,  the  forcing  function  varies  linearly  after  time  j.  For  the  next  time 

station  (j  +  list,  the  value  of  Ati  can  be  Increased  (At2  =  2  At],  for  example) 

without  Introducing  any  errors  In  tne  computations.  In  such  cases,  caution  has 
to  be  exercised  when  using  the  recurrence  formula.  In  table  3-14,  the  value 
of  At  changes  after  the  jth  step,  to  At2  =  2Atj.  In  the  (j  +  list  step,  xn-1 

should  be  the  value  of  x  for  the  (j-l)nd  step  and  xn  will  be  the  value  of  xn+1 

for  the  (j-l)st  step. 

3-19.2.2  A  Two-Degree-Of-Freedom  System  .  A  two-degree-of-f reedom  system  is 
shown  in  figure  3-48  where  the  coordinates  defining  the  configuration  of  the 
system  are  Xj  and  X2.  The  equation  of  motion  for  each  of  the  two  masses  can  be 
expressed  as  follows: 

Fi  +  R2  -  Rj^  =  M^a^  3-76a 

F2  -  Rj>  =  M2a2  3-76b 

Equation  3-76  can  be  rewritten  as: 


ax  =  F1/M1  +  lyMj  -  R1/M1  3-77a 

a2  =  F2/M2  -  R2/M2  3-77b 

The  integration  procedures  described  for  a  single-degree-of-f reedom  system  can 
also  be  applied  to  two-degrees-of-freedom  systems.  In  the  average  acceleration 
procedure.  If  the  equations  are  placed  in  the  same  form  as  equation  3-61,  i.e. , 
If  the  acceleration  at  time  t  are  written  in  terms  of  displacements  and 
velocities  at  time  t  -  At,  then  two  equations  are  obtained  which  must  he  solved 
simultaneously  for  at  and  a2.  In  the  other  numerical  integration  procedure, 
the  recurrence  formula  (equation  3-72)  applies  to  each  of  the  two  displacements 
independently. 

3-19.2.3  Damping  .  The  effects  of  damping  are  hardly  ever  considered  in  blast 
design  because  of  the  following  reasons: 

(1)  Damping  has  very  little  effect  on  the  first  peak  of  response 
which  is  usually  the  only  cycle  of  response  that  is  of  interest. 

(2)  The  energy  dissipated  through  plastic  deformation  is  much  greater 
than  that  dissipated  by  normal  structural  damping. 

(3)  Ignoring  damping  is  a  conservative  approach. 

If  damping  has  to  be  considered  in  an  analysis,  however,  it  should  be  expressed 
as  some  percentage  of  critical  damping.  For  free  vibration,  this  is  the  amount 
of  damping  that  would  remove  all  vibration  from  the  system  and  allow  it  to 
return  to  its  neutral  position.  Critical  damping  is  expressed  as 

Ccr  =  2  kM  3-78 

where  k  =  stiffness  of  the  system 

M  =  mass  of  system 


The  damping  coefficient,  c,  in  equations  3-61  and  3-67  is  expressed  as  a 
percentage  of  critical  damping,  Ccr  For  steel  structures,  c  should  be  taken 
as  0.05Cr  and  0.01  Ccr  for  reinforced  concrete  structures. 

3-19.3  design  Charts  for  Idealized  Loadings. 

3-19.3.1  General  .  The  response  of  single-degree-of-freedom  systems  subjected 
to  Idealized  blast  loadings  is  presented  in  the  form  of  non-dimensional 
curves.  In  order  to  utilize  these  response  charts,  both  the  blast  loads 
(pressure  -  time  history)  and  the  resistance-deflection  curve  of  the  structural 
system  are  idealized  to  linear  or  bilinear  functions.  Methods  for  computing 
these  Idealized  blast  loads  are  given  in  Volume  II  while  the  methods  for 
computing  the  resistance-deflection  functions  as  well  as  converting  the  actual 
system  to  an  elastic  or  elasto-plastic  single-degree-of-freedom  system  have 
been  presented  in  previous  sections  of  this  volume. 

The  response  of  a  structural  system  subjected  of  a  dynamic  load  is  defined  in 
terms  of  its  maximum  deflection  Xm  and  the  time  t^  to  reach  this  maximum 
deflection.  The  dynamic  load  is  defined  by  its  peak  value  P  and  duration  T 
while  the  single-degree-of-freedom  system  is  defined  in  terms  of  its  ultimate 
resistance  r  ,  elastic  deflection  Xc  and  natural  period  Tu.  Response  charts 
relate  the  dynamic  properties  of  trie  blast  load  (P  and  T)  to  those  of  the 
element  (r,  X,  TN ) ,  that  is,  Xm/X^-  and  t^/T  are  plotted  as  a  function  of  ru/P 
and  T/Tn. 

Response  charts  have  been  prepared  for  simplified  loads  as  well  as  for  the  more 
complex  bilinear  loadings.  The  simplified  loadings  include  triangular, 
retangular,  step  load  with  finite  rise  time,  triangular  with  rise  time  and 
sinusoidal  pulse.  The  idealized  triangular  loading  is  utilized  in  the  analyses 
of  acceptor  structures.  These  structures  are  subjected  to  the  effects  of 
unconfined  explosions,  such  as  free-air  or  surface  bursts  or  partially  confined 
explosions  in  donor  structures  which  are  fully  vented.  On  the  other  hand, 
idealized  instantaneously  applied  or  gradually  applied  flat  top  loads  are 
usually  encountered  in  the  design  of  roof  elements  of  acceptor  structures. 
These  structures  are  subjected  to  very  long  duration  loads  produced  by  very 
large  explosive  quantities.  The  triangular  loading  with  a  rise  time  is  similar 
to  the  previously  mentioned  triangular  load  except  that  it  is  an  equivale.  t 
loading  produced  by  the  blast  traversing  an  element.  The  sinusoidal  loadings 
is  usually  associated  with  the  vibrational  type  of  load  rather  than  blast 
induced  input.  The  bilinear  loading  is  utilized  in  the  analysis  of  containment 
structures  which  allow  partial  venting  of  an  explosion  and  in  the  analysis  of 
acceptor  structures  where  reflected  pressures  clear  rapidly  around  the 
structure. 

The  effects  of  damping  have  not  been  included  in  the  preparation  of  the 
response  charts.  In  the  design  of  the  blast  resistant  structures,  the  first 
peak  of  response  is  usually  the  only  cycle  of  resonse  that  is  of  interet  since 
the  maximum  resistance  and  deflection  is  attained  in  that  cycle.  Damping  has 
very  little  effect  on  this  first  peak  and  consequently,  neqlectinq  damping  has 
negligible  effects  on  maximum  response  calculations. 

3-19.3.2  Maximum  Response  Of  Linear  Elastic  System  To  Simplified  Loads  _  T0 
obtain  the  response  of  a  linear  elastic  system,  it  is  convenient  to  consider 
the  concept  of  the  dynamaic  load  factor.  This  factor  is  defined  as  the  ratio 
of  the  maximum  dynamic  deflection  to  the  deflection  which  would  have  resulted 


from  the  static  application  of  the  peak  load  P,  which  Is  used  in  specifying  the 
load-time  variation.  Thus  the  dynamalc  load  factor  (DLF)  is  given  by: 

DLF  =  Xm/X$  3-79a 

where  X$  *  static  deflection  or.  In  other  words,  the  displacement  produced 
In  the  system  when  the  peak  load  is  applied  statically. 

Xm  =  maximum  dynamic  deflection 

Since  deflections,  spring  forces,  and  stresses  in  an  elastic  system  are  all 
proportional,  the  dynamic  load  factor  may  be  applied  to  any  of  these  to 
determine  the  ratio  of  dynamic  to  static  effects.  Therefore,  the  dynamic  load 
factor  may  be  considered  as  the  ratio  of  the  maximum  resistance  attained  to  the 
peak  load  or: 

DLF  =  r/p  3-79b 

where  r  =  maximum  dynamic  resistance 

p  =  peak  load  used  in  specifying  the  load-time  variation. 

For  a  linear  elastic  system  subjected  to  a  simplied  dynamic  load,  the  maximum 
response  is  defined  by  the  dynamic  load  factor,  DLF  and  maximum  response  time, 
tm.  The  dynamic  load  factor  and  time  ratio  WT  are  plotted  versus  the  time 
ratio  T/Tn  for  a  triangular  load,  rectangular  load,  step  load  with  finite  rise 
time  (Tr),  triangular  load  with  rise  time,  and  sinusoidal  pulse  in  figures  3-49 
through  3-53. 

In  many  structural  problems  only  the  maximum  value  of  the  DLF  is  of  interest. 
For  the  most  prevalent  load  case,  namely,  the  triangular  load  as  well  as  the 
retangular  and  step  load  with  rise  time,  the  maximum  value  of  the  DLF  is  2. 
This  immediately  indicates  that  all  maximum  displacements,  forces,  and  stresses 
due  to  the  dynamic  load  are  twice  the  value  that  would  be  obtained  from  a 
static  analysis  for  the  maximum  load  P. 

The  above  response  charts  apply  for  elastic  systems.  However,  the  charts  can 
be  applied  to  the  entire  elasto-plastic  range  if  the  actual  resistance- 
deflection  curve  (two  step  system,  three  step  system,  etc.)  is  replaced  by  the 
equivalent  elastic  system  where  the  equivalent  elastic  stiffness  KE  and  the 
equivalent  elastic  deflection  X^  are  obtained  according  to  previously  explained 
procedures. 

In  a  typical  design  example,  the  pressure-time  loading  is  calculated  and 
idealized  to  one  of  the  simplified  loads  defined  by  P  and  T  or  Tr.  A 
structural  member  is  assumed  and  its  corresponding  dynamic  properties  are 
calculated.  For  a  completely  elastic  response,  the  values  of  r.  and  T\j  are 
calculated  while  for  a  response  in  the  elasto-plastic  range,  ru,  and  are 
obtained.  It  should  be  noted  that  for  a  response  in  the  elasto-plastic  range, 
the  value  of  Tu  is  calculated  using  the  effective  mass  which  is  an  average  of 
the  elastic  ana  elasto-plastic  values.  Knowing  the  ratio  of  T/Tij,  the  dynamic 

load  factor  (DLF)  and  the  time  ratio  t^/T  can  be  read  from  tne  appropriate 

figures.  The  maximum  resistance  r  attained  by  the  structural  member  is 

calculated  from  the  DLF  and  the  response  time  is  obtained  from  the  time 
ratio.  If  the  resistance  r  is  greater  than  re  for  a  completely  elasti.c 

response  in  the  elasto-plastic  range,  then  the  analysis  is  not  valid  and  the 


t/tn 
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igure  3-49  Maximum  response  of  elastic  one-degree-of-f reedom  system  for 
triangular  load 


Dynamic  Load  Factor 


Dynamic  Load  Factor 


Figu 


Figure  3-53  Maximum  response  of  elastic  one-degree-of-freedom  system  for  sinusoidal 
pulse  load 


procedure  Is  repeated.  The  maximum  deflection  Is  obtained  from  the  resistance 
r  and  the  stiffness  Kg  or  KE. 

3-19.3.3  Maximum  Plastic  Response  Of  An  Elasto-Plastic  System  To  Idealized 
Loads.  An  elasto-plastic  system  may  have  an  elastic  or  plastic  response 
depending  upon  the  magnitude  of  the  blast  load.  If  the  response  is  elastic, 
that  Is,  the  member  attains  a  resistance  r  which  Is  less  than  its  ultimate 
resistance  ru,  then  the  charts  of  the  preceding  section  are  used.  The  response 
charts  presented  In  this  section  are  only  for  the  plastic  response  of  members 
where  the  ultimate  resistance  ru  is  attained. 


The  maximum  plastic  response  of  an  elasto-plstic  system  subjected  to  a  blast 
load  is  defined  by  the  maximum  deflection,  Xm  it  attains  and  the  time,  t,,,  it 
takes  to  reach  this  deflection.  The  blast  load  is  defined  by  its  peak  value  P 
and  duration  T  while  the  slnge-degree-of-freedom  system  is  defined  by  its 


ultimate  resistance  ru,  elastic  deflection  Xr  and  natural  period  TN.  A  non- 
dimensional  response  cnart  is  constructed  by  plotting  the  ductility  ratio  Xm/XE 
and  the  time  ratio  t^/T  as  a  function  of  rj?  and  T/TN.  Response  charts  are 
given  for  a  triangular  load,  retangular  load,  step  load  with  finite  rise  time 


Tr  and  triangular  load  with  rise  time  in  figures  3-54  through  3-61.  It  should 
be  noted  that  for  the  step  load  with  finite  rise  time,  the  load  is  defined  by 


the  rise  time  Tr  and,  consequently,  the  response  curves  are  plotted  using  Tr 
rather  than  T. 


In  a  typical  example,  the  pressure-time  loading  is  calculated  and  idealized  to 
one  of  the  simplified  loads  defined  by  either  P  and  T  or  P  and  Tr.  A 
structural  member  is  assumed  and  is  corresponding  dynamic  properties  (r  Xr, 
Tn)  are  calculated.  Knowing  the  ratios  of  ru/P  and  T/TN  (or  Tr/TNJ,  the 
ductility  ratio  Xm/XE  and  the  time  ratio  t_,/T  can  be  read  from  the  appropriate 
figures.  The  maximum  deflection  Xm  ana  response  time  t,-  can  readily  be 
calculated.  If  the  ductility  ratio  X^/Xr-  or  the  maximum  deflection  Xm  and 
corresponding  response  time  ^  are  unsatisfactory,  the  procedure  is  repeated. 

It  should  be  noted  that  the  value  of  the  natural  period  of  vibration  T^  used  in 
conjuntion  with  the  response  charts  varies  according  to  the  magnitude  of  the 
ductility  ratio  Xm/XE.  For  small  plastic  deformations  (Xm/XE  less  than  5),  the 
calculations  of  TN  is  based  on  an  average  of  the  equivalent  elastic  and  plastic 
masses.  Whereas,  for  larger  plastic  deformations  (Xm/XE  greater  than  5),  the 
equivalent  plastic  mass  is  used  to  obtain  T^. 

3-19.3.4  Maximum  Plastic  Response  Of  An  Elasto-Plastic  System  to  Idealized  Bi¬ 
linear  Loads  .  Response  charts  have  been  prepared  for  bilinear  loads  in  much 
the  same  manner  as  for  simplified  loads.  However,  four  parameters  are  required 
to  define  the  bilinear  loading  rather  than  two  parameters  which  are  required 
for  the  simplified  loads.  These  additional  parameters  greatly  increase  the 
number  of  charts  required  for  the  rapid  prediction  of  the  dynamic  response  of 
an  elasto-plastic  system.  However,  the  computational  procedures  remain 
comparatively  simple. 
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Figure  3-54  Maximum  deflection  of  elasto-pl astic ,  one-degree-of-f reedom  system  for 
triangular  load 
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Maximum  deflection  of  elasto-pl astic ,  one-degree-of-freedom  system 
for  rectangular  load 
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Maximum  response  time  of  elasto-plastic,  one-degree-of-freedom  for 
rectangular  load 


t/tn 

<12  as  I  2  5  10  20  40 


2 


5 


10 


20 


j  r,,/P  «  1.02 

Sj 


UJ  2 
X 


USCL  ' 


ru/P* 


Grodualiy  applied 
load 


Resistance 

function 


*m 

Displacement 

function 


L-u. 

I 


0.3  0.8  I 


Tr/TN 


8  10 


Figure  3-58  Maximum  deflection  of  elasto-plastic,  one-degree-of-freedom  system 
for  gradually  applied  load 


Fil 


A  bilinear  load  Is  Illustrated  in  figure  3-62c.  Four  parameters  are  required 
to  define  this  load  shape,  namely: 

P  peak  pressure  of  shock  load  (primary  pulse) 

CjP  =  peak  pressure  of  gas  load  (secondary  pulse) 

T  =  duration  of  shock  load 

C2T  =  duration  of  gas  load 

It  can  be  seen  from  figure  3-62c  that  the  value  of  Cj^  will  always  be  less  than 
1  while  the  value  of  C2  will  always  be  greater  than  1. 


Response  charts  for  the  bilinear  load  are  prepared  in  the  same  way  as  the 
simplified  loads  except  that  each  chart  is  prepared  for  a  given  value  of  and 
C2.  The  ductility  ratio  X_/X£  versus  T/Tu  is  plotted  on  each  chart  for  various 
values  of  P/ru,  t^/T^  and  tE/T  (where  tr  is  the  time  to  reach  the  maximum 
elastic  deflection  X^).  Therefore,  each  chart  contains  three  families  of 
curves.  In  addition,  each  chart  contains  one,  two  or  three  boundaries  which 
are  formed  by  symbols.  These  boundaries  delineate  regions  of  each  chart  where 
certain  approximations  are  preferable  to  chart  interpolation.  These 
approximations  Involve  modification  of  the  bilinear  load  for  the  various 
regions.  The  loads  considered  for  each  of  the  regions  are  given  in  figure  3-62 
while  the  boundary  for  the  various  regions  are  defined  in  figure  3-63. 


Numerous  response  charts  are  required  for  bilinear  loads.  Table  3-15  lists  the 
figure  numbers  of  the  response  charts  prepared  for  the  selected  values  of 
and  C2.  These  response  charts  are  presented  in  figures  3-64  through  3-266. 


To  use  the  response  charts,  the  first  step  is  to  enter  table  3-15  with  the 
required  values  of  and  C2.  Locate  the  points  in  the  table  with  coordinates 
Ci  and  C2.  The  number  in  the  box  containing  the  point  is  the  appropriate 
figure  number  to  use  to  solve  the  problem.  If  the  point  is  not  located  in  a 
box  having  a  number,  the  two  or,  more  frequently,  four  number  bracketing  the 
point  are  the  appropriate  figure  numbers  to  use.  Interoolation  for  the 
required  values  of  Ci  and  C2  must  be  performed.  A  graphical  and  mathematical 
interpolation  proceudre  is  presented  in  subsequent  sections. 


In  a  typical  design,  a  structural  member  is  assumed  and  the  idealized 

resistance  function  defined  by  ru,  X^-  and  can  be  determined  along  with  the 
natural  period.  As  previously  explained  for  simplified  loads;  the  equivalent 
mass  used  to  calculated  TN  varies  according  to  the  amount  of  plastic 

deformation.  Knowing  the  ratio  of  P/ru  and  T/T^,  note  if  the  intersection 
point  corresponding  to  these  parameters,  is  located  in  region  A,8,C  or  D  on  the 
response  chart(s)  for  the  required  values  of  and  C2.  The  boundaries  for 
these  regions  are  represented  by  a  line  of  asterisks,  solid  circles,  and  solid 
squares,  as  illustrated  in  figure  3-63. 


The  final  step  in  obtaining  the  solution  depends  on  which  region  of  the 
response  chart(s)  is  applicable.  These  regions  were  established  to  reduce  the 
amount  of  calculations  required  for  a  solution  by  eliminating,  where  possible, 
interpolation  between  charts.  Except  for  one  of  the  four  regions,  the  actual 
bilinear  blast  load  is  replaced  by  a  simplified  load.  Solution  of  the  response 
for  these  simplified  loads  will  predict  a  solution  which  is  within  10  percent 
of  the  exact  solution.  Of  course,  the  exact  solution  of  the  response  for  the 
structural  member  in  these  three  regions  can  be  obtained  by  using  the  required 
response  charts  and  interpolating,  where  required,  between  charts. 
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REGION  A  -  NEGLECT  CHARTS 
USE  EQUATIONS  3-80  AND  3-81 
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Figure  3-62  Various  bilinear  triangular  loads 
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3-19.3.4.1  Region  A  .  This  region  is  defined  as  the  area  to  the  left  of  the 
solid  circles  on  the  response  charts.  If  a  chart  does  not  have  a  line  of  solid 
circles,  then  region  A  does  not  exist.  Figures  3-221,  3-222  and  3-223 
Illustrate  this  case. 

In  region  A  the  maximum  dynamic  response  depends  primarily  upon  the  total 
impulse  in  the  blast  loading  (area  under  the  pressure-time  curve).  The  actual 
pressure-time  distribution  does  not  significantly  affect  the  maximum  dynamic 
response  because  in  these  regions  the  durations  T  and  C2T  are  small  in 
comparison  to  the  response  time,  ^  of  the  structural  member.  In  this  region 
the  response  charts  are  not  used  and  a  solution  is  obtained  by  considering  the 
modified  loading  shown  in  figure  3-62a.  This  load  shape  yields  the  following 
equations: 


1  +  i? 

7  T~ 


if  fe)(U 


where 


C,  +  C,(l-C,)(l-C9) 

1  (c2  -  q) 


3-19.3.4.2  Region  B  .  This  region  is  defined  as  the  area  in  the  chart  to  the 
left  of  the  asterisks.  If  a  chart  has  no  line  of  asterisks,  then  region  B  does 
not  exist.  Figure  3-65  illustrates  the  case  where  region  B  does  not  exist. 

In  region  B  the  maximum  dynamic  response  depends  primarily  upon  the  gas  load 
which  is  described  by  CjP  and  C2T.  The  shock  load  described  by  P  and  T  is 
neglected.  This  load  condition  is  shown  in  figure  3-62b.  Therefore  the 
solution  is  obtained  from  consideration  of  a  single  triangular  load.  The 
dynamic  response  is  obtained  from  figure  3-64  which  like  figures  3-54  and  3-55 
is  for  a  triangular  load  and  will  yield  the  same  results.  When  using  figure  3- 
64  it  must  be  realized  that  the  P  and  T  used  in  the  chart  are  actually  CjP  and 
C2T,  respectively.  In  a  typical  design,  enter  figure  3-64  with  the  normalized 
parameters,  CjP/ru  and  C2T/TN,  and  read  the  solution,  VXE’  t^/T^j,  and  t^/C2T. 

Figure  3-63  depicts  region  B  as  the  area  between  the  line  of  solid  circles  and 
the  line  of  asterisks.  It  should  be  understood  that  the  solution  technique 
associated  with  region  B  (i.e.,  neglect  the  shock  load  and  use  figure  3-64) 
applies  everywhere  to  the  left  of  the  line  of  asterisks,  including  region  A. 
In  other  words,  two  solution  techniques  are  available  in  region  A. 

3-19.3.4.3  Region  D  .  This  region  is  defined  as  the  area  in  the  charts  to 
the  right  of  the  line  of  solid  squares.  If  a  chart  has  no  line  of  solid 
squares,  then  region  D  does  not  exist.  Figure  3-69  illustrates  the  case  where 
region  D  does  not  exist. 

In  region  D  the  maximum  dynamic  response  depends  primarily  upon  the  shock  load 
which  is  described  by  P  and  T.  The  gas  load  described  by  CiP  and  C2T  is 
neglected.  This  load  condition  is  shown  in  figure  3-62d.  Therefore,  the 
solution  is  obtained  from  consideration  of  a  single  triangular  load.  Similar 


vj 


\  ■  -1 


-  114  - 


to  region  B,  the  dynamic  response  Is  obtained  from  figure  3-64  which  is  for  a 
triangular  load.  Unlike  region  B,  the  parameters,  P  and  T  which  describe  the 
load  are  used  in  the  figure.  Therefore,  In  a  typical  design,  enter  figure  3-64 
with  the  normalized  parameters  P/ru  and  T/TN,  and  read  the  solution,  Xm/XE, 
and  tE/T. 

3-19.3.4.4  Region  C  .  This  is  the  region  In  the  charts  which  do  not  meet  the 
definitions  of  regions  A,  B  and  0.  In  most  charts,  region  C  is  the  area  to  the 
left  of  the  line  of  solid  squares  and  to  the  right  of  the  line  of  asterisks,  as 
Illustrated  in  figure  3-63.  Region  C  does  not  exist  in  some  charts,  such  as 
figure  3-97  which  have  over  lapping  regions  A  and  D. 

In  region  C,  both  the  shock  and  gas  load  must  be  considered.  Replacement  of 
the  actual  bilinear  load  with  a  simplified  load,  as  done  for  the  other  regions, 
will  yield  an  incorrect  solution.  Therefore,  in  this  region  the  response 
charts  must  be  used.  If  the  required  values  of  and/or  C2  do  not  correspond 
to  the  response  chart  values,  interpolation  between  response  charts  will  be 
required  to  obtain  a  solution.  In  a  few  cases  one  or  two  response  charts  are 
needed,  however.  In  general  four  response  charts  are  required  for  a  solution. 

3-19.3.4.5  Reeponae  Chart  Interpolation  .  Interpolation  between  four 
response  charts  will  usually  be  required  for  region  C.  For  regions  A,  B,  or  D, 
if  conditions  warrant  an  exact  solution .  rather  than  the  approximate  solution 
usually  used  in  these  regions,  interpolation  between  charts  may  be  required. 
Either  a  graphical  or  mathematical  interpolation  procedure  may  be  employed. 
The  method  selected  depends  upon  personal  choice.  A  brief  description  of  each 
procedure  is  presented  below  and  an  example  of  each  is  given  in  Appendix  A. 

Graphical  interpolation  requires  a  sheet  of  log-log  graph  paper.  A  convenient 
size  is  2  x  1  cycle  with  the  single-cycle  axis  representing  Cj  and  C2,  and  the 
two-cycle  axis  representing  the  desired  parameter  (Xm/XE,  t^/Tjj  or  tr/T), 
called  Y  for  ease  of  presentation.  The  appropriate  figures  to  be  used  for  a 
solution  are  obtained  form  table  3-15  for  the  required  values  of  and  C2< 
The  procedure  will  illustrate  the  interpolation  between  four  figures  since  this 
is  by  far  the  usual  case.  For  the  values  of  P/ru  and  T/T^  corresponding  to  the 
structural  system  selected,  determine  the  desired  parameter  Y  for  each  of  the 
four  figures.  Organize  a  table  in  the  same  format  as  table  3-16  and  enter  each 
figure  number  and  corresponding  value  of  C^,  C2  and  Y. 


Table  3-16  Response  Chart  Interpolation 


Figure  Number  C1  C2  Desired  Parameter 


In  the  table,  Cji  and  C21  are  the  values  of  and  C2 ,  respectively,  from 
figure  1.  Likewise,  Ci2  anc*  C?2  are  the  va^ues  ^1  an(*  ^2*  respectively, 
from  figure  2,  etc.  Tne  symbol  Yj,  Is  the  desired  parameter,  such  as  Xm/Xr, 
which  Is  read  from  figure  1.  The  symbol  Y?  is  the  value  read  from  figure  2, 
etc.  The  interpolation  Is  first  performed  for  the  required  value  of  and 
then  for  the  required  value  of  C2.  That  Is,  with  C2j  constant,  graphically 
interpolate  on  the  log-log  paper  between  points  (Yj,  Cii)  and  { Yo ,  Ci2)  to 
Yfl  corresponding  to  Cj  and  C2i,  as  shown  on  figure  3-2b7.  Similarly,  with  C23 
constant,  graphically  Interpolate  between  (Yo,  C13)  and  (Y4,  C14)  to  find  Yu 
corresponding  to  and  C23,  also  shown  in  figure  3-267 .  The  values  of  Ya  ana 
*b  are  recorded  In  table  3-16.  Finally,  with  C.  con^tjnt,  graphically 
interpolate  on  the  log-log  paper  between  (Ya,  C2^)  a^nd  (Yb,  C03)  to  find  Y 
corresponding  to  the  required  and  C2,  as  snown  on  figure  3-267.  The  values 
of  Y  Is  the  solution.  Since  there  are  three  parameters  to  define  the  response 
of  a  structural,  namely,  X_/XE,  t-,/TN  and  tE/T,  the  interpolation  procedure  is 
repeated  three  times,  once  for  eacn  parameter. 

Mathematical  interpolation  requires  the  same  initial  steps  as  graphical 
interpolation.  That  is,  the  appropriate  figures  to  be  used  for  a  solution  are 
obtained  from  table  3-15  and  the  required  parameters  are  determined  and  entered 
into  table  3-16.  Logarithmic  equations  are  used  to  obtain  the  values  of  Ya,  Yb 
and  Y.  The  value  of  lnYa  is  obtained  from: 


lnYa  = 


lnYx  + 


ln(Yo/Yi )1n(C, /Ci 1 ) 


1 n ( C 1 2/C  2 j) 

and  the  value  of  lnYb  is  obtained  from: 

lnYb  =  1 nYi  *  1 Y4/^3 ) 1 C1/C13 ) 


'14/U13J 


Using  the  values  of  lnYa  and  lnYb;  the  value  of  InY  is  obtained  from: 

(lnYu  -  lnY.)ln(Co/Coi ) 


lnYa  + 


1 n( C23/C2 1 ) 


The  desired  parameter  Y  is  then  obtained  from: 


3-86a 


The  above  equations  use  natural  logarithms.  Common  logarithms  can  be  used  to 
solve  the  above  equations  and  then  solve  for  Y  using: 

Y  =  101o9Y  3-86b 


It  should  be  noted  that  if  C2  is  represented  by  a  response  chart,  then  only  two 
charts  will  be  involved,  and  only  equation  3-83  will  apply. 


3-19.3.4.6  Accuracy .  The  prediction  error  is  less  than  10  percent  for  the 
approximate  solutions  obtained  in  regions  A,  B  and  D.  This  is  true  provided 
the  recommended  procedures  are  employed,  that  is  ,  the  actual  load  is  replaced 


by  the  approximate  loads  as  shown  on  figure  3-62  and  the  solution  is  obtained 
using  the  equations  provided  for  region  A  and  using  figure  3-64  for  regions  B 
and  D. 

Exact  solutions  are  obtained  from  the  response  charts  if  the  required  values  of 
Cj  and  C£  correspond  to  those  given  in  the  charts.  This  is  true  for  all  four 
regions.  Errors  result  from  interpolation  between  the  response  charts.  The 
prediction  error  in  Xm/XE  for  region  C,  where  interpolation  between  charts  is 
required,  is  less  than  10  percent  provided  response  charts  bounded  by  the 
dashed  or  solid  lines  or.  table  3-15  are  not  used.  The  prediction  error  in 
X_/X£  will  range  from  55  to  100  percent  for  solutions  in  region  C,  if  the  point 
(Cj,  Co)  lies  in  the  area  bounded  by  the  solid  line  in  table  3-15,  and  will 
range  from  10  to  55  percent  for  the  area  hounded  by  the  dashed  line.  If  the 
solution  involves  charts  from  both  sides  of  either  the  dashed  or  solid  lines, 
the  prediction  error  will  range  from  10  to  55  percent.  The  large  interpolation 
error  for  these  charts  result  from  the  big  change  in  Xm/X^  for  a  small  change 
in  P/ru  which  is  peculiar  to  this  chart. 

It  should  be  noted  that  the  large  errors  described  above  will  always  be  on  the 
high  side,  that  is  ,  the  predicted  value  of  Xm/X^  will  always  be  greater  than 
the  exact  value.  Hence,  a  conservative  design  will  always  result.  In 
addition,  it  should  be  noted  that  the  prediction  error  in  t^/T^j  is  about  half 
the  error  in  X(n/XE. 

While  the  errors  produced  from  interpolation  between  the  charts  in  table  3-15 
bounded  by  the  dashed  and  solid  lines  is  large,  the  area  in  which  they  apply  is 
small.  As  can  be  seen,  region  C  is  rather  small  on  these  charts.  It  is  for 
points  in  this  small  area  where  interpolation  must  be  performed.  For  the 
remaining  areas  (A,  B  and  D),  the  approximate  solutions  may  be  used  rather  than 
interpolation.  These  approximate  solutions  will  result  in  a  error  of  less  than 
10  percent. 

Solutions  involving  0.920  <  Cj  <  1.00  and  C2  >_100  can  result  in  very  large 
errors  if  interpolating  procedures  are  employed.  In  this  region,  the  dynamic 
response  is  primarily  due  to  the  gas  load.  Therefore,  figure  3-64  is  used  to 
obtain  X^/X^-  for  the  gas  load.  However,  this  value  is  too  low  and  should  be 
increased  by  20  percent.  The  value  of  tm/Tij  corresponding  to  the  increased 
value  of  Xm/XE  is  then  obtained  from  figure  3-64 . 

3-19.3.5  Determination  of  Rebound  .  In  the  design  of  elements  which  respond  to 
the  pressure  only  and  pressure-time  relationship,  the  element  must  be  designed 
to  resist  the  negative  deflection  or  rebound  which  can  occur  after  maximum 
positive  deflection  has  been  reached.  The  ratio  of  the  required  unit  rebound 
resistance  to  the  ultimate  unit  resistance  r/ru,  such  that  the  element  will 
remain  elastic  during  rebound  is  presented  in  figure  3-268  for  a  single-degree- 
of-freedom  system  subjected  to  a  triangular  loading  function.  Entering  with 
ratios  Xm/XE  and  t/T^  previously  determined  for  design,  the  required  unit 
rebound  resistance  r-  can  be  read  in  terms  of  the  originally  designed  ultimate 
unit  resistance  ry.  To  obtain  the  rebound  resistance  for  an  element  subjected 
to  another  form  of  load  other  than  the  triangular  loading  function,  a  time- 
history  analysis  such  as  the  one  described  in  Section  3-19  has  to  be  performed. 
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LOAD  DURATION  TO  FUNDAMENTAL  PERIOD.  T/TN 

Figure  3-72  Maximum  response  of  el asto-pl asti c ,  one-degree-of-f reedom 
system  for  bilinear-triangular  pulse  (C^  =  0.032,  C?  =  1. 
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LOAD  DURATION  TO  FUNDAMENTAL  PERIOD.  T/TN 

Figure  3-78  Maximum  response  of  elasto-plastic ,  one-degree-of-freedom 
system  for  bilinear-triangular  pulse  (Ci  =  1.215,  C?  =  3.1 


j/37  .'Nouvana 


riON  ro  rUNDAMCNFAL  PCRIOD.  T/TN 

response  of  elasto-plastic ,  one-degree-of-f reedom 


YICLD  TO  LOAD  DURATION,  /j-/T 


Ni/UJ;  .'cioiaua  ~iviN3nvaNnj  oi  Noiiorunci  nnwixvw  jo  uwu 


DURATION  TO  TUNDAMCNTAL  PERIOD,  T/TN 

mum  response  of  elasto-plastic,  one-degree-of-freedom 


Nouvana 


iiKSSi 


EwisainsisaiSiaiiKsawesitssiSbSiKsaeNiiBH 

»^«»!Sisp^s!ea»ii8s»eHaaSiaa 


LOAD  DURATION  TO  FUNDAMENTAL  PERIOD,  T/TN 

Maximum  response  of  elasto-plastic ,  one-degree-of-f reedom 
system  for  bilinear-triangular  pulse  (Cj  =  0.056,  C2  =  3. 


riMC  Of  YILLD  10  LOAD  DURATION.  it/T 


Ni/UJ?  ‘aoiaid  iviNiwvaNnj  oi  noiio33j30  wnwixvw  jo  3wu 


•*“>  —  O  — 


o 

E 

• 

o 

rn 

T3 

OJ 

II 

CD 

V. 

CM 

LJ 

V- 

41 

o 

CM 

1 

m 

a>  o 

a> 

i 

o 

C7> 

a> 

ii 

z 

T3 

1 

r— ♦ 

\ 

U 

o 

►— 

c 

o 

« 

<D 

Q 

4k 

0 

O 

u 

‘f— 

3 

cc 

4-J 

a. 

u 

CO 

CL 

OJ 

i. 

03 

CL 

< 

1 

3 

k— 

O 

O' 

Z 

w 

c 

u 

CO 

OJ 

2 

C_ 

< 

4-i 

Q 

1 

Z 

U 

3 

t_ 

o 

a» 

<D 

c 

O 

CO 

k— 

c 

O 

•4" 

z 

CL  A 

o 

CO 

Of 

c. 

r*_ 

o 

< 

CL 

& 

3 

a 

E 

Q 

E 

'H 

4-> 

o 

X 

CO 

-tj 

>> 

z: 

»o 

o 

-J 

<\J 

03 


a» 

u 


-137- 


TIM C  or  YICLD  TO  LOAD  DURATION,  ^/T 


l/w1  'aoia^d  TviNiwvaNnj  01  NOiiorLoa  nhwixvin  jo  jnl 


leseazwwimm 

musssam 


SalSS§i«?:?aS»K^SB5«s!!?ii 

HHBMaVIMWReSf 


DURATION  TO  FUNDAMENTAL  PERIOD,  T/TN 

imum  response  of  elasto-plastic ,  one-de<iree-of-f  reedcm 


TIME  Of  nLLD  TO  LOAD  DURATION.  lt/T 


Ni/UJ;  'aoiaod  iviNiwvaNrij  01  noiiodudo  wnwxvw  jo  3wll 


—  o  —  "1 

%>  ^  >  fe >  «J-  W  V  V  *5  v  -v  »v  ~  -  O  O  O  O  O 


luiisw^K^saraa^ii^ssMRiiiiiiKViS'SiiHi^iriiiaBH 

I^BSiS9S^;3i[BV^Sil^i'«4R^|llllll 


^■nkxpuimjv< 


mkv.aviii  ■■ 


l^ansaaaa ’timEUMav’*wm>^r^2r^is'»ziizfwnFye,mxmws*m-mKirwnMMwm  1 
|^r«mJK.Wk?9iMK>5HAr^aai»'^KV?»rA'tv;aikiaH] 

hi  91 

■8ag«esaa>M— ■■■■■■— ■— ■— i 


aessRSsainissjissii 


HflBHNfMHH 


BMHB 


■K:aiiesfli!irii!iiR:tf!Mi&aiiii2ii^!U!K«9Vjis«w) 

iSiSSSSiB  I!  SSI*i?^ifi!i!SC!!5i  iJISSSS.SMiiSWfi 

»m^^0igsig0^8i 

Sie;&!lg:3fll^!^g:3Bg»e:!ilg!:i&Si8RS!i!8SB 


o  o  o  o  _ 

q  p*.  .n  -O 


3x/UJx  nouo 3uiG  ai3iA  oi  NOiiDiijia  anwixvn 


%  \  \  '.  '. 
•/’••  V-'  vVW-. 
,  .*,  -\  **.  .*.  /.  / 

-  ’  V  '.L  ■  -  ■  .  ’  .  -  .'■ 


LOAD  DURATION  TO  FUNDAMENTAL  PERIOD.  T/TN 

Figure  3-84  Maximum  response  of  elasto-plastic ,  one-degree-of-f reedom 
system  for  bilinear-triangular  pulse  (Cj  =  0.010,  C 2  *  3.1 
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LOAD  DURATION  TO  FUNDAMENTAL  PERIOD,  T/TN 

Figure  3-100  Maximum  response  of  elasto-plastic,  one-degree-of-f reedom 
system  for  bilinear-triangular  pulse  (Cj  =  0.648,  C2  =  10) 
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Maximum  response  of  el asto-pl asti c ,  one-degree-of-f reedom 
system  for  bilinear-triangular  pulse  (C\  =  0.215,  C">  =  30 


Ni/LU;  'aoiaid  3viN3wvaNnj  01  noiio33j3cj  wnwixvw  jo  3wn 


r\ 


%T  VV  to1'  *5 


•*.  «  o  <o 

v  v  *5'  *>'  V  'V  *V  — 


ni  —  o  —  ~i 

0)  CQ  rs  (A  if) 

oodo  o 


_ IIDM1! _ 


<  0> 
“  0V> 


pcsKsSSSBHSsSsSSfiSiiMS^BtHSSSiHSBHiEBiBB 


|isw*SAi»»aKMiv!aisiiK^KMK:w»Wi»Mna 

S^iir^l!5SSii&lGK:^ 

ss^^iii!s^§5^ssii^»sa 

lg|gpiSi||giii|||i|| 

li?^SteSSiiga5^3!ly>5s”ga£sSS^»ti8| 


Unn 


IKS«I 


4o  a 
o  o 


o  o  o  o  o 

o  *»  W>  *0  IN 


__  rv  .n  o 


'//^X  'NOUD3U3G  G-1]  a  01  N0U033J3C  Nn/HXViN 


-184- 


002  003  00b  oi  02  03  05  .07  .1  .2  3  .5  .7  I  2  3 

LOAD  DURATION  TO  FUNDAMENTAL  PERIOD,  T/TN 

Figure  3-129  Maximum  response  of  elasto-plastic ,  one-degree-of-f reedom 
system  for  bilinear-triangular  pulse  (C]  =  0.178,  C?  =  30 
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LOAD  DURATION  TO  FUNDAMENTAL  PCRIOD,  T/T^ 

Figure  3-192  Maximum  response  of  elasto-plastic ,  one-degree-of-freedom 
system  for  bilinear-triangular  pulse  (C^  =  0.316,  =  300. 
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LOAD  DURATION  TO  FUNDAMENTAL  PERIOD,  T/TM 

Figure  3-200  Maximum  response  of  elasto-plastic,  one-degree-of-freedom 

system  for  bilinear-triangular  pulse  (Cj  =  0.162,  C2  =  300.) 
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LOAD  DURATION  TO  FUNDAMENTAL  PERIOD,  T  /  T  N 

Figure  3-203  Maximum  response  of  elasto- plastic ,  one-degree-of-f reedom 

system  for  bilinear-triangular  pulse  (C^  -  0.121,  C2  =  300.) 
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LOAD  DURATION  TO  FUNDAMENT AL  PCRIOD,  T/TN 

Figure  3- 2D 7  Maximum  response  of  elasto-plastic  ,  one-degree-ot-f reedoin 

system  fur  bilinear-triangular  pulse  (C^  -  0.083,  C?  =  300. 
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LOAD  DURATION  TO  FUNDAMENTAL  PERIOD,  T/TN 

Figure  3-2 22  Maximum  response  of  elasto-plastic ,  one-degree-of-f reedom 

system  for  bi 1 inear-tri angul ar  pulse  ( C j  =  0.825,  C2  =  1000.) 
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LOAD  DURATION  TO  FUNDAMENT AL  PERIOD,  T/T,g 

Figure  3-228  Maximum  response  of  el asto-pl asti c ,  one-degree-of-f reedom 

system  for  bi 1 l near- tri angul ar  pulse  (Ci  =  0.819,  Co  =  1000. 
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LOAD  DURATION  TO  FUNDAMENTAL  PERIOD,  T/ TN 

Figure  3-245  Maximum  response  of  el  a s to-pl ast i c ,  one-degr^e-of-f reedom 

system  for  bilinear-triangular  pulse  (C^  =  0.198,  C2  =  1000. 
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LOAD  DURATION  TO  FUND AMCNT AL  PCRIOD,  T/TN 

Figure  3-2'57  Maximum  response  of  elasto-plastic ,  one-degree-of-f reedom 

system  for  b  1  1  1  nea  r- 1.  r  i  angul  ar  pulse  (Cj  =  0.056,  =  1000. 
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LCaD  duration  to  ruNO amcnt al  period.  t/tn 

Figure  3-359  Maximum  response  of  elasto-plastic ,  one-degree-of-f reedom 

system  for  bilinear-triangular  pulse  (Ci  =  0.042,  C2  -  1000.) 
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LOAD  DURATION  TO  FUNDAMENTAL  PERIOD,  T/TN 

Figure  3-263  Maximum  response  of  el  a sto-pl ast i c ,  one-degree-of-f reedom 

system  fur  bilinear-triangular  pulse  (Cj  =  0.018,  C2  =  1000. 
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LOAD  DURATION  TO  FUNDAMENTAL  PERIOD,  T/TN 

Figure  3-264  Maximum  response  of  elasto-plastic ,  one-degree-of-freedom 

system  for  bilinear-triangular  pulse  (Cj  =  0.015,  C?  =  1000.) 
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Figure  3-268  Elastic  rebound  of  simple  spring-mass  system 


It  may  be  noted  that  if  the  loading  is  applied  in  a  relatively  short  time 
compared  to  the  natural  period  of  vibration  of  the  system,  the  required  rebound 
resistance  can  be  equal  to  the  resistance  in  the  initial  design  direction. 
When  the  loading  is  applied  for  a  relatively  long  time,  the  maximum  deflection 
is  reached  when  the  positive  forces  are  still  large  and  the  rebound  resistance 
is  reduced. 

3-20  Elements  Which  Respond  to  Inpluse 
3_20.1  General  Equations  For  Maximum  Response 

When  an  element  responds  to  the  impulse,  the  maximum  response  depends  upon  the 
area  under  the  pressure  time  curve  (impulse  of  blast  loading).  The  magnitude 
and  time  variation  of  the  pressure  are  not  important.  The  response  charts 
presented  in  Section  3-19,  which  are  based  on  pressure-time  relationship  are 
therefore  not  required  for  these  problems.  Instead,  the  element  resistance 
required  to  limit  the  maximum  deflection  to  a  specific  value  is  obtained 
through  the  use  of  a  semi  graphical  method  of  analysis. 


Consider  the  pressure-time 


resistance-time  functions  shown  in  figure 


3-269.  The  resistance  curve  dep’ctr 
resistance-deflection  function  having 
equation  of  motion  it  can  be  shown  tha* 
area  A  as  positive  and  area  B  as  negati.:- 
time  ta  divided  by  the  corresponding 
instantaneous  velocity  of  that  time: 


s  for  a  two-way  element  with  a 
-t-ultimate  range.  From  Newton’s 
•'  ^nation  of  the  areas  (considering 
\  r  the  load-tir..e  curves  up  to  any 
rive  masses  is  equal  to  the 


■  ( 5a(f^i 

u  m. 


The  displacement  at  time  t,  is  found  by  multiplying  each  differential  area 
divided  by  the  appropriate  effective  mass  by  its  distance  to  ta  and  summing  the 
values  algebraically: 


xa  =  (nd(f-r)  (ta-t)dt 

(J  m„ 


In  each  range  the  mass  is  the  effective  plastic  mass: 


Time  Interval 


0<  t<t{ 


Resistance 


Effective  Mass 


mu  ~  ^KLM^um 
mup  =  (KLM^upm 


Time  ti  is  the  time  at  which  the  partial  failure  deflection  X^  occurs,  and  time 
tj,,  is  the  time  at  which  maximum  deflection  Xm  is  reached  (Xm  <  Xu). 

For  an  element  to  be  in  equilibrium  at  its  maximum  deflection,  its  impulse 
capacity  must  be  numerically  equal  to  the  impulse  of  the  applied  blast  load. 
With  the  use  of  the  foregoing  equations,  the  expressions  which  define  the 
motion  and  capacity  of  elements  subjected  to  impulse  type  loads,  can  be 
defined.  These  expressions  are  presented  for  both  large  and  limited  deflection 
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criteria.  This  criteria  varies  for  the  different  materials  used  in  protective 
design.  The  criteria  for  each  material  is  obtained  from  the  volume  that 
describes  the  design  procedures  for  that  material. 


Case  1  -  larqe  deflections.  Utilizing  equation  3-38  and  by  taking  moments  of 
the  areas  under  the  pressure-time  and  resistance-time  curves  (fig.  3-_39) 


about 


assumi ng 


blast  impulse 


instantaneously  at  time  t=0,  and  that  time  to  reach  yield  t 
zero  the  expression  for  the  maximum  deflection  is 


is  applied 
also  close  to 


Xm  *  Jo. 

m.. 


rupl*m~^l I 


If  moments  of  the  areas  are  taken  about  tj ,  then  the  deflection  at  partial 
failure  Xj ,  is 

X1  =  ihti  -  r))t  i  3-80 

m..  2m,, 


Using  equation  3-87  and  summing  the  areas  of  t^  and  recognizinq  that  the 
instantaneous  velocity  at  t^  equals  zero. 

Ik  -  r„ t,  -  rllf,(  t_  -  t, )  =0  3-91 


Solution  of  the  above  three  simultaneous  equations  is  accomplished  by  solving 
equation  3-91  for  t_  and  substituting  this  expression  into  equation  3-89,  and 
solving  eauation  3-^0  for  t^  and  substituting  this  expression  into  the  modified 
equation  3-89.  After  continuing  and  rearranging  terms,  the  general  response 
equation  becomes 


Sip  ,Xm  *  xl> 


The  left  side  of  this  equation  is  simply  the  initial  kinetic  energy  resulting 
from  the  applied  blast  impulse  and  the  right  side  is  the  modified  potential 
energy  of  the  element.  The  modification  is  required  since  the  above  analysis 
requires  the  use  of  two  equivalent  dynamic  systems  (before  and  after  time 
ti).  The  modification  factor  mu/m  equates  the  two  dynamic  systems.  If  the 
effective  mass  in  each  range  was  thepsame,  mu/m  would  equal  one  and  the  right 
side  of  the  expression  would  be  ruXm  which  is  the  potential  energy. 


For  one-way  elements  which  do  not  exhibit  the  post-ultimate  resistance  range, 
or  for  two-way  panels  where  the  maximum  deflection  Xm  is  less  than  Xj,  equation 
3-92  becomes 


The  above  solutions  are  valid  only  for  wnat  is  considered  large  deflection 
design  since  the  variation  of  resistance  with  deflection  in  the  elasto-plastic 
range  has  been  ignored.  This  limitation  1  s  based  on  the  assumption  that  the 


time  to  reach  yield 
comparison  to  tm. 


duration  of  the  impulse 


are  small 
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Case  2  -  limited  deflections.  For  elements  which  respond  to  the  impulse  with 
limited  deflections  where  tt  a  time  to  reach  maximum  deflection  tm  is  greater 
than  three  times  the  duration  t0  of  the  load,  but  where  the  support  rotations 
are  equal  to  or  less  than  the  established  criteria,  the  elasto-plastic  range 
behavior  of  the  element  must  be  accounted  for  in  determining  the  overall 
response  of  the  element  to  the  applied  blast  load.  For  this  case,  equation  3- 
92  becomes 

jh2  3  A  +  fa.  ru(xm  -  XE)  3-94 

~  mu 

where  ma  is  the  average  of  the  effective  elastic  and  plastic  masses  and  XE  the 
equivalent  deflection. 

When  the  response  time  t^  of  an  element  is  less  than  three  times  the  load 
duration  tQ,  the  element  will  respond  to  the  pressure  pulse  rather  than  to  the 
impulse  alone.  In  this  case  the  response  of  the  element  may  be  obtained 
through  the  use  of  a  response  chart  considering  a  fictitious  pressure  pulse  as 
outlined  in  Volume  II.  However,  if  the  element's  maximum  deflection  is  greater 
than  given  in  the  response  charts  (ductility  ration  Xm/X^  greater  than  100), 
the  response  of  the  element  may  be  obtained  through  use  of  the  semi  graphical 
method  of  analysis  as  outlined  in  this  paragraph  considering  the  fictitious 
pressure  pulse.  It  should  be  noted  that  the  resistance  time  relationship  used 
in  the  analysis  to  express  the  element's  response  should  include  the  elasto- 
plastic  region. 

3-20.2  Determination  of  Time  to  Reach  Maximum  Deflection 

The  time  t^  for  each  of  the  cases  covered  previously  can  be  determined  by 
applying  equations  3-87  and  3-88  to  the  particular  problem.  The  resulting 
equations  are  as  follows: 

Case  1  -  large  deflections. 


3-95 


3-96 


Since  the  variation  of  resistance  with  time  is  not  known  in  the  elasto-plastic 
range,  t^  can  only  be  determined  approximately  by  assuming  a  linear  variation 
of  resistance  with  time. 
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APPENDIX  3A  ILLUSTRATIVE  EXAMPLES 


Probl em: 

Procedure: 
Step  1. 
Step  2. 

Step  3. 

Step  4. 

Step  5. 

Step  6. 
Step  7 

Procedure: 
Step  I. 
Step  2. 
Step  3. 


Problem  3A— 1(A)  Ultimate  Unit  Resistance 


Determine  the  ultimate  unit  resistance  of  a  two-way  structural 
element  using  (1)  general  solution  and  (2)  charts. 

Part  (a)  -  General  Solution 

Establish  design  parameters. 

Assume  yield  line  locations  in  terms  of  x  and/or  y  considering 
support  conditions,  presence  of  openings,  etc. 

Determine  negative  and  positive  moment  capacities  of  sections 
crossed  by  assumed  yield  lines.. 

Establish  distribution  of  moments  across  negative  and  assumed 
yield  lines,  considering  corner  effects  and  those  of  openings. 

Determine  the  ultimate  uni t  .resistance  for  each  sector  in  terms 
of  x  and/or  y  considering  free  body  diagram  of  the  sectors 
(fig.  3-3).  Summation  of  the  moments  about  the  axis  of  rotation 
(support)  of  the  sector  yields  equation  3-3. 

Equate  the  ultimate  unit  resistance  of  the  sectors  and  solve  for 
the  yield  line  location  x  and/or  y. 

With  known  yield  line  location,  solve  for  ultimate  unit 
resistance  of  the  element,  using  equations  obtained  in  Step  5. 

Note:  For  complex  problems  (three  or  more  different  sectors), 

the  solution  for  the  ultimate  unit  resistance  is  most  easily 
accomplished  through  a  tri al -and-error  procedure  by  determining 
r  for  each  sector  for  a  given  (assumed)  yield  line  location  and 
adjusting  the  yield  lines  until  the  several  values  of  ru  agree 
to  within  a  few  percent. 

Part  (b)  -  Chart  Solution 

Same  as  in  step  1  of  part  a. 

Same  as  in  step  2  of  part  a. 

Determine  the  negative  and  positive  ultimate  moment  caoacities 
in  vertical  and  horizontal  directions. 
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Step  4.  For  given  support  conditions  (and  value  of  in  the  case  of 

an  element  with  three  edges  supported  and  fourth  free),  use  the 
appropriate  chart  (figs.  3-4  through  3-20)  to  obtain  yield  line 
location  ratios  x/L  or  y/H  for  value  of  quantity  abtained  in 
step  4.  Then  calculate  x  or  y. 

Step  5.  Using  the  appropriate  equation  from  table  3-2,  determine  the 
ultimate  unit  resistance  of  the  element. 


Example  3A-1  Ultimate  Resistance 


Required:  Ultimate  unit  resistance  of  two-way  structural  steel  element 

shown  below  using  (1)  general  solution  and  (2)  charts. 


Figure  3A  -  1 

(Solution:  Part  (a)  -  General  Solution) 

Step  1.  Given: 

(a)  L  =  240  in  H  =  163  in 

(b)  Fixed  on  three  sides  and  free  at  the  fourth 
Step  2.  Assume  yield  line  location  (fig.  3A  -  2) 

Step  3.  The  negative  and  positive  moment  capacities  in  both  the 
horizontal  and  vertical  directions  are  determined  from  the 
properties  of  the  material.  For  this  example,  it  will  be 
assumed  that  the  moment  capacities  are  equal  to  M  =  20,000  in¬ 
lbs/in. 
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168 


Step  4. 
Step  5 


Step  6. 


Mhn  =  MHp  =  Mvn  =  My p  =  20,000  in  Ibs/in. 

For  distribution  of  moments  across  negative  and  assumed  positive 
lines,  see  figure  3A  -  2(a). 

The  ultimate  unit  resistance  of  each  sector  is  obtained  by 
taking  the  summation  of  the  moments  about  its  axis  of  rotation 
(supports)  so  that 

ZMn  +  sMp  =  Rc  =  ruAc 

a.  Sector  I  (fig.  5A  -  2) 

IMVN  +  ZMVP  =  120(20,000)  +  2(2/3>(20 ,000 ) ( 60 )  + 

'  120(20,000)  +  2(2/3)  (20,000)(60) 

=  8.0x10^  in-lbs. 


ruAc  =  ru 


240M.J  f 


J  ]  =  40  ruy2 


therefore, 

ru  =  400(20,000)/40y2  =  0.2xl06/y2 
b.  Sector  II  (fig.  3A  -  2) 

EMuki  +  EMHp  =  (168-y/.2)  (20,000)  +  2/3(20, 000)(y/2)  + 

( 20 , 000 ) ( 168-y/2)  +  2/3(20,000)  (y/2 ) 

=  336(20,000)  -  y/3(20,Q00) 

ruAc  .  ^  [12Q(168+168-y)]  [lZ0[168^(168-y)]]/,1MM68.y) 

.  4,800rg(252-y) 
therefore, 

336(20,000)  -  y /3 ( 20 , 000 ) 
ru  4800  ( 252-y ) 

Equate  the  ultimate  unit  resistance  of  the  sectors. 

10(20,000)  =  336(20,000)  -  y/3( 20 ,000) 

- ^ -  - 4ff00r252“y-) - 

Simpl i fyi ng: 

y3  -  1008y2  -  144000y  +  36288000  =  0 
and  the  desired  root  is:  y  =  137.6  ins. 
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Step  7.  The  ultimate  unit  resistance  is  obtained  by  substituting  the 
value  of  y  into  either  equation  obtained  in  step  5,  both  of 
which  yiel  d: 

ru  =  (20,000)/ ( 137 . 6)2  =  10.6  psi 

Solution:  Part  (b)  -  Chart  Solution 
Note: 

Element  conforms  to  the  requirements  of  section  3-8  since  it  is  fixed  on 
three  sides  and  free  on  the  remaining  side  and  has  uniform  thickness  in 
the  horizontal  and  vertical  directions. 

Step  1.  Same  as  step  1  in  part  a. 

Step  2.  For  illustrative  purposes,  a  different  yield  pattern  (fig.  3A-3) 
will  be  assumed. 


Figure  3A  -  3 

Step  3.  For  ultimate  moment  capacities,  see  step  3  of  part  a 

Step  4.  For  three  sides  fixed  and  the  fourth  free,  calculate  the 
parameter . 

x2/x1=[(MHn3  +  MHp)/(MHM1  +  MHp ) l1/2  =  [(20,000) (2)7(20,000) ( 2 ) ] 1 /2  =  1.0 
From  figure  3-11,  (X2/X^  =  1.0)  calculate  the  parameters: 

L/H[Mvp/(MuW1  +  Mup  )  ]  =  240/168C20, 000/(2) (20, 000)]1/2  =  1.01 


and 


MVP  20,000 

- =  — -  =  1.0 

^VN2  20,000 

Read  yield  line  location 

X^/L  exceeds  the  maximum  possible  value  of  0.5  therefore,  assumed  yield 
line  pattern  is  wrong.  Assume  alternate  yield  line  pattern  as  shown  in 
figure  3A-2. 

From  figure  3-16  calculate  the  following  parameters: 

L_  _ (Mvn3  +  Mvp)1/2 _ 

H  =  (MHN2+MHP)1/2  +  ^HN1+MHP)1/2 

240  _  (20, 000+20, 000)1/2 

= - C - Tn - - FT?]  =  0.71 

168  (20, 000+20, 000)i/2  +  (20 ,000+20 ,000)1/2 

and 

X/L  =  [(Mhn1  *  MHp)/(MHN2  ♦  MHp>]1/2  /  1+C(Mhn1  ♦  «hp)/(Mhn2  *  M„p)]1/2 
=  [40, 000/40, 000]1/2  /  l+[40, 000/40, 000]1/2  =  1/2 
From  figure  3-16  read  of  yield  line  location: 
y/H  =  0.82;  y  =  0.8(168)  =  137.6  in 
X/L  =  0.50;  X  =  0.5(240)  =  120.0  in 
Step  5.  From  table  3-2 

NOTE:  Both  equations  given  in  the  table  for  each  edge  condition  and 

yield  line  location,  will  provide  identical  values  of  ru. 

5(MVN+MVp)  5(20,000+20,000) 

r  =  - o - =  - P - =  (0.00055)20,000  =  10.5  psi 

u  yZ  137. 62 


Example  3A-1(B)  Ultimate  Unit  Resistance 

Required:  Ultimate  unit  resistance  of  the  element  considered  in  example 

3 A- 1 ( A )  except  there  is  an  opening  as  shown  in  figure  3A-4. 
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Figure  3A  -  4 


Solution: 

Step  1.  Given: 

(a)  L  =  240"  H  =  168"  Two  additional  free  edges  are  formed 
due  to  the  presence  of  the  opening. 

Step  2.  Assumed  yield  line  location,  is  shown  in  figure  3A-5  (three 
different  sectors  are  formed). 


Step  3.  Same  as  step  3  of  Example  3A-KA),  part  a. 


Step  4.  For  distribution  of  moments  across  negative  and  assumed 
positively  yield  lines  see  figure  3 A- 5 .  (Since  opening  is 

located  at  lower  right  corner,  there  is  no  reduced  moment 
capacity  in  this  area.) 

Step  5.  The  ultimate  unit  resistance  is  obtained  from: 

IMjj  +  £Mp  -  Rc  -  ruAc 

a.  Sector  1  (fig.  3A-5) 

ZMHN+  £ Mup  =  (20 ,000 ) ( 168-y /2 )  +  2/3(20,000) (y/2)  + 

(20,000)  ( 168-y /2 )  +  2/3(20, 000) (y/2) 

=  336(20,000)  -  y/3{20,000) 

=  (336-y/3)  (20,000) 

ruAc  =  ru[x( 168+168-y)/2][x( 168+2  ( 168-y ) )/3(168+168-y)] 

=  rux2(252-y)/3 

therefore, 

(1008-y)(20,000) 

^  ‘  x2(252-y) 


b.  Sector  II  (fig.  5A-5) 

£MHK|  +  £Mup  =  (20,000)(80)  +  (20,000)  (80) 
=  160(20,000) 

Note: 


The  sector  is  divided  into  parts  a,  b,  and  c  so  that  the 
centriod  may  be  obtained  (see  table  below). 


Portion  of 
Sector 

Area 

(A1) 

Distance  from  Centroid 
to  axis  of  rotation  (c1 ) 

A '  C ' 

a 

(y-88) ( I80-x) 

(I80-x)  +  60  =  360-x 

(y-88) ( 180- x )  (  360-x ) 

2 

3  3 

6 

b 

(y-88) (60) 

60 

(y-88 ) (60 )2 

T 

- 2 - 

c 

( 168-y) ( 240-x ) 

(240-x) 

( 168-y) (240-x  )2 

- T~ 

- ? - 

(y-88) (180-x) (360-x)  +  (y-88)(60)‘  +  (168-y ) (240-x) 


(y-88)  [(180-x) (360-x)  +  10800]  +  (168-y) (240-x) 

3 - 

Ac 

■  6,400,000 


(y-88) [(180-x) (360-x)  +  10,800]  +  (168-y) (240-x) 

Sector  III  (fig.  3A-5) 

EMvn  +  EMvp  =  (20,000) ( 180-X/2)  +  2/3(20,000) (x/2)  + 


(20,000) ( 180-X/2)  +  2/3(20 ,000) (x/2) 
=  360(20,000)  -  (20,000)x/3 


Portion  of 
Sector 


(180-x) (y-88) 


d 

(180-x) (88) 


Distance  from  Centroid 
to  axis  of  rotation  (c1 ) 

A'c1 

(y-88)  +  88  =  y  +  176 

(180-x) (y-88)(y  +  176) 

J.  3 

88 

6 

(180-x)  (88)2 

T 

y 

2 

2 

xyc 

(y-88)  ( 180-x) (y+176 )/6  +  (180-x)  (88)d/2  +  xyV6 
1/6  ( 180-x) [ (y-88) (y+176 )+23 ,232]  +  xy2 


(2160-2x) (20,000) 


(180-x) [(y-88) (y+176  )+23,232]+xy 

Step  6.  Due  to  the  complexity  of  obtaining  a  direct  solution  for  ultimate 
unit  resistance,  a  trial-and-error  solution  will  be  used  (  see  table 


I 


X 

y 

rI 

rII 

rIII 

125 

130 

9.21 

7.67 

9.33 

125 

135 

9.55 

7.92 

8.77 

125 

140 

9.92 

8.19 

8.25 

125 

145 

10.32 

8.48 

7.78 

125 

150 

10.77 

8.79 

7.35 

130 

130 

8.52 

8.29 

9.50 

130 

135 

8.83 

8.55 

8.92 

131 

135 

8.70 

8.68 

3.85 

> 


» 


Therefore: 


x  =  131  ins 
y  =  135  i ns 
ru  =  8.68  psi 


Problem  3A-2  Resistance  -  Defection  Function 

Problem:  Determine  the  acutal  and  equivalent  resistance  deflection  function  in 
the  elasto-plastic  region  for  a  two-way  structural  element. 

Procedure: 

Step  1.  Establish  design  parameters 

a.  Geometry  of  element. 

b.  Support  conditions 

Step  2.  Determine  ultimate  positive  and  negative  moment  capacities. 

Step  3.  Determine  static  properties: 

a.  Moduls  of  elasticity  for  the  element. 

b.  Moment  of  inertia  of  the  element. 

Step  4.  Establish  points  of  interest  and  their  ultimate  moment 
capacities  (fig.  3-23) 

Step  5.  Compute  properties  at  first  yield. 

a.  Location  of  first  yield 

b.  Resistance  at  first  yield  r0 

c.  Moments  at  remaining  points  consistent  with  r? 

d.  Maximum  deflection  at  first  yield. 

Step  6.  Compute  properties  at  second  yield 

a.  Remaining  moment  capacity  at  other  points 

b.  Location  of  second  yield. 

c.  Change  in  unit  resistance  r  between  first  and  second  yield, 

d.  Unit  resistance  at  second  yield  r  . 

e.  Moment  at  remaining  point  consistent  with  rep. 
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f.  Change  in  maximum  deflection. 

g.  Total  maximum  deflection. 

Note: 

An  element  fixed  on  three  or  four  sides  exhibits  three  or  more 
yields  (either  in  the  horizontal  or  vertical  direction). 
Therefore,  repeat  Step  6  as  many  times  as  necessary  to  obtain 
properties  at  the  various  yield  points. 

Step  7.  Compute  properties  at  final  yield  (ultimate  unit  resistance) 

a.  Ultimate  unit  resistance. 

b.  Change  in  resistance  between  ultimate  unit  resistance  and 
resistance  at  prior  yield. 

c.  Change  in  maximum  deflection. 

d.  Total  maximum  deflection. 

Step  8.  Draw  the  actual  resistance-deflection  curve  (fig.  3-39). 

Step  9.  Calculate  equivalent  maximum  elastic  deflection  of  the  element. 


Example:  3A-2  Resistance-Deflection  Function 


Required:  The  actual  and  equivalent  resistance-deflection  function  (curve) 

in  the  elasto-plastic  region  for  the  two-way  structural  steel 
element. 


Sol ution: 

Step  1.  Given: 

a.  L  =  240  in.  H  =  168  in. 

b.  Fixed  on  three  sides  and  free  at  the  fourth. 

Step  2.  Same  as  step  3  of  example  3A-1(A),  part  a. 

Step  3.  Static  properties. 

a.  Modulus  of  elasticity,  E$  for  steel 
Es  =  29  x  106  psi 

b.  Considering  a  1-inch  strip  (b  =  1  inch) 
Assume  I  =  144  i n^ 

Step  4.  For  points  of  interest,  see  figure  3A-6. 


Step  5.  Properties  at  first  yield. 


341 


/"POINT  2  r POINT  I 


POINT  2 


Mhn 

t  mhp 

n 

V 

V 

V 

k 

/•POINT  3  ' 

)mvn  : 

/»  KUIIN 

f  mhn 


Figure  3A-6 

From  figure  3-27  for  H/L  =  0.7 
B1  =  0.077  =  0.160  e3  =  0.115 

Y !  =  0.012  u  =  1/6 

3.  Mj_) p  -  Mpy  =  Myp  =  My jg  =  20,000  in- lbs/ in 
Mp  =  SrH2 
r  =  M/BH2 

r{  =  20 ,000/ [ (0 .077 ) { 168 )2]  =  9.20  psi 
r2  =  20 ,000/C (0 . 160 ) ( 168 )2]  =  4.43  psi 
r3  =  20 ,000/C (0.115)(168)2]  =  6.16  psi 
First  yield  at  point  2  (smallest  r) 

b.  re  =  4.43  psi 

c.  Mp  (Point  1)  =  (0.077  )  (4.43) ( 168)2  =  9,627  in-lbs/in 
My  (Point  3)  =  (0 . 1 15 ) (4 .43) { 168 )2  =  14,379  in-lbs/in 

d.  0  =  EI/b( 1-v2 ) 

=  29  x  106  x  144/1[1-(0.1667)2]  =  43.0  x  108  in-lbs 
xe  =  Yjr^/O  =  (0.0120)  (4.43)  { 168 )4/43(  108)  =  0.0093  in 


p  6.  Properties  at  second  yield. 


After  first  yield  element  assumes  a  simple-simple-fixed-free 
stiffness,  therefore  from  figure  3-29  for  H/L  =0.7. 

B1  =  0.120  63  =  0.220 

Yj  ~  O.045  u  =  0.3 

a.  Mp  (Point  1)  =  MHp  -  Mp  (at  rQ) 

=  20,000  -  9627  =  10373  in-lbs/in 
(Point  3)  =  MVN  -  Mp  (at  rQ) 

=  20,000  -  14,379  =  5621  in-lbs/in 

b.  Mp  (Point  1  )  =  10373  in-lbs/in  =  B^ArH2 

Ar  =  10373/(0. 120X168)-  =  3.U6  psi 
Mn  (Point  3)  =  5,621  in-lbs/in  =  B3ArH2 

Ar  =  5 ,621 /[ (0.220) ( 168 )2]  =.0.90  psi 
Second  yield  at  Point  3  (smaller  Ar) 

c.  Ar  =  0.90  psi 

d.  rep  =  re  +  Ar  =  .4.43  +  0.90  =  5.33  psi 

e.  Mp  (Point  1)  =  0. 120(0. 90 ) ( 163)2 

=  3,048  in-lbs/in 

f.  D  =  E I /b ( 1 - v2 )  =  (29)(106)(144)/1[1-(0.3)2] 

=  45.9  x  108  in-lbs/in 
Ax  =  YlArH4/D  =  0. 030(0. 90 )  ( 1 68 )  4/45 . 9  { 1 0 ) 8 
=  0.0047  in 

g.  Xep  =  Xe  +  AX  =  0.0098  +  0.0047  =  0.014  in 

Step  7.  Properties  at  final  yield  (ultimate  unit  resistance).  After 
second  yield  element  assumes  a  simple-simple-simple-free 
stiffness,  therefore  from  figure  3-30  for  H/L  =  0.7. 


Y1 

=  0.045 

u  •"  0.3 

a. 

ru  =  10.6 

psi  (part  a,  example  3A-1 ( A ) ) 

b. 

Ar  =  ru  - 

rgp  =  10.6  -  5.33  =  5.27  psi 

1  Q 


c.  D  =  EI/b( 1-v2)  =  29(106)(144)/l[l-(0.3)2] 

=  45.9  x  108  in-lbs 

Ax  =  y1  rH4/D  =  (0.045 ) (5 .27 ) ( 168)4/45 ,9  x  108 
=  0.041  in 

d.  Xp  =  Xep  +  AX  =  0.014  +  0.041  =  0.055  in 

Step  8.  For  actual  resistance-deflection  curve,  see  figure  3A-7. 

Step  9.  XE  =  Xe(rep/ru)  +  Xep[l-(re/ru)]  +  Xp[l-( rep/ru) ] 

Equation  3-35 

■  •  --  0.0098(5.33/10.6)  +  0.014  [  1- (4 .43/10 .6 ) ] 

*  C.0S5[l-( 5 .33/10.6) ] 

=  0.00493  +  0.0081  +  0.0273- 
=  0.04033  in 

The  equivalent  resistance-deflection  curve  is  shown  in  figure  3A- 
7. 

Problem  3A-3  Dynamic  Design  Factors  For  A  One  Hay  Element 

Problem:  Determine  the  plastic  load,  mass  and  load-mass  factors  for  a  one¬ 
way  element. 

Procedure: 

Step  1.  Establish  design  parameters. 

Step  2.  Determine  deflected  shape. 

a.  geometry  of  element 

b.  support  conditions 

c.  type  of  load  and  mass 

Step  3.  Determine  maximum  deflection 

Step  4.  Determine  deflection  function 


a.  For  distributed  load  and/or  continuous  mass  determine 
the  deflection  at  any  point. 


UNIT  RESISTANCE  r  (PSI) 


0  .01  .02  .03  .04  .05  .06  .07  .08  .09  .010 

DEFLECTION  X  (INCHES) 


Figure  3A-7 
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b.  For  concentrated  loads  and  concentrated  mass  determine 
the  deflection  at  the  load. 


Step  5.  Calculate  the  shape  function 

a.  For  distributed  load  and/or  continuous  mass  calculate 
4>(x) ,  equation  3-43. 

b.  For  concentrated  load  and  concentrated  mass  calculate 
<f>r,  equation  3-46. 

Steo  6.  Calculate  the  load  factor,  K|_. 

a.  Use  equations  3-41  and  3-42  for  a  distributed  load. 

b.  Use  equations  3-41  and  3-45  for  a  concentrated  load 

Step  7.  Calculate  the  mass  factor,  KM. 

a.  Use  equations  3-47  and  3-48  for  a  continuous  mass 

b.  Use  equations  3-44  and  3-49  for  concentrated  mass 

Step  8.  Calculate  the  load-mass  factor  KLM,  from  equation  3-53. 

Example  3A-3( A)  Dynamic  Design  Factors  For  A  One-Way  Element 

Required:  The  load,  mass  and  load-mass  factors  for  a  structural  steel  beam 

in  the  elastic  range,  with  a  distributed  load. 

Solution: 

Step  1:  Given  structural  steel  beam  shown  in  figure  3A-8 


Figure  3 A -8 


a.  L  =  120  in. 

b.  Simply  supported  on  both  edges 

c.  p ( x )  =  2,000  Ib/in 

m(x)  =  0.0055(  Ib-sVin^l/in 
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Figure  3A-9 

Step  2:  Assumed  deflected  shape  for  elastic  range  is  shown  in 
figure  3A-9 

Step  3:  The  maximum  deflection  at  the  center  is 

<5ma  v  *  5p(  X  )L4 

max  384  El . 

Step  4:  Determine  deflection  function 

6 ( x )  =  p(x)  /■  3  _  o.  y2  .  x3s 

mr u  1 

Step  5:  Calculate  the  shape  function  using  equation  3-43 

*  *  *  P^i  {L3  -  2Lx2  +  x3)  384EI  „ 

6max  24EI  5p( x  )L4 

=  (L3x  -  2Lx3  +  x4) 

5L4 

Step  6:  a.  Using  equation  3-42,  determine  equivalent  force 
Fe  =  /qp( x) 4>( x)dx  =  / q2^(2,000  lb/in)  A£.  (L3x  -  2Lx3  +  x4)dx 


=  6,400  l3x2  -  2Lx4  +  x5  /j20  =  1 .280L/X20 


~T~  ~T~  TT 


=  153,600  lb. 

b.  From  equation  3-41,  find  the  load  factor 
K.  =  FE  =  153,600  lb. 


=  0.64  in  the  elastic  range 


J4 ; 


Step  7:  a.  Find  the  equivalent  mass  from  equation  3-48 

Me  =  7^m(x)4»2(x)dx  =  .0055  256  /J20(L3x  -  2Lx3  +  4)2dx 

li18 

=  1.408  /^2n(L8x2  -  4L*x4  +  2L3x8  +  4L2x8  -  4Lx7  +  x8)dx 
25L8 

=  1.408  L6x3  -  4L4x5  +  2L3x6  +  4L2x7  -  4Lx8  +  x9  ^20 

— 25 — r~  5  5  ~1~  ~T~  Y~ 

=  .00277L/J20 

=  0.3325  lb  -  s2/in3 

b.  From  equation  3-47,  calculate  the  mass  factor 

Km  =  ME  =  0.3325  lb  -  s2/in3 

(0.0055  lb  -  s2/in4  x  120  in) 

Km  =  0.50  in  the  elastic  range 

Step  8:  Calculate  the  load-mass  factor  as  defined  by  equation  3-51 

KLM  =  KM/KL 

=  0.50/0.64 

Klm  *  0.78  in  the  elastic  range 


Example  3A-3(B)  Dynamic  Design  Factors  For  A  One-Way  Element 

Required:  The  load,  mass  and  load-mass  factors  for  a  structural 

steel  beam  in  the  plastic  range  with  a  distributed  load. 

Solution: 

Step  1.  Given  the  structural  steel  beam  shown  in  figure  3A-8 

a.  L  =  120  in. 

b.  Simply-supported  on  both  ends 

c.  p ( x )  =  2,000  lb/in 

d.  m(x)  =  0.0055  (lb  -  S2/in4)/in 

Step  2.  Assume  deflected  shape  for  the  plastic  range  is  shown  in 
figure  3A-10 


Step  3 

Step  4 

Step  5 


Step  6 
a 


b 


Step  7: 
a. 


Figure  3A-10 

Determine  maximum  deflection 
6max  =  (L/2)tane 

Determine  the  deflection  at  any  point. 

5 ( x )  =  xtanG  x  <  L/2 

Calculate  the  shape  function,  equation  3-43 

<p(x)  =  6 ( x )  =  xtanO  • 

TU^rtane 

=  2x/L  x  <  L/2 

Find  Fe  using  equation  3-42. 

Fe  =  /jjp(xH(x)dx  »  2/q°( 2 ,000  lb/in) (2x/L)dx 

=  4,000  lb/in  x2/L  §° 

=  120,000  lb 

From  equation  3-41 

K,  =  El  =  120,000  lb. 

F  '  12,000'  lb/in)"120~'in 

«L  =  0.5  in  the  plastic  range 

Use  equation  3-48  to  find  the  equivalent  mass 

Me  =  /gm(x)<t>2(  x)dx  =  2/ ^  ( 1 . 005  5 }  ( 4  x2/L2)dx 
y3  60 

=  0.044  — “  0 
3l2 

=  0.22  lb  -  s2/in3 
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KM  - 


As  defined  by  equation  3-47 

0.22  1b  -  s2/in3 
(0.0055  1b  -  s2/in4)120  in 
Km  =  0.33  in  the  plastic  range 

Step  8.  Calculate  using  equation  3-53 

KLM  =  KM/KL 
=  0.33/0.5 

Klm  =  0.66  in  the  plastic  range 


Example  3A-3(C)  Dynamic  Design  Factors  For  A  One-Way  Element 


Required:  The  load,  mass  and  load-mass  factors  for  a  structural  steel  beam 

in  the  elastic  range,  with  a  concentrated  load. 


Solution: 


Step  1:  Given  structural  steel  beam  shown  in  figure  3A-11 


l/2 


Figure  3A-11 


a.  L  =  120  in. 

b.  Simply  supported  on  both  sides 

c.  F  =  240  kips 

m(x)  =  0 . 0055 ( 1  b  -  s2/in4)/i n. 

Step  2:  Assume  deflected  shape  for  elastic  range  is  shown  in  figure 
3A-12 


b.  Using  equation  3-41,  calculate  the  load  factor 


K|_  =  F£_  =  240  kips 


=  1.0  for  the  elastic  range 

Step  7:  a.  Equation  3-48  gives  the  equivalent  mass. 

Me  =  /^m(x)4>2(x)dx  =  /A20(0.0055)  (9!_4x2  .  24L2x4  +  16x6)dx 

L6 

=  0.0055  3L4x3  -  24L2x5  +  16x7  i20 

~  ~ 5 -  T” 

=  0.0027L/q20 
=  0.321b  -  s2/i n3 

b.  From  equation  3-47,  calculate  the  mass  factor 


Km  *  JX 


0.321b  -  s2/i n3 


(0.00551b  -  sVin4  x  120in) 


Km  =  0.49  in  the  elastic  range 

Step  8:  Calculate  the  load-mass  factor,  from  equation  3-53 

klm  =  <Vkl 

=  0.49/1.0 


Ki  M  =  0.49  for  the  elastic  range 


Example  3A-3(D)  Dynamic  Design  Factors  For  A  One-Way  Element 

Required:  Determine  the  load,  mass  and  the  load-mass  factors  for  a 

structural  steel  beam,  in  the  plastic  range,  with  a  concentrated 
1  oad. 

Solution: 

Step  1.  Given  structural  steel  beams  shown  in  figure  3A-11. 

a.  L  =  120  in 

b.  Simply-supported  at  both  edges 
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c.  F  =  240  kips 

m(x)  =  0.0055  (lb  -  s^/in4)/in 

Step  2.  Assumed  deflected  shape  for  the  plastic  range  is  shown  in 
figure  3A-13 


L/2 

L/2 

*  ■ 

*  , 

\ 

' 

8r=Smax^ 

Figure  3A-13 


Step  3:  Determine  maximum  deflection 
5 max  =  (L/2)tane 

Step  4:  Determine  deflection  function. 

a.  for  continuous  mass 

6(x)  =  xtano 

b.  for  a  concentrated  load 

<5  r  =  (L/2)tanG 

Step  5:  Calculate  shape  factors  using 

a.  equation  3-43  for  continuous  mass 

$(x)  =  6 (x)  =  xtanO 

6 max  (L/S»tan* 

=  2x/L 

b,  equation  3-46  for  concentrated  load 

=  <5„  =  (1/2  )tanG 

<ax  (L/Zlt-ne 

=  1.0 


x  <  L/2 


x  <  L/2 
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Step  6:  a.  The  equivalent  force  is  found  using  equation  3-45 


Step  7 : 


FE  -  E  Fr*r  =  Pxl 
r=l 

=  240  kips 

b.  Equation  3-41  gives  the  load  factor 
Kl  =  FE  =  240  kips 


Step  8: 


«L  «  1.0  for  plastic  range 

a.  The  equivalent  mass  is  found  using  equation  3-48 
Me  =  / Qm( x )4> ^ ( x ) dx  3  2/q®(0.0055) (4x2/l2)dx 

=  0.44  x3  &° 

—j 

,  =  0.22  lbL-  s2/in3  - 

b.  Solve  for  KM  using  equation  3-47 

Km  =  mE  =  0.221b  -  s2/in3 

^  (0.00551b  -s2/ in4)  120  in 

Km  *  0.33  in  the  plastic  range 
From  equation  3-53,  calculate  klm 

klm  =  km/kl 

=  0.33  in  the  plastic  range 

Problem  3A-4  Plastic  Load-Mass  Factor 


Problem:  Determine  the  platic  load-mass  factor  klm  for  a  two-way  element 

using  (1)  general  solution  and  (2)  chart  solution. 

Note:  The  determination  of  the  plastic  load-mass  factor  follows 

the  calculations  for  the  ultimate  resistance,  hence  the  structural 
configuration  and  the  location  of  the  plastic  yield  lines  will  be 
known. 

Procedure:  Part  (a)  -  General  Solution 

Step  1.  See  part  a,  problem  3A-1  for  the  structural  configuration 
and  location  of  plastic  yield  lines.  Denote  sectors 
formed  by  yiel d  1  ines. 


Procedure: 


Requi red: 
Solution: 


Step  2.  Determine  the  load-mass  factors  properties  I,  c,  and  L' 
for  all  sectors. 

Step  3.  Determine  the  factor  I/cL'  for  all  sectors. 

Step  4.  Calculate  the  total  area  of  the  element. 

Step  5.  With  values  obtained  above,  calculate  the  plastic  load- 
mass  factor  for  the  element  using  equation  3-57. 

Note:  In  the  above  problem,  an  element  of  uniform 
thickness  was  considered.  For  non-uniform  elements,  the 
load-mass  factor  is  calculated  using  equation  3-53  where 
the  mass  of  the  individual  sectors  must  be  considered. 

Part  (b)  -  Chart  Solution 

Step  I.  See  part  b,  problem  3A-1  for  structural  configuration  and 
location  of  plastic  yield  lines  in  terms  of  x/L  or  y/H. 

Step  2.  For  known  value  of  X/L  or  y/H  and  support  condition, 

determine  the  load-mass  factor  for  the  element  from  figure 
3-44. 

Note:  Chart  solution  may  be  used  only  if  the  element 
conforms  to  the  requirements  listed  in  section  3-17.3 


Example  3A-4  Plastic  Load-Mass  Factor 

Plastic  load-mass  factor  for  the  element  considered  in  example  3A- 
1(A)  using  (1)  general  solution  and  (2)  chart  solution. 

Part  (a)  -  General  Solution 

Step  1.  Given  structural  configuration  and  location  of  yield  lines 
shown 

below  (see  part  a,  example  3A- 1(A))  in  figure  3A-14. 


L  =  240  in 
H  =  168  in 
X  =  120  in 
y  =  137.6  in 
Tc  =  constant 


Figure  3A  -  14 
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Step  2.  Load-mass  factor  properties. 

a.  Sector  I . 

L1  =  y  =  137.6  in 

c  =  y/3  =  137.6/3 

I  =  L(L ' )3/12  =  240(137. 6)3/12 

b.  Sector  II 


-AXIS  OF 
ROTATION 


Figure  3A  -  16 


L 1  =  x  =  120  in 


l  • 


Figure  3A-15 


H  -  y  *  168  -  137.6  =  30.4  in 

L'  [H  *  2  (H-y)]  =  120  [168  +  2  (30.4)] 


C  =  120  (0.384) 


(H-y) (L1 )3  +  y(L')3 


30.4  { 1 20 ) 3  +  137.6(1 20 ) 3  *  21.60(120)3 


Step  3.  Calculate  factor  I/cL'  for  euch  sector: 
240(137. 6)3 


Sector  I. 


Sector  II. 


cL'  ,  133.4 

v  — - —  ) ( 133.4) 

I  21 . 60 ( 120 ) 3 


Sector  III.  =  6,646  in' 


Step  4.  Area  of  panel 
A  =  LH  =  240  (168)  =  40,320  in2 


8,256  in2 


6,646  in2 


%  .% 
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Step  5. 


Load-mass  factor 


klm 

I/cL' 

‘  - 1 - 

(eq.  6-14) 

klm 

_  8,256  +  2(6,646)  _  n  c,„ 

'  - 5U7370 -  '  °*b34 

Solution:  Part  (b) 

-  Chart  Solution 

Step  1.  Given:  Panel  fixed  on  3  edges,  1  free  and  y/H  =  0.803  (see  part 
b,  example  3A-1(A) ) . 

Step  2.  From  figure  3-44,  read  load-mass  factor 

Klm  =  0.543 


Problem  3A-5  Response  of  a  Si ngl e-Degree-of  Freedom  System  subject  to  Dynamic 
Load 

Problem:  Determine  the  maximum  response  and  the  corresponding  time  it  occurs 
of  a  si ngl e-degree-of- freedom  system  subjected  to  dynamic  load  using 
(a)  numerical  methods  and  (b)  design  charts. 

Procedure:  Part  (a)  -  Numer1'.  '•  Methods 

Step  1.  Establish  dimensional  parameters  of  the  system. 

Step  2.  Determine  the  natural  period  of  vibration  and  integration  time 
interval . 

Step  3.  Construct  a  table  similar  to  table  3-14  of  section  3-19.2. 

Note:  For  the  first  interval  n=l,  Equation  3-59  is  used  and 
subsequent  intervals,  the  recurrence  formula  (eqn.  3-56)  is  used. 

Procedure:  Part  (b)  -  Chart  solution 

Step  1.  Same  as  step  1  of  example  3A-5,  part  a. 

Step  2.  Determine  the  non-dimensional  parameters. 

Step  3.  Determine  the  ratio  of  the  maximum  desplacement  to  the  elastic 

displacement  X^/X^-  and  the  ratio  of  the  time  at  which  this  maximum 
displacement,  occurs  to  the  duration  of  the  blast  load. 


Example  3A-5  Maximum  Response  of  Single-Degree-of-Freedom  System  Subjected  to 
a  Triangular  Load. 

Required:  The  maximum  response  and  the  time  it  occurs,  of  a  si  ngl e-degree-of- 
freedom  system  subjected  to  blast  loads,  using  (a)  numerical 
methods  and  (b)  design  charts. 
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Solution:  Part  (a)  -  Numerical  Methods 


Step  1.  Given: 


p 


time  displacement 

(a)  Si ngle-Degree-of-  (b)  External  Load  (c)  Resistance 

Freedom  System  Function 

Figure  3A-17  Slngle-Degree-of-Freedom  System  subjected  to  external  loads. 

m  =  2.5'  Kips-sec2/ft 
K  =  9,860  Kips/ft 
ru  =  750  Kips 
XE  =  0.076  ft 
T  =  0.10  sec 
J>  =  1000  Kips 

Step  2.  Natural  period  of  vibration  and  integration  time  interval. 
Tn  =  2  [m/K]1/2  =  2(  )  [2 .5/9 ,860]1/2  =  0.10  sec 

t  -  Tfj/10  =  0.01  sec 

Step  3.  Construct  table  as  shown  below. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

n 

t 

pn 

«n 

pn-«n 

V(pn*Rn)/m 

2X 

n 

■ 

(sec) 

(Kips) 

(Kips) 

(Kips) 

( ft/sec2) 

(ft) 

0 

n 

1000 

0 

1000 

400 

0.040 

0.0 

0.0 

0.020 

1 

900 

197.200 

702.800 

281.120 

0.028112 

0.0400 

0.0 

0.068112 

2 

800 

671.584 

128.426 

51.366 

0.05137 

0.13622 

0.020 

0.121357 

3 

700 

750 

-50.0 

-20.0 

-0.00200 

0.242714 

0.06811 

0.17261 

4 

0.04 

600 

750 

-150.0 

-60.0 

-0.00600 

0.34522 

0.121357 

0.21786 

5 

0.05 

500 

750 

-250.0 

-100.0 

-0.0100 

0.43573 

0.17261 

0.25312 

0.06 

40 

750 

-350.0 

-140.0 

-0.0140 

0.050623 

0.21786 

0.27437 

0.07 

300 

750 

-450.0 

-180.0 

-0.0180 

0.54874 

0.25312 

0.27762 

8 

0.08 

200 

750 

-550.0 

-220.0 

-0.0220 

0.55525 

0.27437 

0.25880 

9 

0.09 

,  =  f(tn)  =  1000  [1-nUt/T)]  KXn  for  Xn  <  X0 

ru  =  750  for  Xn  >  X0 

xn+l  =  2Xn  -  xn-l  +  *n  ^)2 
Note: 

For  n=0,  Xn+1  (Column  10)  =  XQ+1  *  X1  -  ( 1/2 ) aQ(^t)2 

=  (1/2) (0.040) 

=  0.02  ft 

For  n=l ,  2Xn  (Column  8)  =  2Xj_  =  2[(l/2)a0(^t)2] 

=  2(0.02)  =  0.04  ft 
Xn_j  (Column  9)  =  XQ  =  0.0 

Xn+1  (Column  10)  =  X2  =  2Xj-Xq  +  a^At)2 

=  2(0.02)  -  0  +  0.02811 
=  0.06811  ft 

For  n=2 ,  2Xn  (Column  8)  =  2X2  =  2(0.06811) 


0.13622  ft. 


Xn-1  (Column  9)  *  Xj_  »  0.02  ft. 

Xn+1  (Column  10)  =  X3  =  2X2  -  X1  +  a2  (At)2 

=  (2) (0.06811 )  -  0.02  +  0.005137 
=  0.12:357  ft. 

For  n=3,  4,...,  repeat  the  above  procedure. 

Solution:  Part  (b)  -  Design  Charts 


Step  1.  Same  as  step  1  of  example  3A-5,  part  a 
Step  2.  Non-dimensional  parameters 

a.  Natural  period  of  vibration,  Tn 

Tn  =  2ir[m/K]12  =  2it[2.5/9,860]12  =  0.10  sec 

b.  Ratio  of  duration  of  blast  load  T  to  natural  period 

T/Tn  =  0.10/0.10  =  1.0 

c.  Ratio  of  peak  resistance  ru  to  peak  load  P 

ru/P  *  750/1000  =  0.75 

Step  3.  Using  the  ratios  calculated  in  step  2  and  figures  3-54  and  3-55, 
determine  the  value  of  Xm/XE  and  t^/TN. 

For  T/Tn  =  1  and  ru/P  =  0.75 

Xm/XE  =  3.7  from  figure  3-54 

VT  *  0  .77  from  figure  3-55 

Step  4.  Determine  Xm  and  t^ 

Xm/XE  3.7 

Xm  =  (3.7)XE  =  (3.7(ru/KE) 

=  ( 3 . 7 ) ( 750/9,860 )  =  0.28144  ft. 

VT  =  °-77 

t-  =  (0. 77 ) T  =  0.77(0.10)  =  0.077  sec 


Problem  3A-6  Maximum  Response  of  a  Single-Degree-of-Freedom  System  to 
Bilinear  Blast  Loads 

Problem:  Determine  Xm/Xr,  t^/Tjg  and  t^/T  (when  applicable)  for  a  single- 
degree-of-freeaom  system  subject  to  various  bilinear  blast  loads. 

Procedure:  Part  (a)  -  Solution  in  Region  0 

Step  1.  Establish  normalized  parameters 

Step  2.  Enter  table  3-15  with  the  given  C  parameters  and  determine 
which  figures  have  to  be  used. 

Step  3.  Enter  each  of  the  figures  determined  in  step  2,  with  the  given 

values  of  the  other  two  parameters  and  determine  the  region  where 
the  intersection  points  are  located. 

Step  4.  Based  on  the  region  where  the  intersection  points  are  located, 
enter  the  appropriate  figure  and  find  Xm/XE,  t^/T^  and  t^/T. 


Procedure:  Part  (b)  -  Solution  in  Region.  C  -  Graphical  Interpolation. 


Step  1.  Same  as  step  1  in  part  a. 

Step  2.  Same  as  step  2  in  part  a. 

Step  3.  Same  as  step  3  in  part  a. 

Step  4.  Set  up  a  table  as  shown  in  table  3A-1.  Post  each  figure  number 

and  the  corresponding  values  of  and  C^,  leaving  a  space 
between  each  line  of  information.  Post  in  the  spaces  the 
appropriate  valuesof  Ci  and  Co  needed  for  interpolation.  Enter 
each  of  the  figures  determined  in  Step  2  with  the  given 
parameters  and  find  the  values  of  X  /Xr.  Post  these  values  in 
table  3A-1 . 

Step  5.  Use  log-log  graph  paper  to  plot  the  points  obtained  in  Step  4. 

Post  these  values  in  table  3A-1,  using  linear  interpolation  where 
necessary. 

Step  6.  Plot  on  log-log  graph  paper  the  points  which  represent  (Xm/Xr, 

C2 )  for  the  given  value  of  C.  Use  linear  interpolation  to  find 
Xm/XE  for  given  value  of  C2. 

Procedure:  Part  (c)  -  Solution  in  Region  C  -  Mathematical  Interpol ation 

Step  1.  Same  as  step  1  in  part  a. 

Step  2.  Same  as  step  2  in  part  a. 

Step  3.  Same  as  step  3  in  part  a. 
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Step  4.  Same  as  step  4  in  part  b. 

Step  5.  Solve  for  lnYa  and  lnY^  using  equations  3-83  and  3-84. 

Step  6.  Solve  for  InY  using  equation  3-85. 

Step  7.  Solve  for  Y  using  equation  3-86. 

Example  3A-6  Maximum  Response  of  a  Single-Degree-of-Freedora  System  to 
Bilinear  Blast  Loads 

Required:  Determine  X^/Xp,  and  t^/T  (when  applicable)  for  a  single- 

degree-of- freedom  system  subject  to  various  bilinear  blast  loads. 

Solution:  Part  (a)  -  Solution  in  Region  D 


Step  1.  Given: 


P/ru  =  1.0 
T/Tw  =  3.0 


=  0.66 


=  50 


Step  2.  Enter  table  3-15  with  Ci  =  0.66  and  Co  =  50.  Note  fiqures 
3-119,  3-120,  3-147,  and  3-148  apply. 

Step  3.  Enter  each  of  the  figures  determined  in  step  2,  with  P/ru  =  1.0 
and  T/T»  =  3.0.  Note  that  the  intersection  point  is  located  to 
the  right  of  the  line  of  solid  squares,  defined  as  region  0.  In 
region  D,  the  maximum  dynamic  response  depends  only  on  the  shock 
load  described  by  P/ru  and  T/T the  gas  load  described  by  C j_P/rL 
and  C2T/Tn  does  not  influence  the  maximum  dynamic  response. 
Consequently,  figures  3-54a  and  3-64b  for  a  single  triangular 
lr>-jr.  pulse  apply.  Enter  figure  3-64a  with  P/r  =  1.0  and  T/T».  = 
3.0  and  find  Xm/X£  =  3.55,  tm/TN  =  0.98,  tE/T  =  0.086. 

Solution:  Part  (b)  -  Solution  in  Region  C  -  Graphical  Interpolation 

Step  1.  Given:  P/ru  =  32 
T/Tm  =  0.10 


C2  =  20 

Step  2.  Enter  table  3-15  with  Ci  =  0.06  and  Co  =  20.  Note  figures  3-108, 
3-109,  3-133  and  3-134  apply. 

Step  3.  Enter  each  of  these  figures  with  P/ru  =  32  and  T/Tn  =  0.10.  Note 
that  the  intersection  point  is  not  located  in  regions  A,  B  or 
D.  Therefore  the  intersection  points  lie  in  region  C  and 
interpolation  between  charts  is  required  to  obtain  a  solution. 


Step  4, 


Set  up  table  as  shown  in  table  3A-1  below.  Post 
number  and  the  corresponding  values  of  and  C< 


between  each  line  of 
appropriate  values  of  C1 
figure  3-108  with  P/ru  = 
112.  Post  this  value  in 
=  32  and  T /T^  =  0.10  and 
table, 
values 


each  chart 
leaving  a  space 


information.  Post  in  the  spaces  the 
and  Co  needed  for  interpolation. 
32  and  T/TN 
the  table. 


Enter 

=  0.10  and  find  X^/X^  = 

Enter  figure  3-109  with  P/ru 
86.  Post  this  value  in  the 
Repeat  this  process  for  'figures  3-133  and  3-134,  and  post 


fi  nd 


VXF  = 


for 


Vxe 


in  the  table. 


Step  5.  Use  log-log  graph  paper  to  plot  the  points  (112,0.068)  and 


(86,0.046)  which  represent 


<V*E: 


C^)  for  C2  =  10  as  shown 


figure  3-267.  Use  straight-li'ne“interpolation  to  find 
103  for  Ci  =  0.060.  Post  this  value  in  the  table. 


C1  = 

process  for  Co  =  30, 
in  figure  3-257  . 


1  n 

tr/xE  = 
Repeat  this 


and  find 


V*E 


=  251  for  C^  =  0.06  as  shown 


Step  6.  Plot  on  log-log  graph  paper  the  points  (103,10)  and  (251,30) 
which  represent  { Xm/Xg- , C2 )  for  Cj  =  0.060.  Use  straight-line 


interpolation  for 


VXE 

3-267.  Thus  the  solution  is- X 


for  C2  =  50  as  shown  in  figure 


m 


/XE  =  182. 


Table  3A-1 


Figure  No. 

Cl 

c2 

VXE 

3-108 

0.068 

10 

112 

0.060 

10 

— 

> 

103 

3-109 

0.046 

10 

86 

0.060 

20 

-- 

>  182 

3-133 

0.075 

30 

340] 

0.060 

30 

> 

251 

3-134 

0.056 

30 

230 

Solution: 

Part 

(c) 

-  Solution  in 

region 

C  -  Mathematical  Interpol ation 

Step 

1. 

Same 

as 

step 

1 

of 

part 

(b). 

Step 

2. 

Same 

as 

step 

2 

of 

part 

(b). 

ctep 

3. 

Same 

as 

step 

3 

of 

part 

(b). 

Step 

4. 

Same 

as 

step 

4 

of 

part 

(b). 

Step 

5. 

Us  i  ng 

equation 

3-83  and 

3-84, 

find  1 nY a  and  lnY^. 

Step  6. 


Step  7. 


V 


InY,  =  InYi  +  ln(  Y^/Y,  H  ntC,  /C, , ) 

1  l  ^ 

=  1nll2  +  ln(86/112)ln(0. 060/0. 068) 

- [HIOT7OT - 

1nYa  =  4.6339 

1nYb  =  1nY3  +  I  n(Y4/Y3)  1  n(Cl  /Cul 
ln(C14/C13) 

=  1 n340  +  1 n( 230/340  )1 n{0. 06 /o.075) 

- ln(  O' .'056/07075} - 

1nYa  =  5.5304 

Find  InY  from  equation  3-85 

InY  =  lnYa  +  lnYb  "  1  "V1  n(C2/C21 } 
-rTnTC237^)"- - 

=  4.6339  +  (5.5304  -  4.639)1 n(20/10) 

- TnTOTnN - 

InY  =  5.1995 

Solve  for  Y  using  equation  3-86 

Y  =  e1nY 

=  .5.1995 

Y  =  131 
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APPENDIX  38  -  LIST  OF  SYMBOLS 


a  (1)  acceleration  (in. /ms2) 

(2)  depth  of  equivalent  rectangular  stress  block  (in.) 

A  area  (in.2) 

Aa  area  of  diagonal  bars  at  the  support  within  a  width  b  (in.  ) 

Aq  area  of  openings  (ft2) 

A$  area  of  tension  reinforcement  within  a  width  b  (in.  ) 

Aj  area  of  compression  reinforcement  within  a  width  b  (in.  ) 

A$h  area  of  flexural  reinforcement  within  a  width  b  in  the  horizontal 

direction  on  each  face  (in.2)* 

A$v  area  of  flexural  reinforcement  within  a  width  b  in  the  vertical 

direction  on  each  face  (in.2)* 

Ay  total  area  of  stirrups  or  lacing  rei nforcement  in  tension  within  a 

distance,  ss  or  s^  and  a  width  b$  or  b1  (in.*-). 

Aj,  Aj j  area  of  sector  I  and  II,  respectively  ( 1  n.2) 
b  (1)  width  of  compression  face  of  flexural  member  (in.) 

(2)  width  of  concrete  strip  in  which  the  direct  shear  stresses  at 
the  supports  are  resisted  by  diagonal  bars  (in.) 
b$  width  of  concrete  strip  in  which  the  diagonal  tension  stresses  are 

resisted  by  stirrups  of  area  Ay  (in.) 

b-j  width  of  concrete  strip  in  which  the  diagonal  tension  stresses  are 

resisted  by  lacing  of  area  Ay  (in.) 

B  constant  defined  in  pargraph 
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c 


(1)  distance  from  the  resultant  applied  load  to  tho  axis  of 
rotaton  (in.) 

(2)  damping  coefficient 

Cj.Cjj  distance  from  the  resultant  applied  load  to  the  axis  of  rotation 
for  sectors  I  and  II,  respectively  (in.) 
cs  dilatational  velocity  of  concrete  (ft/sec) 

C  shear  coefficient 

Ccr  critical  damping 

shear  coefficient  for  ultimate  shear  stress  of  one-way  elements 

?  d  ft 

Cf  post-failure  fragment  coefficient  (lb  -ms  /in.  ) 

Cra  peak  reflected  pressure  coefficient  at  angle  of  incidence  a 

C$  shear  coefficient  for  ultimate  support  shear  for  one-way  elements 

C$H  shear  coefficient  for  ultimate  support  shear  in  horizontal 

direction  for  two-way  elements* 

CsV  shear  coefficient  for  ultimate  support  shear  in  vertical  direction 

for  two-way  elements* 

Cq  drag  coefficient 

CQq  drag  pressure  (psi) 

CDq0  peak  drag  pressure  (psi) 

C£  equivalent  load  factor 

CH  shear  coefficient  for  ultimate  shear  stress  in  horizontal 

direction  for  two-way  elements* 

CL  leakage  pressure  coefficient 


See  note  at  end  of  symbols. 


3btf 


maximum  shear  coefficient 


CM 

Cu  impulse  coefficient  at  deflection  Xu  ( psi-ms2/in.2) 

Cy  impulse  coefficient  at  deflection  Xm  (psi-ms2/in.2) 

Cy  shear  coefficient  for  ultimate  shear  stress  in  vertical  direction 

for  two-way  elements* 

Cj  (1)  impulse  coefficient  at  deflection  Xj  (psi-ms2/in.2) 

(2)  parameter  defined  in  figure 

(3)  ratio  of  gas  load  to  shock  load 

Cj  impulse  coefficient  at  deflection  Xm  ( psi-ms2/in.2) 

C£  ratio  of  gas  load  duration  to  shock  load  duration 

d  distance  from  extreme  compression  fiber  to  centroid  of  tension 

reinforcement  (in.) 

d'  distance  from  extreme  compression  fiber  to  centroid  of  compression 

reinforcement  (in.) 

dc  distance  between  the  centroids  of  the  compression  and  tension 

reinforcement  (In.) 

de  distance  from  support  and  equal  to  distance  d  or  dc  (in.) 

d^  Inside  diameter  of  cylindrical  explosive  container  (in.) 

d-|  distance  between  center  lines  of  adjacent  lacing  bends  measured 

normal  to  flexural  reinforcement  (in.) 

dco  diameter  of  steel  core  (in.) 

dj  diameter  of  cylindrical  portion  of  primary  fragment  (in.) 


See  note  at  end  of  symbo 1 s . 
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D 


°o 

°E 

OIF 

DLF 

e 

(2E')1/2 

E 

Ec 

Es 

f 


fc 

fdc 

fds 

fdu 

fdy 

fs 


fu 

fy 


(1)  unit  flexural  rigidity  (lb- in.) 

(2)  location  of  shock  front  for  maximum  stress  (ft) 

(3)  minimum  magazine  separation  distance  (ft) 
nominal  diameter  of  reinforcing  bar  (in.) 

equivalent  loaded  width  of  structure  for  non-planar  wave  front 
(ft) 

dynamic  increase  factor 
dynamic  load  factor 

base  of  natural  logarithms  and  equal  to  2.71828... 

Gurney  Energy  Constant  (ft/sec) 

modulus  of  elasticity 

modulus  of  elasticity  of  concrete  (psi) 

modulus  of  elasticity  of  reinforcement  (psi) 

unit  external  force  (psi) 

static  ultimate  compressive  strength  of  concrete  at  28  days  (psi) 

dynamic  ultimate  compresive  strength  of  concrete  (psi) 

dynamic  design  stress  for  reinforcement  (psi) 

dynamic  ultimate  stress  of  rei nforcement  (psi) 

dynamic  yield  stress  of  reinforcement  (psi) 

static  design  stress  for  reinforcement  (a  function  of  fy,  fu  and  9 
(psi) 

static  ultimate  stress  of  reinforcement  (psi) 
static  yield  stress  of  rei nforcement  (psi) 

(1)  total  external  force  (lbs) 

(2)  coefficient  for  moment  of  inertia  of  cracked  section 
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(3)  function  of  C2  S  Cj  for  bilinear  triangular  load 

force  in  the  reinforcing  bars  (lbs) 

equivalent  external  force  (lbs) 

variable  defined  in  table  4-3 

charge  location  parameter  (ft) 

(1)  span  height  (in.) 

(2)  distance  between  reflecting  surface(s)  and/or  free  edge(s)  in 
vertical  direction  (ft) 

height  of  charge  above  ground  (ft) 

scaled  height  of  charge  above  ground  (ft/lb1/3) 

height  of  structure  (ft) 

scaled  height  of  triple  point  (ft/lb*/3) 

unit  positive  impulse  (psi-ms) 

unit  negative  impulse  (psi-ms) 

sum  of  scaled  unit  blast  impulse  capacity  of  receiver  panel  and 
scaled  unit  blast  impulse  attenuated  through  cocnrete  and  sand 
in  a  composite  element  {psi-ms/lb3/3) 
unit  blast  impulse  (psi-ms) 
scaled  unit  blast  impulse  ( psi -ms/1 b1 ) 

total  scaled  unit  blast  impulse  capacity  of  composite  element 
(psi-ms/lb1/3) 

scaled  unit  blast  impulse  capacity  of  receiver  panel  of  composite 
element  (psi -ms/1 bx ^3 ) 

scaled  unit  blast  impulse  capacity  of  donor  panel  of  composite 
element  (psi-ms/lb^3) 


V 


1-e 

V 

V 
U 

V 

I 

!a 

*g 

j 

k 

K 

Ke 

KeP 

KE 

kl 

klm 

^LNpu 

(KLM^ud 

kM 


unit  excess  blast  impulse  (psi-ms) 

unit  positive  normal  reflected  impulse  (psi-ms) 

9 

unit  negative  normal  reflected  impulse  (psi-ms) 
unit  positive  incident  impulse  (psi-ms) 
unit  negative  incident  impulse  (psi-ms) 
moment  of  inertia  (in.4) 

average  of  gross  and  cracked  moments  of  inertia  of  width  b  (in.4) 
moment  of  inertia  of  cracked  concrete  section  of  width  b  (in.4) 
moment  of  inertia  of  gross  concrete  section  of  width  b  (in.4) 
mass  moment  of  inertia  (lb-ms^-in.) 

ratio  of  distance  between  centroids  of  compression  and  tension 
forces  to  the  depth  d 
constant  defined  in  paragraph 

(1)  unit  stiffness  (psi-in  for  slabs)  (Ib/in/in  for  beams) 

(2)  constant  defined  in  paragraph 

elastic  unit  stiffness  (psi/in  for  slabs)  (Ib/in/in  for  beams) 
elasto-pl astic  unit  stiffness  (psi-in  for  slabs)  (psi  for  beams) 
equivalent  elastic  unit  stiffness  (psi-in  for  slabs)  (psi  for 
beams) 

equivalent  spring  constant 
load  factor 
load-mass  factor 

load-mass  factor  in  the  ultimate  range 
load-mass  factor  in  the  post-ultimate  range 
mass  factor 
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» 


*-wb*  ^wd 


resistance  factor 

factor  defined  in  paragraph 

kinetic  energy 

charge  location  parameter  (ft) 
spacing  of  same  type  of  lacing  bar  (in.) 

(1)  span  length  (In.)  except  in  chapter  4  (ft)* 

(2)  distance  between  reflecting  surface(s)  and/or  free  edge(s)  in 
horizontal  direction  (ft) 

length  of  lacing  bar  required  in  distance  s-|  (in.) 
embedment  length  of  reinforcing  bars  (in.) 
wave  length  of  positive  pressure  phase  (ft) 
wave  length  of  negative  pressure  phase  (ft) 
wave  length  of  positive  pressure  phase  at  points  b  and  d, 
respectively  (ft) 

total  length  of  sector  of  element  normal  to  axis  of  rotation  (in.) 
unit  mass  (psi-ms^/in.) 

average  of  the  effective  elastic  and  plastic  unit  masses 
(ps1-ms2/in.) 

effective  unit  mass  (psi-ms2/in. ) 

effective  unit  mass  in  the  ultimate  range  (psi-ms  /in.) 

p 

effective  unit  mass  in  the  post-ultimate  range  (psi-msVin.) 

(1)  unit  bending  moment  (in-lbs/in.) 

(2)  total  mass  (lb-ms2/in.) 


lee  note  at  eba  of  symbols 


Me  effective  total  mass  (lb-ms^/in.) 

Mu  ultimate  unit  resisting  moment  (in-lbs/in.) 

Mc  moment  of  concentrated  loads  about  line  of  rotation  of  sector 

( in. -lbs) 

Ma  fragment  distribution  parameter 

Me  equivalent  total  mass  (lb-ms^/in.) 

Mhn  ultimate  unit  negative  moment  capacity  in  horizontal  direction 

( in. -lbs/in. )* 

MHp  ultimate  unit  positive  moment  capacity  in  horizontal  direction 

( in. -lbs/in. )* 

Mn  ultimate  unit  negative  moment  capacity  at  supports  (in. -lbs/in.) 

Mp  ultimate  unit  positive  moment  capacity  at  midspan  (in. -lbs/in. ) 

Mym  ultimate  unit  negative  moment  capacity  in  vertical  direction  (in.- 

lbs/in.)* 

My p  ultimate  unit  positive  moment  capacity  in  vertical  direction 

(in. -lbs/in. )* 

n  (1)  modular  ratio 

(2)  number  of  time  intervals 
N  number  of  adjacent  reflecting  surfaces 

Nf  number  of  primary  fragments  larger  than  Wf 

p  reinforcement  ratio  equal  to  As  or  L 

TicT  Fn_ 


See  note  at  end  of  symbo 1 s . 
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I 


psb.  pse 


reinforcement  ratio  equal  to  or  Aj 

bd  Wc 

t 

reinforcement  ratio  producing  balanced  conditions  at  ultimate 
strength 

mean  pressure  in  a  partially  vented  chamber  (psi) 
peak  mean  pressure  in  a  partially  vented  chamber  (psi) 
reinforcement  ratio  in  horizontal  direction  on  each  face* 
reinforcement  ratio  equal  to  +  pv 
reinforcement  ratio  in  vertical  direction  on  each  face* 
distributed  load  per  unit  length 

(1)  pressure  (psi) 

(2)  concentrated  load  (lbs) 
negative  pressure  (psi) 

interior  pressure  within  structure  (psi) 
interior  pressure  increment  (psi) 
fictitious  peak  pressure  (psi) 
peak  pressure  ( psi ) 

peak  positive  normal  reflected  pressure  (psi) 
peak  negative  normal  reflected  pressure  (psi) 
peak  reflected  pressure  at  angle  of  incidence  a  (psi) 
positive  incident  pressure  (psi) 

positive  incident  pressure  at  points  b  and  e,  respectively  (psi) 
peak  positive  incident  pressure  (psi) 
peak  negative  incident  pressure 


.ee  note  at  end  of  symbols. 


^sob*  ^sod' 

P 

rsoe 

q 


%> 

% 

%b’^oe 

r 


peak  positive  incident  pressure  at  points  b,  d,  and  e, 
respectively  (psi) 
dynamic  pressure  (psi) 

dynamic  pressure  at  points  b  and  e,  respectively  (psi) 
peak  dynamic  pressure  (psi) 

peak  dynamic  presure  at  points  b  and  e,  respectively  (psi) 

(1)  unit  resistance  (psi) 

(2)  radius  of  spherical  TNT  (density  equals  95  lb/ft^  charge 


(ft)) 


r-  unit  rebound  resistance  (psi) 

r  change  in  unit  resistance  (psi) 

re  elastic  unit  resistance 

rep  elasto-plastic  unit  resistance  (psi) 

ru  ultimate  unit  resistance  (psi,  for  slabs)  (lb/in  for  beams) 

rup  post-ultimate  unit  resistant  (psi) 

rj  radius  of  hemispherical  portion  of  primary  fragment  (in.) 

R  (1)  total  internal  resistance  (lbs) 

(2)  slant  distance  ( ft) 

Rf  distance  traveled  by  primary  fragment  (ft) 

R1  radius  of  lacing  bend  (in.) 

RA  normal  distance  (ft) 

R£  equivalent  total  internal  resistance  (lbs) 

Rg  ground  distance  (ft) 

Ru  total  ultimate  resistance 

Rj ,Rjj  total  internal  resistance  of  sectors  I  and  II,  respectively  (lbs) 
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spacing  of  stirrups  in  the  direction  parallel  to  the  longitudinal 
reinforcement  (in.) 

spacing  of  lacing  in  the  direction  parallel  to  the  longitudinal 
reinforcement  (in.) 

height  of  front  wall  or  one-half  its  width,  whichever  is  smaller 
(ft) 

strain  energy 
time  (ms) 

time  increment  (ms) 
any  time  (ms) 

time  of  arrival  of  blast  wave  at  points  b,  e,  and  f,  respectively 
(ms) 

(1)  clearing  time  for  reflected  pressures  (ms) 

(2)  container  thickness  of  explosive  charges  (in.) 
rise  time  (ms) 

time  to  reach  maximum  elastic  deflection 
time  at  which  maximum  deflection  occurs  (ms) 
duration  of  positive  phase  of  blast  pressure  (ms) 
duration  of  negative  phase  of  blast  pressure  (ms) 
fictitious  positive  phase  pressure  duration  (ms) 
fictitious  negative  phase  pressure  duration  (ms) 
fictitious  reflected  pressure  duration  (ms) 
time  at  which  ultimate  deflection  occurs  (ms) 
time  to  reach  yield  (ms) 
time  of  arrival  of  blast  wave  (ms) 


t 

T 

T 

T 

T 

T 

T 

T 


1 


c 

c 

N 

r 

s 

s 


u 


U 


V 

va 

vb 


(avg. ) 


time  at  which  partial  failure  occurs  (ms) 

duration  of  equivalent  triangular  loading  function  (ms) 

thickness  of  concrete  section  (in.) 

scaled  thickness  of  concrete  section  (ft/lb^) 

effective  natural  period  of  vibration  (ms) 

rise  time  (ms) 

thickness  of  sand  fill  (in.) 

scaled  thickness  of  sand  fill  (ft/lb^) 

particle  velocity  (ft/ms) 

ultimate  flexural  or  anchorage  bond  stress  (psi) 
shock  front  velocity  (ft/ms) 
velocity  (in. /ms) 

instantaneous  velocity  at  any  time  (in. /ms) 

boundary  velocity  for  primary  fragments  (ft/sec) 

ultimate  shear  stress  permitted  on  an  unreinforced  web  (psi) 

maximum  post- failure  fragment  velocity  (in. /ms) 

average  post- failure  fragment  velocity  (in. /ms) 

velocity  at  incipient  failure  deflection  (in. /ms) 

initial  velocity  of  primary  fragment  (ft/sec) 

residual  velocity  of  primary  fragment  after  perforation  (ft/sec) 

striking  velocity  of  primary  fragment  (ft/sec) 

ultimate  shear  stress  (psi) 

ultimate  shear  stress  at  distance  de  from  the  horizontal  support 
(psi)* 


See  note  at  end  of  symbols. 
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vuy  ultimate  shear  stress  at  distance  de  from  the  vertical  support 

(spl)* 

V  volume  of  partially  vented  chamber  (ft3) 

Vd  ultimate  direct  shear  capacity  of  the  concrete  of  width  b  (lbs) 

VdH  shear  at  distance  dg  from  the  vertical  support  on  a  unit  width 

(lbs. /In.)* 

VdV  shear  at  distance  de  from  the  horizontal  support  on  a  unit  width 

(Ibs/in. )* 

VQ  volume  of  structure  (ft3) 

V$  shear  at  the  support  on  a  unit  width  (lbs/in)* 

V$H  shear  at  the  vertical  support  on  a  unit  width  (Ibs/in. )* 

Y $y  shear  at  the  horizontal  support  on  a  unit  width  (Ibs/in.)* 

Vy  total  shear  on  a  width  b  (lbs) 

w  weight  density  of  concrete  (lbs/ft3) 

ws  weight  density  of  sand  { lbs/ft3) 

W  charge  weight  (lbs) 

Wc  total  weight  of  explosive  containers  (lbs) 

Wf  weight  of  primary  fragment  (oz) 

Wc0  total  weight  of  steel  core  (lbs) 

wcl*wc2  total  we19ht  of  P1ates  1  and  2,  respectively  (lbs) 

W$  width  of  structure  (ft) 

WD  work  done 

x  yield  line  location  in  horizontal  direction  (in.)* 

*  See  note  at  end  of  symbols. 
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a  (2)  angle  of  incidence  of  the  pressure  front  (deg) 

3  (1)  coefficient  for  determining  elastic  and  elasto-plastic 

resistances 

(2)  particular  support  rotation  angle  (deg) 

Y  coefficient  for  determining  elastic  and  elasto-plastic  deflections 

\  increase  in  support  rotation  angle  after  partial  failure  (deg) 

9  support  rotation  angle  (deg) 

9  angular  acceleration  (rad/ms  ) 

9_av  maximum  support  rotation  angle  (deg) 

9h  horizontal  rotation  angle  (deg)* 

9V  vertical  rotation  angle  (deg)* 

2o  effective  perimeter  of  reinforcing  bars  (in.) 

2m  summation  of  moments  (in. -lbs) 

2Mn  sum  of  the  ultimate  unit  resisting  moments  acting  along  the 

negative  yield  lines  (in. -lbs) 

5mp  sum  of  the  ultimate  unit  resisting  moments  acting  along  the 

positive  yield  lines  (in. -lbs) 

U  ductility  factor 

v  Poisson's  ratio 

<D  (1)  capacity  reduction  factor 

(2)  bar  diameter  (in.) 

<Pr  assumed  shape  function  for  concentrated  loads 


See  note  at  end  of  symbols. 
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<p(x)  assumed  shape  function  for  distributed  loads 


free  edge 

simple  support 


fixed  support 


xxxxxXxm  either  fixed,  restrained,  or  simple  support 


*  Note.  This  symbol  was  developed  for  two-way  elements  which  are  used  as 
walls.  When  roof  slabs  or  other  horizontal  elements  are  under  consideration, 
this  symbol  will  also  be  applicable  if  the  element  is  treated  as  being  rotated 
into  a  vertical  position. 
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